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]]