:多変量正規desityに関連した二重積分を解く私は、既知の平均ベクトルと共分散行列を持つ多変量正規密度に関連した二重積分解決しようとしています
library(cubature)
mu1 <- matrix(c(3,3), nrow=2)
sigma1 <- rbind(c(4,-1), c(-1,6))
quadratic <- function(a,b) {
X <- matrix(c(a,b),nrow=2)
Q <- (-1/2)*t(X-mu1)%*%solve(sigma1)%*%(X-mu1)
}
NormalPDF <- function(x1,x2) {
f <- (1/(2*pi))*(1/sqrt(det(sigma1)))*exp(quadratic(x1,x2))
}
# Solving for P(1 < X1 < 3, 1 < X2 < 3)
P <- adaptIntegrate(NormalPDF(x1,x2), c(1,3), c(1,3))
しかし、それは続けて私にエラーを与えました:
Error in matrix(c(a, b), nrow = 2) : object 'x1' not found
私のコードに明らかなエラーはありますか?
は 'adaptIntegrate(NormalPDF、c(1,3)、c(1,3)) 'となります – HubertL