Hi, dear all, I wish to get one sample (2500-d vector) from the truncated multivariate normal distribution, so I choose use the R function rtmvnorm() to do this. But the error information shows that for this function, the dimension should be lower than 1000, So could you help me to find out if there's any solution could do such sampling from a truncated multivariate normal distribution with the sigma dimension is 2500*2500??
Here is parts of my code for sampling: rg[k,]<- rtmvnorm(n=1, as.vector(SG1u), SG1, lower=rep(0,2500), upper=rep(Inf,2500), algorithm="rejection") and the error information is: "Error: mvt(lower = lower, upper = upper, df = 0, corr = corr, delta = mean, : only dimensions 1 <= n <= 1000 allowed" Thanks a lot. [[alternative HTML version deleted]] ______________________________________________ R-help@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.