非中心卡方分布的逆矩

Question

我想计算R中非中心卡方分布的逆矩和截断逆矩。我怎样才能在R中做到？

假设X遵循具有自由度“k”和非中心性参数“t”的非中心卡方分布。我的问题是数值计算以下各种“t”值的期望值，这样我就可以模拟James-Stein型估计量的风险。

(i)   E[X^(-1)] and E[X^(-2)]
(ii)  E[X^(-1)I(A)] where I(A) is an indicator function of set A
(iii) E[1-c{X^(-2)}I(A)] where c is a constant.

Answer 1

Paolella＆＃39; book，中级概率给出了非中心卡方的各种力矩。见等式（10.10）。您可以在sadists package中找到这些代码的R代码。