Displaying 1 result from an estimated 1 matches for "gibbsr".
Did you mean:
gibbs
2011 Feb 07
1
Question about checkTmvArgs function in rtmvnorm (package tmvtnorm)
...owever, it is NOT a correlation matrix).
Thanks a lot!
Dimitri
function (n, mean = rep(0, nrow(sigma)), sigma = diag(length(mean)),
lower = rep(-Inf, length = length(mean)), upper = rep(Inf,
length = length(mean)), algorithm = c("rejection", "gibbs",
"gibbsR"), ...)
{
algorithm <- match.arg(algorithm)
cargs <- checkTmvArgs(mean, sigma, lower, upper)
mean <- cargs$mean
sigma <- cargs$sigma
lower <- cargs$lower
upper <- cargs$upper
if (n < 1 || !is.numeric(n) || n != as.integer(n) || length(n) >...