[R] rtmvt

statfan irene_vrbik at hotmail.com
Mon Apr 11 02:52:24 CEST 2011


I have been using the rtmvt function in the {tmvtnorm} package i'm getting
the warning:

"Acceptance rate is very low and rejection sampling becomes inefficient.
Consider using Gibbs sampling."

but i AM specifying the gibbs algorithm!!:

rtmvt(M, mean=q[,,i,j], sigma=((u[i,j] + nu[i])/(p+nu[i]))*delta[,,i],
			df=ceiling(nu[i]+p), lower=c(0,0), algorithm="gibbs")

Any ideas why I am getting this warning and how I can fix it?  

--
View this message in context: http://r.789695.n4.nabble.com/rtmvt-tp3440751p3440751.html
Sent from the R help mailing list archive at Nabble.com.



More information about the R-help mailing list