1
我正在使用pmvnorm来计算R中的多元正态cdfs 但我在计算中获得负数 概率不能为负数。有人能帮助我找到了解决办法pmvnorm正在返回一个负数
的代码是在这里
library(mvtnorm)
pmvnorm(lower = rep(-Inf, 2), upper = c(-5, -18),
mean = c(0,0), sigma = matrix(c(1, -0.5, -0.5, 1),
nrow = 2, ncol = 2))
输出是
[1] -2.507297e-84
attr(,"error")
[1] 1e-15
attr(,"msg")
[1] "Normal Completion"
这是*勉强*负 – Dason
这是一个数字零。 – Roland
解决方法:'result < - pmax(0,result)'... –