我需要将一个循环与函数mvtnorm结合起来。 使用pmvnorm可以轻松完成：例如，在这里，我根据平均值的不同值创建一个具有不同累积概率值的向量：
a <- c(0, 0, 0)
y <- c(0, 0, 0)
for(i in 1:3)
a[i] <- i*0.5
for(i in 1:3)
y[i] <- pnorm(2, a[i], 1)
y
# [1] 0.9331928 0.8413447 0.6914625
但是，如果我使用mvtnorm怎么办？我想有一个双变量概率的向量，两个分布中的一个的平均值取不同的值。如何编写mu命令？例如
library(mvtnorm)
a <- c(0, 0, 0)
for(i in 1:3)
a[i] <- i*0.5
mu <- c(0, 0)
for(i in 1:2)
mu[i] <- (1, a)
......等等。我被mu命令阻止了。谁能帮我？谢谢