Dear all members The following function is truncated normal distribution function i wrote it by using winbugs so now i want to write it by using R - program.
for(j in 1:P){y[i,j]~djl.dnorm.trunc(mu[i,j],psi[j],thd[j,z[i,j]],thd[j,z[i,j]+1])} where y:underlying continuous variables and treated as a missing variable and augmented wit the posterior analysis. mu: mean psi: variance z:ordered categorical variables the : preassigned values for thresholds and represent matrix(10*6) dimention. I hope anyone help me Many thanks in advance Thanoon [[alternative HTML version deleted]] ______________________________________________ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.