確率変数の期待値をRで求めたいのですが、、、

質問

今学校で出されている課題で、

ｆ（ｘ）＝（3/10*x^2+1/5*x)*exp(-x) (x>0)

の期待値をRで求めよという課題が出されているのですが、どう書いたら出せるか教えていただけないでしょうか。これまでは標本平均をmeanでしか出したことがなく、また、ありふれた分布であれば、公式から出していたのですが、見たことない分布についてRでどのように出せるかという基本的なことについて教えていただければありがたいです。

回答No.1

noname#227064 · ベストアンサー

Rにはintegrateという関数がありますので、これを使えば
f <- function(x) ifelse(x > 0, (3/10*x^2+1/5*x)*exp(-x), 0)
integrate(function(x) x*f(x), -Inf, Inf)
とするだけですが、そのf(x)は確率密度関数ではありません。
式の写し間違い等はしていませんか？
（1/5ではなくて2/5だったとか）

お礼

quaestioさん

写し間違いしてて、quaestioさんの言うように2/5でした。
これでは、積分して１にならないですね。
迅速にまた丁寧に回答いただいて本当感謝してます。この分野は得意でないので、またわからない点あったら教えていただけると有難いです。お蔭様で助かりました。