CDOoDocuments.StdDocumentDescDocuments.DocumentDescContainers.ViewDescViews.ViewDescStores.StoreDescUTDocuments.ModelDescContainers.ModelDescModels.ModelDescStores.ElemDescS STextViews.StdViewDescTextViews.ViewDescESTextModels.StdModelDescTextModels.ModelDescOQGQTextModels.AttributesDesc'* 0Arial% e$Courier New1VQT|.StdFolds.FoldDesc  b@ $Courier NewPSTransition probabilities[r2|. model{ #Screen detected polips / cancers #========================= for(l in 1:4){ for(k in 1:num.screens){ pol.l[k,l]~dbin(p.pol.l[k,l], n.scr[k,l]) pol.h[k,l]~dbin(p.pol.h[k,l], n.scr[k,l]) can.a[k,l]~dbin(p.can.a[k,l], n.scr[k,l]) can.b[k,l]~dbin(p.can.b[k,l], n.scr[k,l]) can.d[k,l]~dbin(p.can.d[k,l], n.scr[k,l]) p.pol.l[k,l]<-s[2]*rho[2,k,l] p.pol.h[k,l]<-s[3]*rho[3,k,l] p.can.a[k,l]<-s[4]*rho[4,k,l] p.can.b[k,l]<-s[6]*rho[6,k,l] p.can.d[k,l]<-s[8]*rho[8,k,l] } } #Interval cancers #============ for(l in 1:4){ #Rate parameters la[l]<-laa[l]+lba[l]+lda[l] lb[l]<-lba[l]+lda[l] ld[l]<-lda[l] for(t in 1:2){ can.int.a[t,l]~dpois(lambda.a[t,l]) can.int.b[t,l]~dpois(lambda.b[t,l]) can.int.d[t,l]~dpois(lambda.d[t,l]) lambda.a[t,l] <-py[t,l]*la[l]*(1-exp(-(l4[l]+l5[l])*(t-0.5))) +sum(can.a[,l])*(1-s[4])/s[4] * (exp(-(l4[l]+l5[l])*(t-1)) -exp(-(l4[l]+l5[l])*t) ) lambda.b[t,l] <-py[t,l]*lb[l]*(1-exp(-(l6[l]+l7[l])*(t-0.5))) +sum(can.b[,l])*(1-s[6])/s[6] * (exp(-(l6[l]+l7[l])*(t-1)) -exp(-(l6[l]+l7[l])*t) ) lambda.d[t,l] <-py[t,l]*ld[l]*(1-exp(-(l8[l])*(t-0.5))) +sum(can.d[,l])*(1-s[8])/s[8] * (exp(-(l8[l])*(t-1)) -exp(-(l8[l])*t) ) } } #Recursive equations for observed and true state probabilities #=============================================== #Prevalance and incidence screen probability of being observed in well state #theta[k,l]=theta[screening round, age agroup]=P(observed W at round k for age group l) for(l in 1:4){ for(k in 1:num.screens){ theta[k,l]<-(1-s[2])*rho[2,k,l]+(1-s[3])*rho[3,k,l]+(1-s[4])*rho[4,k,l]+(1-s[6])*rho[6,k,l]+(1-s[8])*rho[8,k,l] } } #Prevalance screen probability of being in true state j (j=2,3,4,6,8) #rho[j,k,l]=rho[state, screening round, age group]=P(true state = j at round k for age group l) for(l in 1:4){ p.can[l]<-pl[5,l]+pl[7,l]+pl[9,l] rho[2,1,l]<-pl[2,l]/(1-p.can[l]) rho[3,1,l]<-pl[3,l]/(1-p.can[l]) rho[4,1,l]<-pl[4,l]/(1-p.can[l]) rho[6,1,l]<-pl[6,l]/(1-p.can[l]) rho[8,1,l]<-pl[8,l]/(1-p.can[l]) rho[1,1,l]<-0 rho[5,1,l]<-0 rho[7,1,l]<-0 } #Incidence screen probability of being in true state j (j=2,3,4,5,6) for(l in 1:4){ for(k in 2:num.screens){ rho[2,k,l]<-(1-s[2])*rho[2,k-1,l]*p22[l,2]/theta[k-1,l] rho[3,k,l]<-((1-s[2])*rho[2,k-1,l]*p23[l,2]+(1-s[3])*rho[3,k-1,l]*p33[l,2])/theta[k-1,l] rho[4,k,l]<-((1-s[2])*rho[2,k-1,l]*p24[l,2]+(1-s[3])*rho[3,k-1,l]*p34[l,2]+(1-s[4])*rho[4,k-1,l]*p44[l,2])/theta[k-1,l] rho[6,k,l]<-((1-s[2])*rho[2,k-1,l]*p26[l,2]+(1-s[3])*rho[3,k-1,l]*p36[l,2]+(1-s[4])*rho[4,k-1,l]*p46[l,2]+(1-s[6])*rho[6,k-1,l]*p66[l,2])/theta[k-1,l] rho[8,k,l]<-((1-s[2])*rho[2,k-1,l]*p28[l,2]+(1-s[3])*rho[3,k-1,l]*p38[l,2]+(1-s[4])*rho[4,k-1,l]*p48[l,2]+(1-s[6])*rho[6,k-1,l]*p68[l,2]+(1-s[8])*rho[8,k-1,l]*p88[l,2])/theta[k-1,l] rho[1,k,l]<-0 rho[5,k,l]<-0 rho[7,k,l]<-0 } } #Sensitivities s[1]<-0 #P(observed in non-well state when well) s[2]~dunif(0,1) #P(observed in state L when true state is L) s[3]<-s[2] #P(observed in state H when true state is H) #s[3]~dunif(0,1) s[4]~dunif(0,1) #P(observed in state A when true state is A) s[6]<-s[4] #P(observed in state B when true state is B) s[8]<-s[4] #P(observed in state D when true state is D) s[5]<-0 s[7]<-0 #ONS data #======== for(l in 1:4){ ons.r[l]~dbin(clincaninc[l], ons.n[l]) clincaninc[l]<- (pl[1,l]*(p15[l,1]+p17[l,1]+p19[l,1]) + pl[2,l]*(p25[l,1]+p27[l,1]+p29[l,1]) + pl[3,l]*(p35[l,1]+p37[l,1]+p39[l,1]) + pl[4,l]*(p45[l,1]+p47[l,1]+p49[l,1]) + pl[6,l]*(p67[l,1]+p69[l,1]) + pl[8,l]*(p89[l,1])) } #Wessex audit on stage distribution #================================== for(l in 1:4){ #age for(stage in 1:3){ #stage wessex.r[stage,l]~dbin(p.stage[stage,l],wessex.n[l]) } p.stage[1,l]<-pl[5,l]/p.can[l] p.stage[2,l]<-pl[7,l]/p.can[l] p.stage[3,l]<-pl[9,l]/p.can[l] } #Polyp necrospy study #==================== for(l in 1:4){ #age nec.r[l]~dbin(p.pol23[l],nec.n[l]) p.pol23[l] <- pl[2,l]+pl[3,l] #prob of low or high risk polyp by age l } #Transition rates #============ #Interpolate the transition rates across age groups for (l in 2:3){ l1[l]<- ((4-l)*l1[1] + (l-1)*l1[4]) / 3 l2[l]<- ((4-l)*l2[1] + (l-1)*l2[4]) / 3 l3[l]<- ((4-l)*l3[1] + (l-1)*l3[4]) / 3 l4[l]<- ((4-l)*l4[1] + (l-1)*l4[4]) / 3 l5[l]<- ((4-l)*l5[1] + (l-1)*l5[4]) / 3 l6[l]<- ((4-l)*l6[1] + (l-1)*l6[4]) / 3 l7[l]<- ((4-l)*l7[1] + (l-1)*l7[4]) / 3 l8[l]<- ((4-l)*l8[1] + (l-1)*l8[4]) / 3 } l1[1]~dunif(0.001,5) l2[1]~dunif(0.001,5) l3[1]~dunif(0.001,5) l4[1]~dunif(0.001,5) l5[1]~dunif(0.001,5) l6[1]~dunif(0.001,5) l7[1]~dunif(0.001,5) l8[1]~dunif(0.001,5) l1[4]~dunif(0.001,5) l2[4]~dunif(0.001,5) l3[4]~dunif(0.001,5) l4[4]~dunif(0.001,5) l5[4]~dunif(0.001,5) l6[4]~dunif(0.001,5) l7[4]~dunif(0.001,5) l8[4]~dunif(0.001,5) #transition probabilities #================== #pij[l,t]=pij[age, time]=probability of being in state j from state i after time t if age l #set up matrix for these transition probabilities #P[i, j, l ] = probability of being in state j friom state i after time 5 years (time[3]) for age l. for(l in 2:4){ P[1,1,l]<-p11[l,3] P[1,2,l]<-p12[l,3] P[1,3,l]<-p13[l,3] P[1,4,l]<-p14[l,3] P[1,5,l]<-p15[l,3] P[1,6,l]<-p16[l,3] P[1,7,l]<-p17[l,3] P[1,8,l]<-p18[l,3] P[1,9,l]<-p19[l,3] P[2,2,l]<-p22[l,3] P[2,3,l]<-p23[l,3] P[2,4,l]<-p24[l,3] P[2,5,l]<-p25[l,3] P[2,6,l]<-p26[l,3] P[2,7,l]<-p27[l,3] P[2,8,l]<-p28[l,3] P[2,9,l]<-p29[l,3] P[3,3,l]<-p33[l,3] P[3,4,l]<-p34[l,3] P[3,5,l]<-p35[l,3] P[3,6,l]<-p36[l,3] P[3,7,l]<-p37[l,3] P[3,8,l]<-p38[l,3] P[3,9,l]<-p39[l,3] P[4,4,l]<-p44[l,3] P[4,5,l]<-p45[l,3] P[4,6,l]<-p46[l,3] P[4,7,l]<-p47[l,3] P[4,8,l]<-p48[l,3] P[4,9,l]<-p49[l,3] P[5,5,l]<-1 P[5,6,l]<-0 P[5,7,l]<-0 P[5,8,l]<-0 P[5,9,l]<-0 P[6,6,l]<-p66[l,3] P[6,7,l]<-p67[l,3] P[6,8,l]<-p68[l,3] P[6,9,l]<-p69[l,3] P[7,7,l]<-1 P[7,8,l]<-0 P[7,9,l]<-0 P[8,8,l]<-p88[l,3] P[8,9,l]<-p89[l,3] P[9,9,l]<-1 for(i in 2:9){ for(j in 1:(i-1)){ P[i,j,l]<-0 } } } #pl[j,l] = probabilities for being in stage j from stage 1 for age group l years at risk #l = 22.5, 27.5, 32.5, 37.5, cancers start at age 30. #uses the age specific transition rates. pl[1,1]<-p11[1,4] pl[2,1]<-p12[1,4] pl[3,1]<-p13[1,4] pl[4,1]<-p14[1,4] pl[5,1]<-p15[1,4] pl[6,1]<-p16[1,4] pl[7,1]<-p17[1,4] pl[8,1]<-p18[1,4] pl[9,1]<-p19[1,4] for(l in 2:4){ for(j in 1:9){ pl[j,l]<-pl[1,l-1]*P[1,j,l]+pl[2,l-1]*P[2,j,l]+pl[3,l-1]*P[3,j,l]+pl[4,l-1]*P[4,j,l]+pl[5,l-1]*P[5,j,l]+pl[6,l-1]*P[6,j,l]+pl[7,l-1]*P[7,j,l]+pl[8,l-1]*P[8,j,l]+pl[9,l-1]*P[9,j,l] } } #the relationship between lambda and the transition probabilities for(l in 1:4){ #l = index for age for(t in 1:4){ #t = index for time p11[l,t]<-(exp(-l1[l]*time[t])) p12[l,t]<-((((exp(-l1[l]*time[t])-exp(-l2[l]*time[t])))*l1[l])/(-l1[l]+l2[l])) p13[l,t]<-((l1[l]*l2[l]*((exp(-l3[l]*time[t])*((l1[l]-l2[l]))+exp(-l1[l]*time[t])*((l2[l]-l3[l]))+exp(-l2[l]*time[t])*((-l1[l]+l3[l])))))/(((l1[l]-l2[l]))*((l1[l]-l3[l]))*((l2[l]-l3[l])))) p14a[l,t]<-exp((-((l4[l]+l5[l])))*time[t])*((l1[l]-l2[l]))*((l1[l]-l3[l]))*((l2[l]-l3[l]))-exp(-l3[l]*time[t])*((-l1[l]+l2[l]))*((l2[l]-l4[l]-l5[l]))*((-l1[l]+l4[l]+l5[l])) p14b[l,t]<-exp(-l2[l]*time[t])*((-l1[l]+l3[l]))*((l3[l]-l4[l]-l5[l]))*((-l1[l]+l4[l]+l5[l]))-exp(-l1[l]*time[t])*((-l2[l]+l3[l]))*((l3[l]-l4[l]-l5[l]))*((-l2[l]+l4[l]+l5[l])) p14[l,t]<-(((l1[l]*l2[l]*l3[l]*((p14a[l,t]+p14b[l,t]))))/((((l1[l]-l2[l]))*((l1[l]-l3[l]))*((l2[l]-l3[l]))*((l1[l]-l4[l]-l5[l]))*((l2[l]-l4[l]-l5[l]))*((l3[l]-l4[l]-l5[l]))))) p15a[l,t]<-exp(-l3[l]*time[t])*l1[l]*((l1[l]-l2[l]))*l2[l]*((l1[l]-l4[l]-l5[l]))*((l2[l]-l4[l]-l5[l]))*((l4[l]+l5[l])) p15b[l,t]<- -exp(-l2[l]*time[t])*l1[l]*((l1[l]-l3[l]))*l3[l]*((l1[l]-l4[l]-l5[l]))*((l3[l]-l4[l]-l5[l]))*((l4[l]+l5[l]))+exp(-l1[l]*time[t])*l2[l]*((l2[l]-l3[l]))*l3[l]*((l2[l]-l4[l]-l5[l]))*((l3[l]-l4[l]-l5[l]))*((l4[l]+l5[l]))-((-l1[l]+l2[l]))*((-l1[l]+l3[l]))*((-l2[l]+l3[l]))*((l3[l]-l4[l]-l5[l]))*((-l1[l]+l4[l]+l5[l]))*((-l2[l]+l4[l]+l5[l])) p15[l,t]<-(((l4[l]*(((-exp((-((l4[l]+l5[l])))*time[t]))*l1[l]*((l1[l]-l2[l]))*l2[l]*((l1[l]-l3[l]))*((l2[l]-l3[l]))*l3[l]+p15a[l,t]+p15b[l,t]))))/((((l1[l]-l2[l]))*((l1[l]-l3[l]))*((l2[l]-l3[l]))*((l1[l]-l4[l]-l5[l]))*((l2[l]-l4[l]-l5[l]))*((l3[l]-l4[l]-l5[l]))*((l4[l]+l5[l]))))) p16a[l,t]<-(-exp((-((l6[l]+l7[l])))*time[t]))*((-l1[l]+l2[l]))*((-l1[l]+l3[l]))*((-l2[l]+l3[l]))*((l3[l]-l4[l]-l5[l]))*((-l1[l]+l4[l]+l5[l]))*((-l2[l]+l4[l]+l5[l])) p16b[l,t]<-exp((-((l4[l]+l5[l])))*time[t])*((-l1[l]+l2[l]))*((-l1[l]+l3[l]))*((-l2[l]+l3[l]))*((l3[l]-l6[l]-l7[l]))*((-l1[l]+l6[l]+l7[l]))*((-l2[l]+l6[l]+l7[l])) p16c[l,t]<- -exp(-l3[l]*time[t])*((-l1[l]+l2[l]))*((-l1[l]+l4[l]+l5[l]))*((-l2[l]+l4[l]+l5[l]))*((l4[l]+l5[l]-l6[l]-l7[l]))*((-l1[l]+l6[l]+l7[l]))*((-l2[l]+l6[l]+l7[l])) p16d[l,t]<-exp(-l2[l]*time[t])*((-l1[l]+l3[l]))*((-l1[l]+l4[l]+l5[l]))*((-l3[l]+l4[l]+l5[l]))*((l4[l]+l5[l]-l6[l]-l7[l]))*((-l1[l]+l6[l]+l7[l]))*((-l3[l]+l6[l]+l7[l])) p16e[l,t]<- -exp(-l1[l]*time[t])*((-l2[l]+l3[l]))*((-l2[l]+l4[l]+l5[l]))*((-l3[l]+l4[l]+l5[l]))*((l4[l]+l5[l]-l6[l]-l7[l]))*((-l2[l]+l6[l]+l7[l]))*((-l3[l]+l6[l]+l7[l])) p16[l,t]<-(((l1[l]*l2[l]*l3[l]*l5[l]*((p16a[l,t]+p16b[l,t]+ p16c[l,t]+ p16d[l,t] + p16e[l,t]))))/((((l1[l]-l2[l]))*((l1[l]-l3[l]))*((l2[l]-l3[l]))*((l1[l]-l4[l]-l5[l]))*((l2[l]-l4[l]-l5[l]))*((l3[l]-l4[l]-l5[l]))*((l1[l]-l6[l]-l7[l]))*((l2[l]-l6[l]-l7[l]))*((l3[l]-l6[l]-l7[l]))*((l4[l]+l5[l]-l6[l]-l7[l]))))) p17a[l,t]<-(-exp((-((l6[l]+l7[l])))*time[t]))*l1[l]*((l1[l]-l2[l]))*l2[l]*((l1[l]-l3[l]))*((l2[l]-l3[l]))*l3[l]*((l1[l]-l4[l]-l5[l]))*((l2[l]-l4[l]-l5[l]))*((l3[l]-l4[l]-l5[l]))*((l4[l]+l5[l])) p17b[l,t]<-exp((-((l4[l]+l5[l])))*time[t])*l1[l]*((l1[l]-l2[l]))*l2[l]*((l1[l]-l3[l]))*((l2[l]-l3[l]))*l3[l]*((l1[l]-l6[l]-l7[l]))*((l2[l]-l6[l]-l7[l]))*((l3[l]-l6[l]-l7[l]))*((l6[l]+l7[l])) p17c[l,t]<-exp(-l2[l]*time[t])*l1[l]*((l1[l]-l3[l]))*l3[l]*((l1[l]-l4[l]-l5[l]))*((l3[l]-l4[l]-l5[l]))*((l4[l]+l5[l]))*((l1[l]-l6[l]-l7[l]))*((l3[l]-l6[l]-l7[l]))*((l4[l]+l5[l]-l6[l]-l7[l]))*((l6[l]+l7[l])) p17d[l,t]<- -((-l1[l]+l2[l]))*((-l1[l]+l3[l]))*((-l2[l]+l3[l]))*((-l1[l]+l4[l]+l5[l]))*((-l2[l]+l4[l]+l5[l]))*((-l3[l]+l4[l]+l5[l]))*((l4[l]+l5[l]-l6[l]-l7[l]))*((-l1[l]+l6[l]+l7[l]))*((-l2[l]+l6[l]+l7[l]))*((-l3[l]+l6[l]+l7[l])) p17e[l,t]<-exp(-l3[l]*time[t])*l1[l]*((l1[l]-l2[l]))*l2[l]*((l1[l]-l4[l]-l5[l]))*((l2[l]-l4[l]-l5[l]))*((l4[l]+l5[l]))*((l1[l]-l6[l]-l7[l]))*((l2[l]-l6[l]-l7[l]))*((l6[l]+l7[l]))*((-l4[l]-l5[l]+l6[l]+l7[l])) p17f[l,t]<-exp(-l1[l]*time[t])*l2[l]*((l2[l]-l3[l]))*l3[l]*((l2[l]-l4[l]-l5[l]))*((l3[l]-l4[l]-l5[l]))*((l4[l]+l5[l]))*((l2[l]-l6[l]-l7[l]))*((l3[l]-l6[l]-l7[l]))*((l6[l]+l7[l]))*((-l4[l]-l5[l]+l6[l]+l7[l])) p17[l,t]<-(((l5[l]*l6[l]*(( p17a[l,t]+ p17b[l,t]+ p17c[l,t] + p17d[l,t] + p17e[l,t]+p17f[l,t]))))/((((l1[l]-l2[l]))*((l1[l]-l3[l]))*((l2[l]-l3[l]))*((l1[l]-l4[l]-l5[l]))*((l2[l]-l4[l]-l5[l]))*((l3[l]-l4[l]-l5[l]))*((l4[l]+l5[l]))*((l1[l]-l6[l]-l7[l]))*((l2[l]-l6[l]-l7[l]))*((l3[l]-l6[l]-l7[l]))*((l4[l]+l5[l]-l6[l]-l7[l]))*((l6[l]+l7[l]))))) p18a[l,t]<-(-exp(-l8[l]*time[t]))*(l2[l]-l1[l])*(l3[l]-l1[l])*(l3[l]-l2[l])*(((-l1[l])+l4[l]+l5[l]))*(((-l2[l])+l4[l]+l5[l]))*(((-l3[l])+l4[l]+l5[l]))*((l4[l]+l5[l]-l6[l]-l7[l]))*(((-l1[l])+l6[l]+l7[l]))*(((-l2[l])+l6[l]+l7[l]))*(((-l3[l])+l6[l]+l7[l])) p18b[l,t]<-exp((-l3[l])*time[t])*(((-l1[l])+l2[l]))*(((-l1[l])+l4[l]+l5[l]))*(((-l2[l])+l4[l]+l5[l]))*((l4[l]+l5[l]-l6[l]-l7[l]))*(((-l1[l])+l6[l]+l7[l]))*(((-l2[l])+l6[l]+l7[l]))*((l4[l]+l5[l]-l8[l]))*((l6[l]+l7[l]-l8[l]))*(((-l1[l])+l8[l]))*(((-l2[l])+l8[l])) p18c[l,t]<- -exp((-l2[l])*time[t])*(((-l1[l])+l3[l]))*(((-l1[l])+l4[l]+l5[l]))*(((-l3[l])+l4[l]+l5[l]))*((l4[l]+l5[l]-l6[l]-l7[l]))*(((-l1[l])+l6[l]+l7[l]))*(((-l3[l])+l6[l]+l7[l]))*((l4[l]+l5[l]-l8[l]))*((l6[l]+l7[l]-l8[l]))*(((-l1[l])+l8[l]))*(((-l3[l])+l8[l])) p18d[l,t]<-exp((-l1[l])*time[t])*(((-l2[l])+l3[l]))*(((-l2[l])+l4[l]+l5[l]))*(((-l3[l])+l4[l]+l5[l]))*((l4[l]+l5[l]-l6[l]-l7[l]))*(((-l2[l])+l6[l]+l7[l]))*(((-l3[l])+l6[l]+l7[l]))*((l4[l]+l5[l]-l8[l]))*((l6[l]+l7[l]-l8[l]))*(((-l2[l])+l8[l]))*(((-l3[l])+l8[l])) p18e[l,t]<- -exp((-((l6[l]+l7[l])))*time[t])*(((-l1[l])+l2[l]))*(((-l1[l])+l3[l]))*(((-l2[l])+l3[l]))*(((-l1[l])+l4[l]+l5[l]))*(((-l2[l])+l4[l]+l5[l]))*(((-l3[l])+l4[l]+l5[l]))*(((-l1[l])+l8[l]))*(((-l2[l])+l8[l]))*(((-l3[l])+l8[l]))*(((-l4[l])-l5[l]+l8[l])) p18f[l,t]<-exp((-((l4[l]+l5[l])))*time[t])*(((-l1[l])+l2[l]))*(((-l1[l])+l3[l]))*(((-l2[l])+l3[l]))*(((-l1[l])+l6[l]+l7[l]))*(((-l2[l])+l6[l]+l7[l]))*(((-l3[l])+l6[l]+l7[l]))*(((-l1[l])+l8[l]))*(((-l2[l])+l8[l]))*(((-l3[l])+l8[l]))*(((-l6[l])-l7[l]+l8[l])) p18g[l,t]<-((((l1[l]-l2[l]))*((l1[l]-l3[l]))*((l2[l]-l3[l]))*((l1[l]-l4[l]-l5[l]))*((l2[l]-l4[l]-l5[l]))*((l3[l]-l4[l]-l5[l]))*((l1[l]-l6[l]-l7[l]))*((l2[l]-l6[l]-l7[l]))*((l3[l]-l6[l]-l7[l]))*((l4[l]+l5[l]-l6[l]-l7[l]))*((l1[l]-l8[l]))*((l2[l]-l8[l]))*((l3[l]-l8[l]))*((l4[l]+l5[l]-l8[l]))*((l6[l]+l7[l]-l8[l])))) p18[l,t]<-(((l1[l]*l2[l]*l3[l]*l5[l]*l7[l]*((p18a[l,t]+p18b[l,t]+p18c[l,t]+p18d[l,t]+p18e[l,t]+p18f[l,t]))))/p18g[l,t]) p22[l,t]<-(exp((-l2[l])*time[t])) p23[l,t]<-((((exp((-l2[l])*time[t])-exp((-l3[l])*time[t])))*l2[l])/((-l2[l])+l3[l])) p24[l,t]<-((l2[l]*l3[l]*((exp((-((l4[l]+l5[l])))*time[t])*((l2[l]-l3[l]))+exp((-l2[l])*time[t])*((l3[l]-l4[l]-l5[l]))+exp((-l3[l])*time[t])*(((-l2[l])+l4[l]+l5[l])))))/(((l2[l]-l3[l]))*((l2[l]-l4[l]-l5[l]))*((l3[l]-l4[l]-l5[l])))) p25[l,t]<-(((l4[l]*(((-exp((-((l4[l]+l5[l])))*time[t]))*l2[l]*((l2[l]-l3[l]))*l3[l]+exp((-l3[l])*time[t])*l2[l]*((l2[l]-l4[l]-l5[l]))*((l4[l]+l5[l]))-exp((-l2[l])*time[t])*l3[l]*((l3[l]-l4[l]-l5[l]))*((l4[l]+l5[l]))+(((-l2[l])+l3[l]))*((l3[l]-l4[l]-l5[l]))*(((-l2[l])+l4[l]+l5[l]))))))/((((l2[l]-l3[l]))*((l2[l]-l4[l]-l5[l]))*((l3[l]-l4[l]-l5[l]))*((l4[l]+l5[l]))))) p26a[l,t]<-(-exp((-((l6[l]+l7[l])))*time[t]))*(((-l2[l])+l3[l]))*((l3[l]-l4[l]-l5[l]))*(((-l2[l])+l4[l]+l5[l]))+exp((-((l4[l]+l5[l])))*time[t])*(((-l2[l])+l3[l]))*((l3[l]-l6[l]-l7[l]))*(((-l2[l])+l6[l]+l7[l])) p26b[l,t]<- -exp((-l3[l])*time[t])*(((-l2[l])+l4[l]+l5[l]))*((l4[l]+l5[l]-l6[l]-l7[l]))*(((-l2[l])+l6[l]+l7[l]))+exp((-l2[l])*time[t])*(((-l3[l])+l4[l]+l5[l]))*((l4[l]+l5[l]-l6[l]-l7[l]))*(((-l3[l])+l6[l]+l7[l])) p26[l,t]<-((-(((l2[l]*l3[l]*l5[l]*((p26a[l,t]+p26b[l,t]))))/((((l2[l]-l3[l]))*((l2[l]-l4[l]-l5[l]))*((l3[l]-l4[l]-l5[l]))*((l2[l]-l6[l]-l7[l]))*((l3[l]-l6[l]-l7[l]))*((l4[l]+l5[l]-l6[l]-l7[l]))))))) p27a[l,t]<-(-exp((-((l6[l]+l7[l])))*time[t]))*l2[l]*((l2[l]-l3[l]))*l3[l]*((l2[l]-l4[l]-l5[l]))*((l3[l]-l4[l]-l5[l]))*((l4[l]+l5[l]))+exp((-((l4[l]+l5[l])))*time[t])*l2[l]*((l2[l]-l3[l]))*l3[l]*((l2[l]-l6[l]-l7[l]))*((l3[l]-l6[l]-l7[l]))*((l6[l]+l7[l])) p27b[l,t]<-exp((-l2[l])*time[t])*l3[l]*((l3[l]-l4[l]-l5[l]))*((l4[l]+l5[l]))*((l3[l]-l6[l]-l7[l]))*((l4[l]+l5[l]-l6[l]-l7[l]))*((l6[l]+l7[l]))-(((-l2[l])+l3[l]))*(((-l2[l])+l4[l]+l5[l]))*(((-l3[l])+l4[l]+l5[l]))*((l4[l]+l5[l]-l6[l]-l7[l]))*(((-l2[l])+l6[l]+l7[l]))*(((-l3[l])+l6[l]+l7[l])) p27c[l,t]<-exp((-l3[l])*time[t])*l2[l]*((l2[l]-l4[l]-l5[l]))*((l4[l]+l5[l]))*((l2[l]-l6[l]-l7[l]))*((l6[l]+l7[l]))*(((-l4[l])-l5[l]+l6[l]+l7[l])) p27[l,t]<-(((l5[l]*l6[l]*((p27a[l,t]+p27b[l,t]+p27c[l,t]))))/((((l2[l]-l3[l]))*((l2[l]-l4[l]-l5[l]))*((l3[l]-l4[l]-l5[l]))*((l4[l]+l5[l]))*((l2[l]-l6[l]-l7[l]))*((l3[l]-l6[l]-l7[l]))*((l4[l]+l5[l]-l6[l]-l7[l]))*((l6[l]+l7[l]))))) p28a[l,t]<-(-exp((-l8[l])*time[t]))*(((-l2[l])+l3[l]))*(((-l2[l])+l4[l]+l5[l]))*(((-l3[l])+l4[l]+l5[l]))*((l4[l]+l5[l]-l6[l]-l7[l]))*(((-l2[l])+l6[l]+l7[l]))*(((-l3[l])+l6[l]+l7[l]))+exp((-l3[l])*time[t])*(((-l2[l])+l4[l]+l5[l]))*((l4[l]+l5[l]-l6[l]-l7[l]))*(((-l2[l])+l6[l]+l7[l]))*((l4[l]+l5[l]-l8[l]))*((l6[l]+l7[l]-l8[l]))*(((-l2[l])+l8[l])) p28b[l,t]<- -exp((-l2[l])*time[t])*(((-l3[l])+l4[l]+l5[l]))*((l4[l]+l5[l]-l6[l]-l7[l]))*(((-l3[l])+l6[l]+l7[l]))*((l4[l]+l5[l]-l8[l]))*((l6[l]+l7[l]-l8[l]))*(((-l3[l])+l8[l]))-exp((-((l6[l]+l7[l])))*time[t])*(((-l2[l])+l3[l]))*(((-l2[l])+l4[l]+l5[l]))*(((-l3[l])+l4[l]+l5[l]))*(((-l2[l])+l8[l]))*(((-l3[l])+l8[l]))*(((-l4[l])-l5[l]+l8[l])) p28c[l,t]<-exp((-((l4[l]+l5[l])))*time[t])*(((-l2[l])+l3[l]))*(((-l2[l])+l6[l]+l7[l]))*(((-l3[l])+l6[l]+l7[l]))*(((-l2[l])+l8[l]))*(((-l3[l])+l8[l]))*(((-l6[l])-l7[l]+l8[l])) p28[l,t]<-(((l2[l]*l3[l]*l5[l]*l7[l]*((p28a[l,t]+p28b[l,t]+p28c[l,t]))))/((((l2[l]-l3[l]))*((l2[l]-l4[l]-l5[l]))*((l3[l]-l4[l]-l5[l]))*((l2[l]-l6[l]-l7[l]))*((l3[l]-l6[l]-l7[l]))*((l4[l]+l5[l]-l6[l]-l7[l]))*((l2[l]-l8[l]))*((l3[l]-l8[l]))*((l4[l]+l5[l]-l8[l]))*((l6[l]+l7[l]-l8[l]))))) p33[l,t]<-(exp((-l3[l])*time[t])) p34[l,t]<-((((exp((-l3[l])*time[t])-exp((((-l4[l])-l5[l]))*time[t])))*l3[l])/((-l3[l])+l4[l]+l5[l])) p35[l,t]<-((exp((-l3[l])*time[t])*l4[l]*((exp(((l3[l]-l4[l]-l5[l]))*time[t])*l3[l]-l4[l]-l5[l]+exp(l3[l]*time[t])*(((-l3[l])+l4[l]+l5[l])))))/(((l4[l]+l5[l]))*(((-l3[l])+l4[l]+l5[l])))) p36[l,t]<-(((exp((-l3[l])*time[t])*l3[l]*l5[l]*((l4[l]+exp(((l3[l]-l6[l]-l7[l]))*time[t])*((l3[l]-l4[l]-l5[l]))+l5[l]-l6[l]-l7[l]+exp(((l3[l]-l4[l]-l5[l]))*time[t])*(((-l3[l])+l6[l]+l7[l]))))))/((((l3[l]-l4[l]-l5[l]))*((l3[l]-l6[l]-l7[l]))*((l4[l]+l5[l]-l6[l]-l7[l]))))) p37a[l,t]<-(-exp(((l3[l]-l6[l]-l7[l]))*time[t]))*l3[l]*((l3[l]-l4[l]-l5[l]))*((l4[l]+l5[l]))+exp(((l3[l]-l4[l]-l5[l]))*time[t])*l3[l]*((l3[l]-l6[l]-l7[l]))*((l6[l]+l7[l]))-((l4[l]+l5[l]))*((l4[l]+l5[l]-l6[l]-l7[l]))*((l6[l]+l7[l])) p37b[l,t]<-exp(l3[l]*time[t])*(((-l3[l])+l4[l]+l5[l]))*((l4[l]+l5[l]-l6[l]-l7[l]))*(((-l3[l])+l6[l]+l7[l])) p37[l,t]<-((-(((exp((-l3[l])*time[t])*l5[l]*l6[l]*((p37a[l,t]+p37b[l,t]))))/((((l4[l]+l5[l]))*(((-l3[l])+l4[l]+l5[l]))*((l3[l]-l6[l]-l7[l]))*((l4[l]+l5[l]-l6[l]-l7[l]))*((l6[l]+l7[l]))))))) p38a[l,t]<-exp(l3[l]*time[t])*(((-l3[l])+l4[l]+l5[l]))*((l4[l]+l5[l]-l6[l]-l7[l]))*(((-l3[l])+l6[l]+l7[l]))-exp(l8[l]*time[t])*((l4[l]+l5[l]-l6[l]-l7[l]))*((l4[l]+l5[l]-l8[l]))*((l6[l]+l7[l]-l8[l])) p38b[l,t]<-exp(((l3[l]-l6[l]-l7[l]+l8[l]))*time[t])*(((-l3[l])+l4[l]+l5[l]))*(((-l3[l])+l8[l]))*(((-l4[l])-l5[l]+l8[l]))-exp(((l3[l]-l4[l]-l5[l]+l8[l]))*time[t])*(((-l3[l])+l6[l]+l7[l]))*(((-l3[l])+l8[l]))*(((-l6[l])-l7[l]+l8[l])) p38[l,t]<-(((exp((((-l3[l])-l8[l]))*time[t])*l3[l]*l5[l]*l7[l]*((p38a[l,t]+p38b[l,t]))))/((((l3[l]-l4[l]-l5[l]))*((l3[l]-l6[l]-l7[l]))*((l4[l]+l5[l]-l6[l]-l7[l]))*((l3[l]-l8[l]))*((l4[l]+l5[l]-l8[l]))*((l6[l]+l7[l]-l8[l]))))) p44[l,t]<-(exp((((-l4[l])-l5[l]))*time[t])) p45[l,t]<-((l4[l]-exp((((-l4[l])-l5[l]))*time[t])*l4[l])/(l4[l]+l5[l])) p46[l,t]<-(((((-exp((((-l4[l])-l5[l]))*time[t]))+exp((((-l6[l])-l7[l]))*time[t])))*l5[l])/(l4[l]+l5[l]-l6[l]-l7[l])) p47[l,t]<-((-((l5[l]*l6[l]*(((((-1)+exp((((-l6[l])-l7[l]))*time[t])))*l4[l]+(((-1)+exp((((-l6[l])-l7[l]))*time[t])))*l5[l]+l6[l]+l7[l]-exp((((-l4[l])-l5[l]))*time[t])*((l6[l]+l7[l])))))/(((l4[l]+l5[l]))*((l4[l]+l5[l]-l6[l]-l7[l]))*((l6[l]+l7[l])))))) p48[l,t]<-(((l5[l]*l7[l]*((exp((-l8[l])*time[t])*((l4[l]+l5[l]-l6[l]-l7[l]))-exp((-((l6[l]+l7[l])))*time[t])*((l4[l]+l5[l]-l8[l]))+exp((-((l4[l]+l5[l])))*time[t])*((l6[l]+l7[l]-l8[l]))))))/((((l4[l]+l5[l]-l6[l]-l7[l]))*((l4[l]+l5[l]-l8[l]))*((l6[l]+l7[l]-l8[l]))))) p66[l,t]<-(exp((((-l6[l])-l7[l]))*time[t])) p67[l,t]<-((l6[l]-exp((((-l6[l])-l7[l]))*time[t])*l6[l])/(l6[l]+l7[l])) p68[l,t]<-((exp((((-l6[l])-l7[l]))*time[t])*(((-1)+exp(((l6[l]+l7[l]-l8[l]))*time[t])))*l7[l])/(l6[l]+l7[l]-l8[l])) p88[l,t]<-(exp((-l8[l])*time[t])) p19[l,t]<-1-p11[l,t]-p12[l,t]-p13[l,t]-p14[l,t]-p15[l,t]-p16[l,t]-p17[l,t]-p18[l,t] p29[l,t]<-1-p22[l,t]-p23[l,t]-p24[l,t]-p25[l,t]-p26[l,t]-p27[l,t]-p28[l,t] p39[l,t]<-1-p33[l,t]-p34[l,t]-p35[l,t]-p36[l,t]-p37[l,t]-p38[l,t] p49[l,t]<-1-p44[l,t]-p45[l,t]-p46[l,t]-p47[l,t]-p48[l,t] p69[l,t]<-1-p66[l,t]-p67[l,t]-p68[l,t] p89[l,t]<-1-p88[l,t] } #end of t loop } #end of l loop # }TextControllers.StdCtrlDescTextControllers.ControllerDescContainers.ControllerDescControllers.ControllerDesc TextRulers.StdRulerDescTextRulers.RulerDescTextRulers.StdStyleDescTextRulers.StyleDescZTextRulers.AttributesDesc$ ZGo * ,[ @Documents.ControllerDesc Ws8 [h