我想解一个方程如下,其中X
是唯一的未知变量,函数f()
是一个多元学生t
分配。更准确地说,我有一个学生密度函数的多 k 维积分,它给了我们一个概率,我知道这个概率是 q。所有积分的下限是 -Inf 并且我知道最后一个 k-1 维度的上限(如给定的),唯一未知的变量是第一个积分的上限。它应该有一个变量和一个方程的解。我试图在 R 中解决它。我做了动态条件相关以获得相关矩阵,以指定我的 t 分布。因此,将此相关矩阵插入我的多 t 分布“dmvt”,并使用“cubature”包中的“adaptIntegral”函数构造一个函数作为命令“uniroot”的参数,以求解第一个积分的上限。但是我有一些困难来实现我想要的。(我希望我的问题很清楚)我之前提供了我的代码,有人告诉我有问题,但找不到那里有问题的原因。非常感谢您的帮助。
我现在如何用一维积分来处理它,但我不知道如何在 R 中求解多维积分方程?(例如二维案例)
\int_{-\infty}^{X}
\int_{-\infty}^{Y_{1}} \cdots
\int_{-\infty}^{Y_{k}}
f(x,y_{1},\cdots y_{k})
d_{x}d_{y_{1},}\cdots d_{y_{k}} = q
此代码失败:
require(cubature)
require(mvtnorm)
corr <- matrix(c(1,0.8,0.8,1),2,2)
f <- function(x){ dmvt(x,sigma=corr,df=3) }
g <- function(y) adaptIntegrate(f,
lowerLimit = c( -Inf, -Inf),
upperLimit = c(y, -0.1023071))$integral-0.0001
uniroot( g, c(-2, 2))