现代统计图形 - 科学网—博客

5.31 生存函数图 131
1 > library(survival)
2 > leukemia.surv = survfit(Surv(time, status) ~ x, data = aml)
3 > plot(leukemia.surv, lty = 1:2, xlab = "time")
4 > legend("topright", c("Maintenance", "No Maintenance"),
5 + lty = 1:2, bty = "n")
0.0 0.2 0.4 0.6 0.8 1.0
0 50 100 150
time
图 5.37: 急性髓细胞白血病病人生存函数图
Maintenance
No Maintenance
注意图5.36的QQ图中我们用了纵横比参数asp = 1,这样使得图中横
纵坐标的单位长度表示的数值大小相同,为我们比较点的横纵坐标值大小
提供了更好的视觉辅助。
QQ图的用途不仅在检查数据是否服从某种特定理论分布,它也可以
推广到检查数据是否来自某个位置参数分布族。例如,若数据来自正态分
布 5 N(5, 1),我们拿它和标准正态分布N(0, 1)画QQ图的话,数据点仍然会
大致排列在直线上,此时直线的斜率仍然是1,但截距就不是0了,而是5。
5.31 生存函数图
在很多医学研究中,我们主要关心的变量是病人的某种事件发生的时
间,例如死亡、疾病复发等。事实上,以“生存时间”为研究对象的领域
并不仅限于医学,例如在金融领域,我们可能需要了解信用卡持有者的
5 正态分布是位置参数分布族中的一种分布:若X ∼ N(µ, σ 2 ),那么X + δ ∼ N(µ + δ, σ 2 )仍然是
正态分布。均值µ是位置参数。