R plot.cox.zph 比例风险的图形测试

说明

显示缩放舍恩菲尔德残差的图形以及平滑曲线。

用法

## S3 method for class 'cox.zph'
plot(x, resid=TRUE, se=TRUE, df=4, nsmo=40, var,
        xlab="Time", ylab, lty=1:2, col=1, lwd=1, hr=FALSE, ...)

参数

x	cox.zph 函数的结果。
resid	如果 TRUE 残差包含在绘图中，则为逻辑值以及平滑拟合。
se	逻辑值，如果 TRUE ，将添加两个标准误差的置信带。
df	拟合自然样条的自由度 df=2 导致线性拟合。
nsmo	用于线的点数
var	需要绘图的变量集。默认情况下，依次为模型的每个变量生成绘图。选择单个变量允许将其他特征添加到图中，例如零处的水平线或主标题。 这已被下标方法取代；请参见下面的示例。
hr	如果为 TRUE，则使用估计的风险比而不是估计的系数来标记 y 轴。 (绘图没有改变，只有轴标签改变。)
xlab	绘图 x 轴的标签
ylab	绘图 y 轴的可选标签。如果缺少，则提供默认标签。这可以是标签向量。
lty , col , lwd	重叠曲线的线型、颜色和线宽。其中每个元素都可以是长度为 2 的向量，在这种情况下，第二个元素用于置信区间。
...	传递给 plot 函数的附加图形参数。

副作用

在当前图形设备上生成绘图。

例子

vfit <- coxph(Surv(time,status) ~ trt + factor(celltype) + 
              karno + age, data=veteran, x=TRUE) 
temp <- cox.zph(vfit) 
plot(temp, var=3)      # Look at Karnofsy score, old way of doing plot 
plot(temp[3])     # New way with subscripting 
abline(0, 0, lty=3) 
# Add the linear fit as well  
abline(lm(temp$y[,3] ~ temp$x)$coefficients, lty=4, col=3)  
title(main="VA Lung Study") 

也可以看看

coxph、cox.zph。