Hello,
I think there's a bug in this line:
Instead of
precX <- rWishart(1,paramX+1,diag(nparamX))[,,1]
it should be, I guess:
precX <- rWishart(1,nparamX+1,diag(nparamX))[,,1]
(paramX is a matrix, not a number)
The same for this line also.
Again, I forgot to answer this comment following my correction. I also just realized that I had forgotten to make the correction in the iniParamXTr function. Now it is corrected
Hello, I think there's a bug in this line: Instead of
precX <- rWishart(1,paramX+1,diag(nparamX))[,,1]
it should be, I guess:precX <- rWishart(1,nparamX+1,diag(nparamX))[,,1]
(paramX is a matrix, not a number) The same for this line also.
I hope I'm seeing it right. Best, Miguel