[R] R-Hierarchical binomial proportion modelling

[zaid bhatti]

Hi, 
i am running a model without covariates, along with shape file and
neighbourhood matrix,
i am facing difficulty to edit this code for 34 provinces, can any help me
to edit this below given code.



FinalModels_WINScriptToRunModel_AllMetrics v2 final_20141205

          inits1<-list(precu = 10,precv = 10,alpha =-1,
          v=rep(1,Ndistricts), Bprov=c(NA,/0.1,0.2,0.3,0.4,0.5,0.6,0.7/),
u=rep(1,Ndistricts))                   
          
          inits2<-list(precu = 20,  precv = 20, alpha =-0.9,
          v=rep(-1,
Ndistricts),Bprov=c(NA,/0.01,0.02,0.03,0.04,0.05,0.06,0.07/),
u=rep(-1,Ndistricts))
          inits<-list(inits1,inits2)



Model district dbin_3

Bprov[1] <- 0
  for(j in 2:34){
        Bprov[j] ~ dnorm(0, 0/.0000001/)   ### how should we edit this ?
  }





--
View this message in context: http://r.789695.n4.nabble.com/R-Hierarchical-binomial-proportion-modelling-tp4701239.html
Sent from the R help mailing list archive at Nabble.com.