24小时热门版块排行榜    

查看: 1145  |  回复: 4

幻魇

新虫 (初入文坛)

[交流] 求教一个利用累计概率函数拟合分布参数的问题 已有1人参与

http://cn.mathworks.com/help/stats/examples/fitting-a-univariate-distribution-using-cumulative-probabilities.html#zmw57dd0e2006
我在看上边matlab的文档里边,用CDF估计参数的部分。这里边实际上用了一个最小化平方误差的方法,但是里边有个加权wgt = 1 ./ sqrt(pEmp.*(1-pEmp));不是很懂,请大神们指教是不是出自某篇文章里的?具体为什么这样做会更好。

求教一个利用累计概率函数拟合分布参数的问题
QQ图片20151103144338.png
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

feixiaolin

荣誉版主 (文坛精英)

优秀版主


小木虫: 金币+0.5, 给个红包,谢谢回帖
To fit a Weibull distribution to these data, notice that the CDF for the Weibull is p = Pr{X <= x} = 1 - exp(-(x/a)^b). Transforming that to log(a) + log(-log(1-p))*(1/b) = log(x) again gives a linear relationship, this time between log(-log(1-p)) and log(x). We can use least squares to fit a straight line on the transformed scale using p and x from the ECDF, and the slope and intercept of that line lead to estimates of a and b.
人家已经讲清楚啦
2楼2015-11-03 20:26:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

幻魇

新虫 (初入文坛)

引用回帖:
2楼: Originally posted by feixiaolin at 2015-11-03 20:26:57
To fit a Weibull distribution to these data, notice that the CDF for the Weibull is p = Pr{X <= x} = 1 - exp(-(x/a)^b). Transforming that to log(a) + log(-log(1-p))*(1/b) = log(x) again gives a li ...

多谢回复,不过我指的是wgt = 1 ./ sqrt(pEmp.*(1-pEmp))这个权重,为什么这么取?
3楼2015-11-03 21:12:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

feixiaolin

荣誉版主 (文坛精英)

优秀版主

引用回帖:
3楼: Originally posted by 幻魇 at 2015-11-03 21:12:27
多谢回复,不过我指的是wgt = 1 ./ sqrt(pEmp.*(1-pEmp))这个权重,为什么这么取?...

老眼昏花了

[ 发自小木虫客户端 ]
4楼2015-11-03 23:08:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

幻魇

新虫 (初入文坛)

引用回帖:
4楼: Originally posted by feixiaolin at 2015-11-03 23:08:53
老眼昏花了
...

你觉得那个权重为什么这么取呢?
5楼2015-11-04 13:52:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 幻魇 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[论文投稿] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +5 gy1nBQXYQJqL 2026-08-29 6/300 2026-08-31 01:37 by hZiFeudZyoGR
[博后之家] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +7 gy1nBQXYQJqL 2026-08-29 9/450 2026-08-31 01:15 by hZiFeudZyoGR
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +4 gy1nBQXYQJqL 2026-08-29 8/400 2026-08-31 00:45 by hZiFeudZyoGR
[考博] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +7 ASdOkHsho7FD 2026-08-28 12/600 2026-08-31 00:34 by hZiFeudZyoGR
[考博] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +4 ASdOkHsho7FD 2026-08-28 7/350 2026-08-31 00:21 by hZiFeudZyoGR
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +6 ASdOkHsho7FD 2026-08-28 8/400 2026-08-31 00:21 by hZiFeudZyoGR
[找工作] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 4FFAWE8HcgUD 2026-08-29 5/250 2026-08-30 20:29 by hZiFeudZyoGR
[考博] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 4FFAWE8HcgUD 2026-08-29 3/150 2026-08-30 20:19 by hZiFeudZyoGR
[博后之家] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +4 jCd0dEvKHShX 2026-08-29 5/250 2026-08-30 19:46 by hZiFeudZyoGR
[考博] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +4 gy1nBQXYQJqL 2026-08-29 6/300 2026-08-30 18:20 by hZiFeudZyoGR
[硕博家园] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +5 gy1nBQXYQJqL 2026-08-29 6/300 2026-08-30 18:17 by hZiFeudZyoGR
[论文投稿] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +6 ASdOkHsho7FD 2026-08-28 7/350 2026-08-30 17:03 by hZiFeudZyoGR
[找工作] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +6 ASdOkHsho7FD 2026-08-28 9/450 2026-08-30 11:19 by l0VvVHGBGRLv
[基金申请] 2026年叶企孙基金 +4 bud_bud 2026-08-27 7/350 2026-08-29 07:23 by foolishmani
[基金申请] 哪位高人中了,把查询到的截图贴出来让我看看,让我长长见识 +5 yuleib84 2026-08-26 6/300 2026-08-28 00:02 by yudaoqian88
[文学芳草园] 梦想 +3 myrtle 2026-08-26 3/150 2026-08-27 10:01 by angelyueyi
[基金申请] 为什么国自然不能直接公布 +4 bjdxyxy 2026-08-26 4/200 2026-08-26 13:12 by qingmu1201
[基金申请] 国合现在查不到了吗? +10 chengyan1220 2026-08-24 20/1000 2026-08-26 08:57 by peasantsprig
[基金申请] 某些机构,以效率低为荣,以效率低作为存在感 +9 yuleib84 2026-08-25 10/500 2026-08-25 17:14 by alexon
[基金申请] 如果此刻你正在为国基感到焦虑,不妨来听听这首《基金之外》 +8 scalable 2026-08-24 8/400 2026-08-25 12:52 by jnhyjjm
信息提示
请填处理意见