作者kria5304 (XenoMegaREENovaSaga)
看板Statistics
标题[问题] 连续型机率密度函数抽样算均数
时间Wed Feb 6 16:27:02 2019
一个观念问题 感觉有点基本但太久没碰统计有点生疏(遮脸...
题目要从λ=2指数分配抽样10000笔 然後算其均数 变异数
答案应该是趋近於(λ,λ^2)吧
原本我是直接取 x = (1,100,step=0.01) 丢进机率密度函数转成机率f(x)
然後 sum(x*f(x))
结果就炸了 f(x)自己累加就超过1了...
自己想想 上面的作法是离散吧 连续感觉不能这样
可是就一直卡在这跳不出去XDD
可以请大神开导一下吗 谢谢!!
--
恩?你说去去武器走之类的很弱 想学强大一点的咒语?
霹雳卡霹雳拉拉波波莉纳贝贝鲁多
这个你觉得怎样?
--
※ 发信站: 批踢踢实业坊(ptt.cc), 来自: 118.168.175.173
※ 文章网址: https://webptt.com/cn.aspx?n=bbs/Statistics/M.1549441628.A.87E.html
1F:推 tnlomtt: 机率积分转换? 02/06 18:10
2F:→ kria5304: 是说*J吗 好像有可能诶 回去再算算看 谢谢 02/06 18:57
3F:推 tnlomtt: 用PIT从uniform生成指数乱数 02/06 21:36
4F:→ Pieteacher: 用 -log(U) 为 Exp(1) where U ~ U(0,1) 在调一下 sca 02/06 22:08
5F:→ Pieteacher: le parameter 02/06 22:08
6F:→ recorriendo: 绍面两楼说的是蒙地卡罗法 你做的是直接积分但是你忘 02/07 06:04
7F:→ recorriendo: 了成dx 可以想想如果你的step更精细 得到更多组f(x) 02/07 06:07
8F:→ recorriendo: 全部加起来不就更大? 所以显然是缺一个normalization 02/07 06:08
9F:推 goshfju: 用excel就可以啦 先生成U(0,1)乱数:u1,...,u10000 02/07 12:28
10F:推 goshfju: u=F(x)=1-e^(-2x) , x = -ln(1-u)/2,可算出x1,...,x10000 02/07 12:30
11F:→ goshfju: 计算xbar,Sx^2 然後我看不懂你在干吗 你的命题不是要算 02/07 12:31
12F:→ goshfju: 样本平均数跟样本变异数? 02/07 12:31
13F:→ goshfju: 如果你要用pdf去算,根本不用模拟阿 02/07 12:32
14F:→ goshfju: E(X)=∫xf(x)dx=1/λ, E(X^2)=∫x^2f(x)dx 02/07 12:33
15F:→ goshfju: Var(X)=E(X^2)-(E(X))^2=1/λ^2 02/07 12:33