post47 <- MCMCordfactanal(~a1+a2+a3+a4+a5+a6+a7+a8+a9+a10+a11+a12+a13+a14+a15+a16+a17+a18+a19+a20+a21+a22+a23+a24+a25+a26+a27+a28+a29+a30+a31+a32+a33+a34+a35+a36+a37+a38, data=un47,lambda.constraints=list(a8=list(2,"+")), factors=1, burnin=5000, mcmc=200000, thin=40, verbose=1000, L0=0.05, store.lambda=TRUE, store.scores=TRUE, tune=1.2) I have a question regarding setting lambda.constraints in the code above. What is the criteria one should use for setting lambda constraint? For example, why would you chose to set the 2nd factor loading of a8 as + instead of choosing a different variables and factor loadings and setting it to negative? Above is a code for ordinal MCMC factanal model. Thank you very much in advance and I hope to hear from you soon! -- View this message in context: http://r.789695.n4.nabble.com/Setting-Lambda-constraints-for-MCMC-ordinal-model-tp3313777p3313777.html Sent from the R help mailing list archive at Nabble.com.