[R] rtmvt

Tóth Dénes tdenes at cogpsyphy.hu
Mon Apr 11 11:16:24 CEST 2011

Hi there,

Since you failed to provide us with data and sessionInfo(), I can only
guess that for some reason you call the rtmvt.rejection function instead
of rtmvt.gibbs.
Just look at the code of rtvmt by typing:

There you can see that it is a wrapper for rtmvt.rejection or rtmvt.gibbs.
You can run them directly by typing:

... shall be the arguments as pre-processed by rtmvt().

But first, you might upgrade your base R and the installed packages.


> 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.
> ______________________________________________
> R-help at 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.

More information about the R-help mailing list