Computational Algebra Assignment3 Michaelgboneh
Computational Algebra Assignment3 Michaelgboneh
Question 1
Given the following reaction system
1k
A −→X
2 k
2X + Y −→3X
3 k
B + X −→Y +D
4k
X −→E
A (-) is when a species is used up while a (+) is used when a species is gained.The differential equation
for the mass action is
d[X]
=k1 [A] + k2 [X]2 [Y ] − k3 [B][X] − k4 [X]
dt
d[Y ]
= − k2 [X]2 [Y ] + k3 [B][X]
dt
1
Question 2
Explicit definition of the propensity function.
Given a reaction
k
s1 x1 + s2 x2 + · · · + sn xn −
→Y
where
n
k Y
cj = sj !,
(nA V )n−1 j=1
k1
C1 =
(nA V )1−1
=k1
n
Y A
=A
i=1
1
a1 = k1 A Propensity function
B.
2 k
2X + Y −→3X
2
k2
C2 = 2!1!
(nA V )3−1
2k2
=
(nA V )2
n
Y X Y X(X − 1)Y
=
i=1
2 1 2
k2 (X)(X − 1)Y
= Propensity function
(nA V )2
C.
3 k
B + X −→Y +D
k3
C3 = 1!1!
(nA V )2−1
k3
=
nA V
n
Y B X
= BX
i=1
1 1
k3 BX
a3 = Propensity function
nA V
4.
4k
X −→E
k4
C4 =
(nA V )1−1
=k4
3
n
Y X
=X
i=1
1
a4 = k4 X Propensity function
4
Question 3
Probability that no reaction fires over time interval [s,t]
where
The sum of all propensities in state x(s) is given by a0 (x(s)).
The probability that a reaction will fire on the time interval (t, t + dτ ) is given byP (τ, j).
The probability that no reaction occurs on the interval (s, t)is given by P0 (τ )
The probability that the reaction will occur on the interval (t, t + dτ ) is given by aj .
This gives
P0 (τ + dτ ) − P0 (τ ) = − P0 (τ )a0
P0 (τ + dτ ) − P0 (τ )
= − P0 (τ )a0
dτ
dP0 (τ ) = − P0 (τ )a0 dτ
dP0 (τ )
= − a0 dτ
P0 (τ )
lnP0 (τ ) =a0 τ
P0 (τ ) =ea0 τ , where τ =t−s
5
Question 4
Proportional relationship for the likelihood of a path in the complete data case:
( M ) ( Z T )
Y
L(c, x) ∝ crj i exp −ci gi (x(t))dt
i=1 0
Solution
From slide 14 of week three and day three the lecture notes,we have the following:
( Z )
Y r T
L(x) = apj (x(tj−1 )) exp − a0 (x(t))dt (1)
j=1 0
Then ( Z )
Yr M
T X
L(c; x) = apj (x(tj−1 ), cj ) exp − ai (x, ci )dt (3)
j=1
0 i=1
Let ai (x, ci ) = ci gi (x), i = 1, ..., M Substituting this into equation (3) and simplifying, we have that
( )
Yr M Z T
X
L(c; x) = cpj gpj (x(tj−1 ), cj ) exp − ci gi (x, ci )dt (3)
j=1
0 i=1
( M
) ( M Z
)
Y X T
∝ crj i exp − ci gi (x(t))dt
i=1 i=1 0
M
( )
Y Z T
= crj i exp −ci gi (x(t))dt simplifying further we generate the below poduct
i=1 0
M
Y
= Li (ci , x) QED
i=1
Where c = (c1 , c2 , ..., cM ) is the vector of stochastic rate constant and x a full sample path of the model.
The component likelihoods are defined by
( )
Z T
Li (ci , x) = crj i exp −ci gi (x(t))dt i = 1, 2, ..., M
0
6
Question 5
The maximum likelihood estimator for each rate constant.
ri
cbi = Pr
g
j=0 i (x(tj ))(tj+1 − tj )
for i = 1,· · · ,M
M is the number of reactions in the model.
From question 4,
RT
Li (ci , x) = (ci )ri eci 0
gi (x(t))dt
Li (ci , x) =0 or
Z T
ri
− gi (x(t))dt =0
ci 0
Z T
ri
= gi (x(t))dt
ci 0
ri
ci = R T
0
gi (x(t))dt