Dear Jags users,
For a dinterval() example shown below, since jags gives the wrong deviance estimate (always 1), is it correct to use pnorm function to calculate the deviance? The example:
y[i]~dinteral(t[i],limit) # y is an ordinal variable with categories 0,1,2; limit has two values -10,10.
To calculate the deviance, we need to calculate p(t[i]<-10),p(-10<t<span>[i]<10) and p(t[i]>10). Since t[i] follows a normal distribution, we could calculate these probabilities using pnorm function in jags. Correct?
Then what funciton should we use for a multivariate case? an example:
In this case, how are we going to calculate all the probabilities? Is there a multivariate pnorm function, like pmvnorm in R package 'mvtnorm'?
Any suggestions will be appreciated.
Log in to post a comment.
Sign up for the SourceForge newsletter:
You seem to have CSS turned off.
Please don't fill out this field.