回复 7楼 fuckit 的帖子
辛苦了,thanksalot!But,你的公式是不是写错了?怎么算出来和列出来的不一样?
Ps,我是为了探讨Weibull3P的RRX方法,数据只是来举例的。
[本帖最后由qinyu于2009-7-2410:17编辑]
Median Rank
MedianRank的计算见下图:下表是用在线秩计算工具(http://www.chinarel.com/onlineTool/tool2.html)得到的结果供参考:
OrderMedianRank10.094320.228530.36414.250.53406.1250.7884
==========================================================================
2009-7-26
最后给出了Weibull++的秩计算结果,以验证计算方法。
[本帖最后由Jack315于2009-7-2602:37编辑] 上面的计算看着不是很明白,我提取威布尔三参数是用商用软件,即便对于censoreddata
Weibull 3P结果
(1)DistributionWizard结果表明应优先采用Weibull3P分布。(2)估计出的参数如下表所示:
参数值Beta2.2685Eta26.6618Gamma17.3800Rho0.9987
(3)分析方法:RRX(RankRegressiononX);SRM(StandardRegressionMethod)
附件给出了Weibull++文档供参考。
没时间玩了,最后的回归方法留待大侠出手了。楼下继续...
[本帖最后由Jack315于2009-7-2602:34编辑]
请教:censored data 的中位秩中order i 的计算方法(通式)。
在用RRX进行截尾数据的威布尔三参数估计过程中,mediarank法F(x)=(i-0.5)/(n+0.4)的i计算方法有点麻烦,如上面几楼的过程。当然也可以用专业的计算软件,但是老板都是想让我们来做出这个程序,这样就可以既省钱又省力了。我研究了两天,还是没参透编程实现方法,也许自己水平太有限吧。有没有人有兴趣一起来研究一下或者指教一下,非常感谢!
[本帖最后由qinyu于2009-7-3011:00编辑] 原帖由qinyu于2009-7-3010:52发表http://www.kekaoxing.com/club/images/common/back.gif
在用RRX进行截尾数据的威布尔三参数估计过程中,mediarank法F(x)=(i-0.5)/(n+0.4)的i计算方法有点麻烦,如上面几楼的过程。当然也可以用专业的计算软件,但是老板都是想让我们来做出这个程序,这样就可以既省钱又...
Order的计算公式确实有点难...。
若从实用角度出发的话,建议LZ不要用RRX(RankRegressiononX/Y),因为当截尾样本数增多时,模型精度较差。可改用极大似然法(MLE-MaximumLikelihood)来算--模型精度高(因为截尾样本也参与计算),且计算较RRX难度低(用Excel可以搞定)。
具体公式请参看ALTA帮助文件:AppendixB:ParameterEstimation->MLEofAcceleratedLifeData
[本帖最后由Jack315于2009-8-302:07编辑] 重开了个新贴FormularofOrderNumberforMedianRankCalculationwithCensoredData来讨论Order公式的问题。
VBA实现的威布尔三参数估计方法(RRX)
谢谢大侠fuckit指点。我已经解决了这个问题,因为该方法用程序实现会很方便,故我用VBA实现了RRX的威布尔三参数估计。以供参考。另:我算出的三参貌似正确的,但是和Weibull++并不是完全一样,不能保证该方法完全正确。
[本帖最后由qinyu于2009-8-313:11编辑] 很好好好学习 不错,好好学习一下威布尔分布