r_user 发表于 2012-7-23 17:50 
谢谢你的回复。此外还有2个疑点想你帮我解答下,谢谢了
sigma 指 covariance matrix of dimension n
me ...
额 这个很好理解啊。。你把对这个函数的input定义好
如 mean=c(1,2); 这个不是计算出来的 是作为input给他的
sigma 指的是 the covariance matrix of dimension n. Either corr or sigma can be specified. If sigma is given, the problem is standardized. If neither corr nor sigma is given, the identity matrix is used for sigma.
例如我上面回复给的link的一个例子就是:
a <- pmvnorm(lower=-Inf, upper=c(2,2), sigma = diag(2)*2)
这里sigma不是为一的对角阵了吧
所以这里就是随便什么对角阵就好 不一定是1;
如果用corr作为input的话就是