Does someone know how to do part c and d using the code I have provided? I would

Question

Does someone know how to do part c and d using the code I have provided? I would really appreciate it!

#beta start
betastart=1/37.0
#beta stop
betastop=1/6.0

#beta for synthesis
betasynth = 5*betastop
gamma=np.log(2)/50
tmax=200

def step(n):
r_gene = betastop+ betastart
xi_gene = betastart/r_gene
if (np.random.rand()<xi_gene):
state=1 #on
else:
state=0 #off

#mrNA synthesis
   r=betasynth+gamma*n # rate that "something" will happen
dt=-np.log(np.random.rand())/r # time increment
   xi=betasynth/r
if (np.random.rand()<xi):
dn=1 * state #birth
else:
dn=-1 #death
return dn,dt

#average mRNA molecules
def average(storedata):
means = np.zeros(tmax)
for time in range(tmax):
n = 0
for i in range(2000):
for pair in storedata[i]:
# check number of mRNA for each time instance
if pair[0] >= time:
n = n + pair[1]
break
means[time] = n
return means

storedata = []
for i in range (2000):  
nc=0 #molecules
tc=0 #time
pairs = []
while (tc<tmax):
dn,dt=step(nc)
tc+=dt
nc+=dn
pairs.append([tc,nc]) #stores ti,ni
storedata.append(pairs)
means = average(storedata)
t=np.linspace(0, tmax, num=tmax)
plt.plot(t,means)
plt.show()

#numpy.var()

As an application of the birth-death process studied in a previous assignment, consider the synthesis of mRNA with rate Bs, and its clearance from the cell with rate ne. This picture, however, is too simple: the gene that controls the mRNA production makes spontaneous transitions between active ("on") and inactive ("off") states at mean rates Bstart and Bstop Only the active state can be transcribcd, lcading to bursts of m production intersperscd with quiet periods. In effcct this means that whenever the gone is in the "off" state, tho synthesis ratc Bs must be set to zero. The situation is visualized below stg t tivc Nene clear Synthesis O (a) Modify your computer program to implement the transcriptional bursting proccss. Ini- tially the genc is "off" and there are no mRNA molecules in the system. Set ne ln(2)/ (50 min) Bstart 1 37 min Bstop 1/(6 min), BS 3 5Bstop (b) Run the simulation 2000 times for 200 min. Instead of plotting all 2000 trajectories, plot the average number of mRNA moleculcs (m) over all 2000 runs as a function of time. (Hint Storc thc data for cach run in pairs of (ti, ni), where i is a step index. For cvery run, and for every valuc of t bctween 0 and 200 min, find the value of i for which ti t for the first time. The valuc of n corresponding to time t is then ni (c) Tabulate the number of trials in which, for a given timet, the number of mRNA molecules is still zero. Convert this result into a graph of the probability of observing zero mRNA molecules at time t. hat would you have expected for a simple birth-death process (without bursts), and how does the result differ? (d) Compute the ratio "Variance (ngnal)/Mean(nfinal)" for all final valus of n. What do you expect for a simplc Poisson process?

Explanation / Answer

The administration of albuterol, arformoterol, budesonide, cromolyn,    formoterol, ipratropium, levalbuterol, or metaproterenol for the management of obstructive pulmonary disease (ICD-10; J41.0-J70.9); or

b) The administration of dornase alpha to a member with cystic fibrosis (ICD10; E84.0); or

c) The administration of tobramycin to a member with cystic fibrosis or bronchiectasis (ICD-10; E84.0), (ICD-10; J47.9), (ICD-10; J47.1), (ICD-10; Q33.4) or (ICD-10; A15.0-A15.9); or

d) The administration of pentamidine to a members with HIV, (ICD-10; B20.), pneumocystosis (ICD-10; B59), or complications of organ transplants (ICD10; T86.890; T86.891; T86.899); or

e) The administration of acetylcysteine for persistent thick or tenacious pulmonary secretions (ICD-10; J12.0; J70.9); (ICD 10; R09.3).

3. A large volume nebulizer, related compressor, and water or saline are covered when it is medically necessary to deliver humidity to a member with thick, tenacious secretions who has cystic fibrosis, (ICD 10; R09.3), bronchiectasis (ICD-10; J47.9), (ICD-10; J47.1), (ICD-10; A15.0) or (ICD-10; Q33.4), a tracheostomy (ICD-10; Z93.0 or V55.0), or a tracheobronchial stent (ICD 10; J39.8 and J98.09).

4. An E0565 or E0572 compressor and filtered nebulizer (A7006) are also covered when it is reasonable and necessary to administer pentamidine to members with HIV (ICD-10; B20), pneumocystosis (ICD 10; B59); or complications of organ transplants (ICD 10; T86.90; T86.91; T86.92; T86.99) and, (ICD 10; T86.890; T86.89; T86.899).

Related Questions

Navigate

• Browse (All)
• Browse D
• Subjects

---

---

Integrity-first tutoring: explanations and feedback only — we do not complete graded work. Learn more.