假设我有来自分布为Y = X ^ 2的群体的100个观察值的随机样本,其中X是正常的。如何使用R计算P(Y <= 2)的样本估计值?这里P()表示概率。以下代码是否有帮助?
X=rnorm(100)
Y=X^2
prob(Y<=2)
事先谢谢
答案 0 :(得分:0)
如果y = x^2
,则P(y <= 2) = P(x^2 <= 2) = P(x <= sqrt(2))
所以答案是:
pnorm(sqrt(2)) - pnorm(-sqrt(2))
# [1] 0.8427008
答案 1 :(得分:0)
感谢@markdly所需的代码是 总和(Y&LT; = 2)/长度(Y)
计算样本观察值小于或等于2的概率,该值随样本而变化。 @akond给出的上述答案基本上可以计算出确切的概率 pchisq(2,1)