censored data的威布尔三参数估计（RRX)

qinyu

辛苦了，thanksalot！
But,你的公式是不是写错了？怎么算出来和列出来的不一样？
Ps,我是为了探讨Weibull3P的RRX方法，数据只是来举例的。

[本帖最后由qinyu于2009-7-2410:17编辑]

Jack315

MedianRank的计算见下图：

下表是用在线秩计算工具(http://www.chinarel.com/onlineTool/tool2.html)得到的结果供参考：

Order	MedianRank
1	0.0943
2	0.2285
3	0.3641
4.25	0.5340
6.125	0.7884

==========================================================================
2009-7-26
最后给出了Weibull++的秩计算结果，以验证计算方法。

[本帖最后由Jack315于2009-7-2602:37编辑]

zhangxw

上面的计算看着不是很明白，我提取威布尔三参数是用商用软件，即便对于censoreddata

Jack315

(1)DistributionWizard结果表明应优先采用Weibull3P分布。
(2)估计出的参数如下表所示：

参数	值
Beta	2.2685
Eta	26.6618
Gamma	17.3800
Rho	0.9987

(3)分析方法：RRX(RankRegressiononX);SRM(StandardRegressionMethod)

附件给出了Weibull++文档供参考。

没时间玩了，最后的回归方法留待大侠出手了。楼下继续...

[本帖最后由Jack315于2009-7-2602:34编辑]

qinyu

在用RRX进行截尾数据的威布尔三参数估计过程中，mediarank法F(x)=（i-0.5）/(n+0.4)的i计算方法有点麻烦，如上面几楼的过程。当然也可以用专业的计算软件，但是老板都是想让我们来做出这个程序，这样就可以既省钱又省力了。我研究了两天，还是没参透编程实现方法，也许自己水平太有限吧。

有没有人有兴趣一起来研究一下或者指教一下，非常感谢！

[本帖最后由qinyu于2009-7-3011:00编辑]

Jack315

原帖由qinyu于2009-7-3010:52发表
在用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编辑]

Jack315

重开了个新贴FormularofOrderNumberforMedianRankCalculationwithCensoredData来讨论Order公式的问题。

qinyu

谢谢大侠fuckit指点。我已经解决了这个问题，因为该方法用程序实现会很方便，故我用VBA实现了RRX的威布尔三参数估计。以供参考。

另：我算出的三参貌似正确的，但是和Weibull++并不是完全一样，不能保证该方法完全正确。

[本帖最后由qinyu于2009-8-313:11编辑]

nanguachi2009

很好好好学习

deadxiaoh

不错，好好学习一下威布尔分布