Shantanu MULLICK
2012-Feb-05 14:29 UTC
[R] Simulating from a Normal Inverted Wishart distribution
Hello everyone I was wondering how would one simulate from a Normal Wishart Distribution in R. A normal inverted Wishart distribution is denoted by NIW (M,C,d,S), where X/(Sigma) ~ N( M,C,(Sigma) ) -> a matrix normal distribution, (Sigma) -> column dispersion matrix (Sigma) ~ IW (d,S) -> inverted Wishart distribution Thanks a lot ! Best Shantanu [[alternative HTML version deleted]]
Ranjan Maitra
2012-Feb-05 17:49 UTC
[R] Simulating from a Normal Inverted Wishart distribution
You can first simulate Sigma from the inverse Wishart distribution and then simulate from the matrix normal given the realized Sigma. As per a DuckDuckGo search, InvWishart function in the contributed R package MCMCpack may be what you need for the first step. HTH, Ranjan On Sun, 5 Feb 2012 15:29:33 +0100 Shantanu MULLICK <b00295766 at essec.edu> wrote:> Hello everyone > > I was wondering how would one simulate from a Normal Wishart Distribution > in R. > > A normal inverted Wishart distribution is denoted by NIW (M,C,d,S), where > > X/(Sigma) ~ N( M,C,(Sigma) ) -> a matrix normal distribution, (Sigma) -> > column dispersion matrix > > (Sigma) ~ IW (d,S) -> inverted Wishart distribution > > Thanks a lot ! > > Best > Shantanu > > [[alternative HTML version deleted]] > > ______________________________________________ > 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. >
Possibly Parallel Threads
- Parameterization of Inverse Wishart distribution available in MCMCpack and bayesm libraries
- Simulating from "matrix variate normal distribution"
- Wishart Density
- sample a matrix with one element to be 1 from wishart distribution
- R utilizing 25% of CPU for Dual core i3 370M processor