CDOoDocuments.StdDocumentDescDocuments.DocumentDescContainers.ViewDescViews.ViewDescStores.StoreDescA,Documents.ModelDescContainers.ModelDescModels.ModelDescStores.ElemDesc+ +TextViews.StdViewDescTextViews.ViewDesc1+TextModels.StdModelDescTextModels.ModelDesc;)3)TextModels.AttributesDesc'* x Courier`h^ *`h $clj^X^# -------------------------------------------------------------------------------- # --------------- MODEL3: ZEROS TRICK ---------------------------------- #---full model --- model BivariatePoisson; { for (i in 1:n) {miny[i]<-min( y1[i], y2[i] )+1 } # # Likelihood C<- 0 for (i in 1:n){ u[i]~dunif(0, miny[i]); z3[i]<-trunc( u[i] ) z1[i]<-y1[i]-z3[i]; z2[i]<-y2[i]-z3[i]; zeros[i] <- 0 zeros[i] ~ dpois( zeros.mean[i] ) zeros.mean[i] <- -l[i] + C l[i]<- -lambda[i,1]+z1[i]*log(lambda[i,1])-loggam(z1[i]+1) -lambda[i,2]+z2[i]*log(lambda[i,2])-loggam(z2[i]+1) -lambda[i,3]+z3[i]*log(lambda[i,3])-loggam(z3[i]+1); for (k in 1:3){ log( lambda[i,k] ) <- beta[k,1] + beta[k,2]*x1[i] + beta[k,3]*x2[i] + beta[k,4]*x3[i] + beta[k,5]*x4[i] + beta[k,6]*x5[i] } } # # Priors for (k in 1:3){ for (j in 1:6) { beta[k,j]~dnorm(0.0, 0.01) } } } INITS list( beta = structure(.Data=c(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0), .Dim = c(3, 6)) ) # -------------------------------------------------------------------------------- # --------------- MODEL3: ZEROS TRICK ---------------------------------- #---full model --- model BivariatePoisson; { for (i in 1:n) {miny[i]<-min( y1[i], y2[i] )+1 } # # Likelihood C<- 0 for (i in 1:n){ u[i]~dunif(0, miny[i]); z3[i]<-trunc( u[i] ) z1[i]<-y1[i]-z3[i]; z2[i]<-y2[i]-z3[i]; zeros[i] <- 0 zeros[i] ~ dpois( zeros.mean[i] ) zeros.mean[i] <- -l[i] + C l[i]<- -lambda[i,1]+z1[i]*log(lambda[i,1])-loggam(z1[i]+1) -lambda[i,2]+z2[i]*log(lambda[i,2])-loggam(z2[i]+1) -lambda[i,3]+z3[i]*log(lambda[i,3])-loggam(z3[i]+1); log( lambda[i,1] ) <- beta[1] + beta[2]*x1[i] + beta[3]*x2[i] log( lambda[i,2] ) <- beta[4] + beta[5]*x2[i] + beta[6]*x5[i] log( lambda[i,3] ) <- beta[7]*x4[i] x3[i]~dnorm(0,1) } # # Priors for (j in 1:7){ beta[j]~dnorm(0.0, 0.01) } } INITS list( beta=c(0, 0, 0, 0, 0, 0, 0) ) DATA list(n=200, y1 = c(2, 0, 4, 2, 16, 1, 0, 3, 20, 1, 4, 2, 2, 5, 1, 4, 7, 6, 2, 6, 8, 2, 0, 1, 1, 2, 1, 3, 2, 2, 0, 3, 3, 2, 1, 5, 3, 4, 1, 1, 2, 1, 1, 6, 10, 4, 1, 0, 4, 0, 1, 0, 1, 2, 2, 3, 3, 3, 1, 5, 3, 7, 0, 5, 2, 2, 4, 2, 4, 2, 3, 1, 3, 1, 1, 0, 3, 1, 4, 1, 6, 1, 2, 3, 1, 4, 3, 2, 3, 2, 5, 3, 4, 3, 1, 5, 5, 3, 3, 3, 2, 1, 2, 3, 1, 1, 2, 5, 4, 3, 2, 4, 2, 2, 3, 17, 4, 3, 3, 2, 4, 3, 0, 8, 1, 1, 0, 4, 5, 6, 3, 2, 5, 4, 2, 3, 4, 9, 7, 1, 3, 5, 0, 3, 0, 1, 1, 4, 5, 2, 9, 4, 3, 3, 2, 4, 7, 3, 6, 4, 1, 5, 3, 3, 3, 4, 12, 8, 2, 3, 2, 5, 9, 1, 3, 2, 1, 2, 5, 5, 3, 2, 6, 1, 1, 2, 2, 0, 1, 17, 3, 6, 2, 3, 4, 3, 1, 2, 1, 5), y2 = c(5, 0, 3, 2, 16, 10, 0, 10, 19, 1, 5, 4, 5, 1, 2, 3, 7, 5, 9, 4, 7, 0, 0, 0, 0, 4, 2, 7, 0, 1, 1, 2, 1, 4, 2, 5, 2, 4, 2, 4, 4, 5, 1, 11, 10, 7, 0, 3, 4, 2, 2, 0, 4, 1, 0, 3, 4, 4, 1, 5, 2, 6, 3, 8, 2, 7, 3, 1, 2, 1, 3, 1, 1, 2, 3, 2, 2, 1, 5, 5, 4, 1, 2, 4, 1, 4, 2, 5, 2, 2, 7, 3, 7, 1, 2, 4, 3, 4, 4, 2, 5, 4, 1, 2, 1, 1, 1, 10, 3, 2, 2, 0, 4, 1, 5, 16, 1, 4, 1, 5, 5, 5, 5, 8, 5, 3, 3, 5, 5, 5, 5, 2, 3, 4, 6, 5, 6, 6, 7, 0, 6, 8, 1, 4, 1, 1, 3, 3, 4, 14, 4, 4, 4, 2, 11, 2, 5, 1, 4, 4, 1, 5, 5, 2, 2, 7, 19, 4, 4, 3, 6, 3, 10, 6, 7, 4, 3, 3, 4, 3, 2, 5, 10, 2, 0, 1, 2, 12, 2, 27, 5, 4, 3, 1, 3, 1, 2, 4, 8, 6), x1 = c(1.23, -0.16, -0.33, 0.97, 0.08, -0.05, 0.27, 0.32, 0.05, 0.22, 0.09, 0.17, -0.44, -1.26, 0.11, -0.36, 0.65, 0.5, -0.42, 0.14, -0.82, -0.43, 0.09, 0.73, 0.17, 1.05, -1.49, -0.83, 1.92, -0.33, -1.44, 0.06, 0.09, 0.19, -0.05, 0.17, -0.25, -0.58, 0.59, -0.12, 1.1, 1, -1.93, -0.12, -2.06, -0.19, -1.45, 1.72, -0.74, 0.52, -0.01, -0.03, -0.03, 1.26, 0.25, 1.39, 0.62, -0.56, 0.13, -0.85, 0.16, -1.29, -1.08, 1.24, -0.51, 0.28, 0.17, -1.01, 0.15, 1.38, 0.1, -0.05, -0.64, -1.5, -1.64, -0.51, 1.55, 1, -0.43, -0.88, 1.54, 0.55, -0.1, 0.49, -1.2, -0.39, 0.91, -1.21, -0.16, -0.76, -1.31,1.07, 0.21, 0.69, -0.23, -0.72, 1.63, -1.98, 1.21, -0.94, -2.89, 0.28, -0.07, 1.88, 1.52, 0.09, -0.47, 0.54, -0.01, 0.5, -0.53, 0.01, -0.3, 1.45, -0.25, -1.18, 1.59, -0.25, 0.13, 1.31, 1.84, -0.97, -0.56, 1.49, 0.26, -0.7, 1.38, -0.4, -0.87, 0.9, 0.18, 0.63, -1.05, 1.6, -1.74, 1.46, 0.54, 0.82, -1.2, -1.74, 0.5, 0.36, -1.15, 0.05, 2.14, -1.33, 0.65, 0.01, -0.24, 1.83, 2.39, 0.43, 0.22, -0.97, 0.01, 1.49, 0.16, -0.77, -0.89, 0.09, -0.11, 0.12, -0.03, -0.91, 0.81, -0.39, -2.67, 0.37, 0.06, -2.27, 1.5,0.47, 0.11, 0.01, 1.28, -1.17, -0.79, 0.01, 0.59, -0.81, 1.23, 0.66, -0.06, 0.06, -0.15, 0.52, -0.99, 0.06, 1.06, -1.37, 0.19, -1.27, 1.37, -1.09, 2.17, -0.24, -0.71, 0.55, -1.78, -0.24), x2 = c(-0.32, -0.03, -0.1, -0.39, 0.2, -0.98, 0.96, 0.08, 0.06, 0.67, 0.73, -1.64, 0.07, -1.43, -0.65, 1.68, -1.03, -0.33, -0.94, -0.63, -0.06, 0.16, 1.34, -0.97, -0.37, 0.26, -0.74, -0.54, 1.51, -0.06, -0.12, 0.38, 0.02, -0.29, -2.53, 0.29, 0.27, -0.53, -0.8, -1.14, -2.27, -1.88, 0.73, -1.21, -0.05, -0.54, 1.28, -0.52, 0.93, -1.48, 0.08, -0.02, -0.34, 0.53, 0.54, 0.75, -0.93, -0.53, 0.15, -0.52, 0.13, -0.5, 0.56, -1.12, -2.33, -1.11, -0.06, 1.62, 0.29, 0.82, -1.49, 1.17, -1.41, 0.52, 1.65, 0.84, -1.1, -0.46, 0.41, -0.04, 1.45, -0.46, -1.66, -1.45, 0.25, 0.96, -0.64, -0.47, 1.15, -0.54, -1.98, -0.16, -0.23, -0.42, 0.92, -0.09, 0.51, -0.22, -0.47, 1.77, 1.16, 1.03, -0.02, -1, 0.62, 0.94, 1.47, -1.53, -0.91, -1.3, -0.51, 0.58, -1.46, 1.59, -0.44, 1.2, -0.56, -0.43, -0.29, -0.7, -1.5, 0.74, 0.54, -0.7, 0.36, 0.78, 1.28, -0.32, 0.57, 0.9, -0.55, -0.57, 1.1, -0.71, -1.13, -1.88, 0.28, 0.45, 0.15, 0.24, -0.31, 0.78, 0.34, 0.73, -0.95, 1.06, 1.51, -2.15, 2.46, -2.73, -1.28, 1.31, -0.73, 0.06, -0.28, -0.35, -0.16, 0.81, 1.04, 0.99, 2.3, 0.32, -0.18, 1.36, -0.62, 0.78, -1.54, -0.68, -0.48, -1.15, -2.15, -0.86, 1.38, -1.28, -1.39, -0.46, -0.68, -0.44, 2.07, -1.37, -1.01, -1.53, 0.56, 0.49, -0.68, 0.71, -2.25, -0.54, 0.13, -1.55, 0.56, 1.05, -0.53, 1.62, -1.45, 2.23, 0.69, -0.44, -0.58, 0.53), x3 = c(1.02, -0.86, -0.95, 0.06, 2.14, 0.66, -0.62, 1.12, 1.48, 0.53, 0.85, -0.11, -0.38, 0.44, 1.79, 1.79, -0.07, -0.32, 0.46, 0.09, -0.16, 0.7, 1.35, -0.87, -0.85, 1.94, -0.56, 0.03, -0.87, -0.75, -0.07, 1.58, -1.63, 1.16, 0.84, -0.82, 1.35, 1.26, -0.71, 2.37, -0.27, 1.05, -1.55, 0.5, -0.91, -0.22, 1.1, 0.09, -0.56, -1.13, -0.88, -0.3, 1.04, -0.5, 1.16, 0.53, 0.36, 0.45, -0.22, 0.44, 0.46, 1.01, -0.35, 1.26, 0.89, -0.18, 2.05, -1.55, 1.1, -0.13, 0.04, -1.65, -0.84, -0.76, 1.28, 0.12, -0.28, 0.05, -0.49, -2.08, 0.3, -0.92, -1.18, 1.05, -1.37, -0.89, 0.99, -1.18, -0.26, -0.52, 0.26, -0.59, 0.86, 1.82, -0.43, 1.17, 0.43, 1.39, -0.63, 0.11, -0.48, 1.7, -2.43, 0.47, 0.36, -0.87, -1.3, -0.56, -1.18, -0.34, -0.04, 1.12, -0.26, 0.39, -0.78, -0.42, 0.88, 1.21, 1.37, -1.36, 0.12, -0.59, 1.35, -1.98, -0.65, 0.41, 1.02, 0.85, -0.83, -0.49, 0.47, 0.63, 0.59, -0.27, 0.13, 1.29, -0.08, 1.03, -0.86, 0.72, 2.25, -0.4, 0.19, -1.15, -0.91, 0.48, 0.89, 0.27, 0.47, 1.04, -1.53, -1.13, 1.8, -0.87, -1.68, -0.43, 0.29, 0.53, -0.15, 0.76, -1.51, -0.43, 0.63, 0.09, 0.83, -0.01, -1.58, 1.3, -1.07, 0.11, -0.63, 0.39, -1.14, 0.2, 1.16, 0.42, -0.41, -0.08, 1.04, -0.16, 1.09, 1.16, -0.05, 1.02, -1.48, -0.53, 1.42, 0.8, 0.49, -1.12, 0.02, 0.97, 1.33, 0.13, 0.07, 0.16, 0.95, -0.1, -0.86, -1.01), x4 = c(-0.45, -0.26, 1.26, -2.08, 2.44, -0.72, 0.09, -0.19, 3, -2.4, 0.75, 0.43, -1.12, -0.76, -0.06, 1.26, 1.71, -0.39, 1.51, 0.15, 0.78, -0.04, 0.15, -0.76, -2.13, -0.65, -0.44, -3.27, -0.15, 0.65, -2.29, -0.85, -0.34, 1, 0.91, 0.37, 1.12, -0.6, -1.36, -1.26, -0.7, -0.08, 0.56, 1.62, 1.67, 0.86, 0.46, -0.18, 0.2, 0.66, -0.78, -0.77, 0.07, -1.13, -0.89, 0.03, 0.56, 1.07, 0.97, 1.5, -0.23, 0.64, -1.27, 0.62, -1.13, -0.48, 0.64, -1.29, 0.97, -1.14, 0.24, -0.44, 0.16, 0.29, 0.05, -0.3, -0.22, 0.74, 0.67, 0.69, 1.33, -0.13, -0.26, 0.12, 0.6, 0.66, -0.18, -0.29, 0.31, -0.19, 1.15, 0.5, 0.87, -0.61, -1.56, -1.34, 0.92, 0.77, -0.19, -0.37, 0.64, -0.65, -1.31, -0.9, -0.08, -0.59, -0.54, 0.62, 0.11, -0.53, 0.68, -2.56, -0.18, 0.01, 1.23, 2.86, -0.51, -0.23, -0.76, -0.09, -0.36, 1.16, -0.25, 1.35, 0.3, 0.14, -0.06, 0.69, 0.8, 0.55, -0.28, -0.77, 1.53, -0.53, -1.2, -0.85, 1.01, 1.17, 0.57, 0.35, -0.03, 1.52, 0.23, -0.95, -2.03, -0.41, -0.8, -1.37, 0.8, 0.31, 0.27, 0.65, -0.52, 0.44, 0.36, 0.48, 1.2, -0.04, 0.85, 1.62, -0.16, -0.15, 0.44, -0.42, 0.48, 0.38, 1.79, 1.08, -0.29, 1.24, 1.09, -0.68, 2.25, -0.34, -0.68, 0.3, -0.34, -0.5, 1.56, -1.44, -2.06, -0.47, 0.49, -1.72, -0.75, -1.02, 0.73, -1.14, -0.68, 3.13, 0.39, 1.57, 1.3, -0.73, 0.93, -1.11, 0.33, -0.62, -1.25, 0.56), x5 = c(-0.58, 0.96, -0.04, -0.78, 0.8, -1.39, 0.28, -1.53, 1.27, 0.25, -0.89, -0.12, -0.27, 0.19, -0.51, -0.29, 0.65, -1.2, -1.36, 0.8, 0.48, 0.15, 0.35, 1.69, 2.11, -0.76, -0.22, 0.31, 2.12, 0.31, 0.43, 0.8, 0.28, -1.64, 0.09, 0.57, 1.75, -0.22, 0.8, -0.41, 2.65, -0.93, -0.13, -0.64, -0.54, -0.49, 1.14, -1.1, -0.91, 1.58, 0.52, 0.06, -1.32, 0.13, 1.44, -0.37, -0.89, 0.52, -0.41, 0.59, -0.88, -0.41, -0.63, 1.01, 0.04, -2.15, 0.04, 1.27, 1.35, 0.38, 2.87, 0.36, 1.58, 0.81, -0.66, -1.41, -0.05, 1.94, -0.86, -3.16, 1.26, -0.77, 1.21, -0.07, 0.4, -1.36, 0.78, -0.68, 0.23, -0.63, -0.95, -0.96, -1.26, 1.91, 0.45, -1.08, -0.84, -0.56, -0.96, 0.41, 0.01, -0.69, -1.53, 2.5, -0.99, 0.07, -0.41, -0.19, 1.27, -0.9, 0.46, 0.08, -0.86, -0.2, 1.38, 0.04, 0.74, -0.74, 1.33, -1.79, -0.45, 0.85, 0.02, 0.35, -1.34, -0.47, -0.72, 0.18, -0.75, -0.29, -1.76, -0.04, 0.15, -0.26, -0.74, -0.42, -0.15, 0.5, 0.99, 0.88, -0.74, 0.22, 1.5, -1.56, 1.23, -1.53, -1.17, 0.47, 2.12, -1.26, -0.24, -1.76, -1.59, -0.9, -3.28, -0.86, -0.15, -0.01, 0.72, 1.09, 0.64, -1.37, -0.15, -1.7, 0.09, -1.4, -0.7, 0.77, -1.54, 1.26, -0.14, 0.24, -1.24, -1.02, -0.77, -0.47, 1.47, -0.31, 0.94, -0.94, 0.32, -0.56, -0.87, 0.64, 1.21, 0.52, 1.03, -2.44, 0.07, -1.13, -0.56, 0, 0.53, -0.1, 0.11, 1.01, 0.84, 0.74, -0.97, -1.21) ) TextControllers.StdCtrlDescTextControllers.ControllerDescContainers.ControllerDescControllers.ControllerDesc TextRulers.StdRulerDescTextRulers.RulerDescTextRulers.StdStyleDescTextRulers.StyleDescZTextRulers.AttributesDesc$ ZGo * ,[ @Documents.ControllerDesc Ks ' Oh*