你也许疑惑在某些SAS提供的生存分析过程里,比如说PROC LIFEREG, PROC RELIABILITY (QCpackage), 如何绘制和阅读概率分布图.基本上讲,该图通过比较实际概率分布和理论概率分布,进而得出拟和的相似程度,以判断模型的优劣.所谓的理论分布是你给出的参数模型在某些应变量条件上的取值.具体我将通过我的一片小文进行阐述.

In the following examplethe data is from Example 50.5Probability Plotting—Right Censoring.

      SAS提供不同的方法来估计实际的分布函数,但是方法和结果大同小异. 以下展示的方法为Expected Rank Plotting Position(mi).而对mi进行变换将得到PPLOT的Y轴.X的取值可以是时间或者log(时间).而与其对比的理论分布直线来源于模型分布本身.对于伪不分布,Y =log(1/(1-cdf)).

      cdf是分布参数的函树.其中更加详细的计算我将在我blog里给出.如果计算合适,得出的结果应该和SAS给出的吻合一致.

复制代码

京剧