CDOoDocuments.StdDocumentDescDocuments.DocumentDescContainers.ViewDescViews.ViewDescStores.StoreDescDocuments.ModelDescContainers.ModelDescModels.ModelDescStores.ElemDescg _TextViews.StdViewDescTextViews.ViewDesc TextModels.StdModelDescTextModels.ModelDesc TextModels.AttributesDesc'*  */w # --------------- ZI-NegBin: GSS1990 DATA ------------------------------------ # --------------- GLM BASED APPROACH ------------------------------------ # --------------- LATENT VAR IMPLEMENTATION ------------------------------------ # 1/2/2012 Variance expression were corrected; mistake was reported by Gheorghe Doros # model{ for(i in 1:n){ y[i] ~ dnegbin( p.zinb[i], r.ind[ i ] ) p.zinb[i] <- r.ind[i]/( r.ind[i]+lambda.ind[i]*(1-u[i]) ) log(lambda.ind[i]) <- beta[1] + beta[2] * gender[i] r.ind[i] <- r[ gender[i] + 1 ] u[i]~dbern( p0[ gender[i] + 1 ] ) l[i] <- log( p0[ gender[i]+1 ] * equals( y[i], 0 ) + (1-p0[ gender[i]+1 ])*fd[i] ) # set the probability function of each distribution # this is the actual p.ind p.ind[i] <- r.ind[i]/( r.ind[i]+lambda.ind[i] ) lfd[i] <- loggam( y[i]+r.ind[i] ) - loggam( r.ind[i] ) - loggam( y[i]+1 ) + r.ind[i]*log( p.ind[i] ) + y[i]*log( 1-p.ind[i] ) fd[i] <- exp( lfd[i] ) cond.l[i] <- lfd[i]*(1-u[i]) } lambda[1] <- exp(beta[1]) lambda[2] <- exp(beta[1]+beta[2]) expbeta2 <- exp(beta[2]) for (j in 1:2){ p0[j] ~ dbeta(1,1) beta[j] ~ dnorm( 0.0, 0.001 ) r[j]~dgamma(0.001,0.001) p[j] <- r[j]/( r[j]+lambda[j] ) y.mean[j] <- (1-p0[j])*lambda[j] y.var[j] <- (1-p0[j])*( lambda[j]*( lambda[j]+r[j] )/r[j] + p0[j]*lambda[j]*lambda[j] ) di[j] <- y.var[j]/y.mean[j] } Deviance <- -2*sum(l[1:n]) cond.dev <- -2*sum(cond.l[1:n]) } # --------------- ZI-NegBin: GSS1990 DATA ------------------------------------ # --------------- GLM BASED APPROACH ------------------------------------ # --------------- ZEROS IMPLEMENTATION ------------------------------------ model{ C<-0 for(i in 1:n){ zeros[i] <- 0 zeros[i] ~ dpois( zeros.means[i] ) zeros.means[i] <- -l[i] + C # log-likelihood for i individual l[i] <- log( p0[ gender[i]+1 ] * equals( y[i], 0 ) + (1-p0[ gender[i]+1 ])*fd[i] ) # set the probability function of each distribution lfd[i] <- loggam( y[i]+r.ind[i] ) - loggam( r.ind[i] ) - loggam( y[i]+1 ) + r.ind[i]*log( p.ind[i] ) + y[i]*log( 1-p.ind[i] ) fd[i] <- exp( lfd[i] ) p.ind[i] <- r.ind[i]/( r.ind[i]+lambda.ind[i] ) r.ind[i] <- r[ gender[i] + 1 ] log(lambda.ind[i]) <- beta[1] + beta[2] * gender[i] } lambda[1] <- exp(beta[1]) lambda[2] <- exp(beta[1]+beta[2]) expbeta2 <- exp(beta[2]) for (j in 1:2){ p0[j] ~ dbeta(1,1) beta[j] ~ dnorm( 0.0, 0.001 ) r[j]~dgamma(0.001,0.001) p[j] <- r[j]/( r[j]+lambda[j] ) y.mean[j] <- (1-p0[j])*lambda[j] y.var[j] <- (1-p0[j])*( lambda[j]*( lambda[j]+r[j] )/r[j] + p0[j]*lambda[j]*lambda[j] ) di[j] <- y.var[j]/y.mean[j] } Deviance <- -2*sum(l[1:n]) } INITS list( beta=c(0,0), p0=c(0.1,0.1), r=c(1,1) ) DATA list( n =550, y = c(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 7, 7, 7, 7, 7, 7, 7, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 9, 9, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 12, 12, 12, 12, 12, 12, 13, 13, 13, 15, 15, 15, 16, 16, 16, 20, 20, 20, 20, 20, 20, 20, 24, 25, 30, 30, 30, 50, 60, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 7, 7, 7, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 9, 9, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 13, 13, 13, 14, 15, 15, 15, 15, 15, 15, 15, 15, 15, 15, 16, 17, 18, 20, 20, 20, 20, 20, 20, 22, 23, 25, 25, 25, 27, 30), gender = c(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1) ) TextControllers.StdCtrlDescTextControllers.ControllerDescContainers.ControllerDescControllers.ControllerDesc TextRulers.StdRulerDescTextRulers.RulerDescTextRulers.StdStyleDescTextRulers.StyleDescZTextRulers.AttributesDesc$ Zo * ,[ @Documents.ControllerDesc Ws,! [h$