如何计算R

Question

我编写了以下代码，从文件中读取数据（有四个参数），并计算多变量分布的均值和协方差矩阵。参数mu和sigma使用MLE（最大似然法）计算。

library(mvtnorm)
df <- read.csv("../data/dataset.train", header=FALSE)
mu = colMeans(df)
sigma <- matrix(0, nrow = ncol(df), ncol = ncol(df))
for(row in 1:nrow(df)) {
  temp = df[row,]-mu
   sigma = sigma + as.matrix(t(temp)) %*% as.matrix(temp)
}
sigma = sigma / nrow(df)

现在我想找到以下概率$ P（a1

我的努力： 我尝试过使用

pmvnorm(lower=-Inf, upper=Inf, mean=rep(0, length(lower)),corr=NULL, sigma=NULL, algorithm = GenzBretz(), ...) method

将下限和上限放在

之后

lower = [a1,b1,c1,d1] and upper = [a2,b1,c1,d2]

但它给出了错误＆＃34; lower == upper＆＃34;。我认为这是因为我的第二和第三个值在低位和高位是相同的。 但我不知道如何在R中计算出来。

请帮帮我。提前谢谢。