CDOoDocuments.StdDocumentDescDocuments.DocumentDescContainers.ViewDescViews.ViewDescStores.StoreDesc Documents.ModelDescContainers.ModelDescModels.ModelDescStores.ElemDesc TextViews.StdViewDescTextViews.ViewDesc| TextModels.StdModelDescTextModels.ModelDesc~TextModels.AttributesDesc1$Courier New  $Courier NewL #------------------------------------------------------------------------------------------------------------------------------ Version 1 #------------------------------------------------------------------------------------------------------------------------------ model{ # model's likelihood for (a in 1:L){ for (b in 1:L){ y[a,b] ~ dlnorm( mu[a,b], tau ) mu[a,b] <- log(inf[a+b-1]) + mu0 + alpha[a] + beta[b] }} #### CR Constraints alpha[1] <- 0.0 beta[1] <- 0.0 # priors mu0~dnorm( 0.0, 1.0E-04) for (a in 2:L){ alpha[a]~dnorm( 0.0, 1.0E-04) beta[a] ~dnorm( 0.0, 1.0E-04) } tau ~dgamma( 0.01, 0.01) s <- sqrt(1/tau) # precision for ( i in (L+1):(2*L-1) ){ for ( a in 1:(i-L) ){Y.T[i-L,a] <- 0.0} for ( a in (i-L+1):L ){Y.T[i-L,a] <- y[ a, i+1-a ] } T[i-L]<-sum( Y.T[i-L,1:L] ) } Total <- sum(T[]) } INITS list( mu0=1.0, alpha=c(NA, 0,0,0,0,0,0), beta=c(NA, 0,0,0,0,0,0), tau=1.0 ) DATA (LIST) list( L=7, inf=c(1.000 , 1.204 , 1.439 , 1.666 , 1.906 , 2.142 , 2.356, 2.484 ,2.618, 2.736, 2.794, 2.875, 2.981), y = structure( .Data=c( 1546.6 , 647.5 , 382.2 , 246.8 , 211.8 , 63.6 , 146.3 , 2099.0 , 1001.8 , 487.2 , 293.0 , 318.8 , 269.9 , NA , 3422.2 , 1257.1 , 488.4 , 456.0 , 562.4 , NA , NA , 4948.8 , 1899.7 , 984.0 , 1253.3 , NA , NA , NA , 8161.3 , 2820.3 , 1304.8 , NA , NA , NA , NA , 10622.0 , 3897.7 , NA , NA , NA , NA , NA , 11744.9 , NA , NA , NA , NA , NA , NA ), .Dim = c( 7,7 ) ) ) TextControllers.StdCtrlDescTextControllers.ControllerDescContainers.ControllerDescControllers.ControllerDesc TextRulers.StdRulerDescTextRulers.RulerDescTextRulers.StdStyleDescTextRulers.StyleDescZTextRulers.AttributesDesc$ ZGo * ,[ @Documents.ControllerDesc t]s ' `h*