data modeling - WinBUGS error: vector valued relation z must involve consecutive elements of variable -


i trying model multivariate probit model binary data. have been trying winbugs in return gives me error. ideas or suggestion warmly welcomed.

model{ (i in 1:ns){ ## loop on studies

for (k in 1:2){   ### loop on arm   (j in 1:2){  ###  loop on outcomes      r[i,k,j] ~ dbin(p[i,k,j],n[i,k,j]);             ## likelihood function       p[i,k,j] <-  phi(z[i,k,j])     z[i,k,1:2] ~  dmnorm(theta[i,1:2],with[i,,])i(-5, 5)    #latent variable (z<0) or probit link            theta[i,1] <- alpha[i,k,1] + beta[i,k,1]    theta[i,2] <- alpha[i,k,2] + beta[i,k,2]   }   ###close  loop on outcomes     }   ###close  loop on arms  alpha[i,2,1] <- 0    alpha[i,2,2] <- 0  alpha[i,1,1:2] ~ dnorm(0,.0001) beta[i,2,1:2] ~  dmnorm(d[1:2],prec[,])     beta[i,1,1] <- 0    beta[i,1,2] <- 0    ## priors on within study cov matrix with[i,1:2,1:2] <- inverse(cov.mat[i,1:2,1:2])        #define elements of within-study covariance matrix   cov.mat[i,1,1] <-  1   cov.mat[i,2,2] <-  1   ### prior ipd data  ######   cov.mat[i,1,2] ~   dbeta(a[i],b[i])   cov.mat[i,2,1] <-  cov.mat[i,1,2]    a[i]<-31.97   b[i]<- 4.52                     }#### close loop on studies         (i in 1:2) {   d[i]  ~ dnorm(0.0000e+00, 0.0001)    # overall treatment effects } ## priors on between study cov matrix prec[1:2,1:2]<-inverse(tau[1:2,1:2]) pi<-3.14/2 a1~dunif(0, pi) rho.tau<-cos(a1) sd[1]~dunif(0,2) sd[2]~dunif(0,2) tau[1,1]<-pow(sd[1],2) tau[2,2]<-pow(sd[2],2) tau[2,1]<-tau[1,2] tau[1,2]<-sd[1]*sd[2]*rho.tau }        #end model

these data:

list(ns=2) t[,1,1] t[,1,2] t[,2,1] t[,2,2]  r[,1,1]  n[,1,1] r[,2,1] n[,2,1] r[,1,2]     n[,1,2] r[,2,2] n[,2,2] 1   0   1   0   19  77  23  77  60  82  70  82 1   0   1   0   27  199 54  199 231 393 318 393 end

the model syntactically correct , allows me load data. once compile error in title. thank given

it looks if inputting 2 2 matrix mean of multivariate normal distribution right here.

z[i,1:2,k] ~  dmnorm(theta[i,,],with[i,,])i(-5, 5)    #latent variable (z<0) or probit link

however, appears if z vector of length 2. need input vector mean of dmnorm , give associated variance covariance matrix (i.e. if supply vector of length 3, needs have 3 3 variance covariance matrix). right have 2 2 matrix input mean (4 parameters) , 2 2 variance covariance matrix. not know motivation behind model can't provide suggestions on how fix per se, seems me need index theta bit more in order prevent putting matrix dmnorm.


-

Comments

Post a Comment

Popular posts from this blog

many to many - Django Rest Framework ManyToMany filter multiple values -

i have 2 models, 1 defining users, other defining labels on these users. using django rest framework create api. able query users @ least labels ids 1 , 2. for instance users' labels are: [(1,2), (1,2,3), (2,3), (1,3)] want query return [(1,2), (1,2,3)] . so far, i've managed query users given label (let's id=1) doing: /api/users/?labels=1 , unable query users labels 1 , 2. i've tried /api/users/?labels=1,2 or /api/users/?labels=1&labels=2 return invalid users, i.e. users without labels 1 or 2... any welcome. thanks, dimitry github test repo: https://github.com/thedimlebowski/drf-m2m-filter code: models.py class label(models.model): name = models.charfield(max_length = 60) class user(models.model): labels = models.manytomanyfield(label) filters.py class userfilter(django_filters.filterset): labels = django_filters.filters.baseinfilter( name='labels', lookup_type='in', ) class meta:...
Read more

java - Jasper subreport showing only one entry from the JSON data source when embedded in the Title band -

Image
i having issues populating data in subreport in title band of main report. data populated correctly in subreport in detail band of main report. trying find wrong in json query. appreciated. json sample data source (expenses.json) { "expenses": { "date": "8 sep 2016", "accounts": [ { "title": "xyz corp (111)", "accountname": "xyz corp", "accountnumber": "111", "transactions": [ { "date": "21 jun 2016", "name": "gas", "price": "17.50" }, { "date": "12 may 2016", "name": "shopping", "price": "111.99" }, { "date": "30 apr 2016", "name...
Read more

Java Entity Manager - JSON reader was expecting a value but found 'db' -

i have code working: entitymanager entitymanager = getentitymanager(); entitymanager.gettransaction().begin(); string query = "db.band.find({})"; list<band> list = (list<band>) entitymanager.createnativequery(query, band.class).getresultlist(); entitymanager.close(); return list; it returns list no problem. want sort list date : entitymanager entitymanager = getentitymanager(); entitymanager.gettransaction().begin(); string query = "db.band.find({something__id:objectid(\"" + myid + "\")}).sort({\"somethingelse.date\":-1})"; list<band> list = (list<band>) entitymanager.createnativequery(query, band.class).getresultlist(); entitymanager.close(); return list; my console gives me message: org.bson.json.jsonparseexception: json reader expecting value found 'db' i checked string working in mongo console , did. ideas? edit: tried putting "something__id" between quotes because would...
Read more