主要内容
生存
Cox比例风险模型的计算生存
描述
例子
计算生存
执行一个Cox比例风险回归的灯泡数据集,其中包含模拟灯泡的寿命。灯泡的第一列数据包含的寿命(小时)两种不同类型的灯泡。第二列包含一个二进制变量指示是否荧光灯或白炽灯泡;0表示荧光灯泡,白炽灯和1表示。第三列包含审查信息,0表示观察灯泡,直到失败,和1表明观察被审查。
适合Cox比例风险模型的生命周期灯泡,占审查。预测变量是灯泡的类型。
负载灯泡coxMdl = fitcox(灯泡(:,2),灯泡(:1),…“审查”灯泡(:3));
计算基线生存函数作为时间的函数t后,这意味着一个灯泡失败的概率t。默认情况下,基线计算均值的预测,在这种情况下意思是(灯泡(:,2))=0.5。返回时间吹捧的生存函数计算。
[年代,兜售]=生存(coxMdl);
情节的生存作为一个梯级的图。(《纽约时报》吹捧也在coxMdl.Hazard (: 1))。
持有在;楼梯(吹捧,年代,“b -”)包含“时间\ t”ylabel失败的概率又一次\ t '
覆盖情节的生存函数荧光和白炽灯泡。
s_fluorescent =生存(coxMdl, 0);s_incandescent =生存(coxMdl, 1);楼梯(s_fluorescent吹捧,的r -)楼梯(s_incandescent吹捧,“k -”)传说(“基线”,“荧光”,“白炽灯”)举行从
创建块没有首先创造生存数据,使用plotSurvival。
生存对于分层模型
加载coxModel数据。(模拟数据生成的例子Cox比例风险模型对象)。模型的命名coxMdl有三个分层水平(1、2和3)和预测X有三个分类值(1、1/20和1/100)。
负载coxModel
计算的生存函数X = 1在这三个分层水平。
c1 =分类(1);X = (c1, c1, c1);分层= [1,2,3];s =生存(coxMdl X,分层);
情节三个生存函数。首先,找到三个分层次的水平。
t1 =找到(coxMdl.Hazard (:, 3) = = 1);t1 = coxMdl.Hazard (t1, 1);t2 =找到(coxMdl.Hazard (:, 3) = = 2);t2 = coxMdl.Hazard (t2, 1);t3 =找到(coxMdl.Hazard (:, 3) = = 3);t3 = coxMdl.Hazard (t3, 1);
三个层次的生存。视图的阴谋通过30 * 1。
情节(t1, {1}, t2, {2}, t3, {3}) xlim([1, 30])传说(“分层级别1”,“分层级别2”,“分层三级”,“位置”,“东北”)包含(“t”)ylabel (“过去的生存概率t”)
另外,评估的生存* 1通过指定30时间论点。
t = linspace (30300);圣=生存(coxMdl X,分层,“时间”t);图绘制(圣{1}t, t,圣{2}t, st{3})传说(“分层级别1”,“分层级别2”,“分层三级”,“位置”,“东北”)包含(“t”)ylabel (“过去的生存概率t”)
输入参数
输出参数
介绍了R2021a
