R - 正则化不完全 beta 函数 w.r.t 的逆函数的导数。形状参数
R - Derivatives of the inverse of the regularized incomplete beta function w.r.t. shape parameters
正则化不完全 beta 函数 I(x,a,b) 是随机变量分布 Beta(a,b) 的 CDF。我有一个问题需要计算其反函数(分位数函数)w.r.t 的偏导数。参数 a 和 b。你知道怎么做吗?谢谢!
不知道有没有闭式解析表达式,不过你可以用数值来做(numDeriv::grad()
默认使用Richardson外推法,见文档):
x <- 0.5
## dummy (lambda) function with x from above and specified shape parameters
f <- function(y) qbeta(x, y[1], y[2])
library(numDeriv)
grad(f,
## shape parameters at which to evaluate partial derivatives
x = c(1,1)
)
## [1] 0.3465736 -0.3465736
?
正则化不完全 beta 函数 I(x,a,b) 是随机变量分布 Beta(a,b) 的 CDF。我有一个问题需要计算其反函数(分位数函数)w.r.t 的偏导数。参数 a 和 b。你知道怎么做吗?谢谢!
不知道有没有闭式解析表达式,不过你可以用数值来做(numDeriv::grad()
默认使用Richardson外推法,见文档):
x <- 0.5
## dummy (lambda) function with x from above and specified shape parameters
f <- function(y) qbeta(x, y[1], y[2])
library(numDeriv)
grad(f,
## shape parameters at which to evaluate partial derivatives
x = c(1,1)
)
## [1] 0.3465736 -0.3465736
?