Displaying 1 result from an estimated 1 matches for "sgamma".
Did you mean:
gamma
2011 May 31
0
rtmvt
I want to use the rtmvt from the {tmvtnorm} package using the "gibbs"
algorithm but how to i specify the nested function rtmvnorm to use gibbs as
well?
Right now I am using the code:
for (i in 1:g){
for (j in 1:n){
sgamma[,,i,j] = rtmvt(n=50, mean=mu[i,j], sigma[i,j],
df=nu[i], lower=rep(0,2),algorithm="gibbs")
}
}
heres an example of one iteration:
> mu[1,1]
-0.09734357 0.51578628
> sigma[1,1]
[,1] [,2]
[1,] 0.4250681 0.0253649
[2,] 0.0253649 0.4250681
when I run this i g...