循环和mvtnorm?

时间:2016-01-16 17:14:40

标签: r cycle

我需要将一个循环与函数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命令阻止了。谁能帮我?谢谢

0 个答案:

没有答案