文档帮助中心文档
计算Cox比例风险模型的生存
s =生存(coxMdl)
s =生存(coxMdl X)
s =生存(coxMdl X,分层)
s =生存(___“时间”,T)
(年代,兜售)=生存(___)
年代=生存(coxMdl)估计Cox比例风险模型的基线生存函数coxMdl.生存函数t估计存活到时间的概率是多少t.这个词基线指在预测指标确定基线时的生存函数。该值存储在coxMdl。基线,默认值为用于训练的数据集的平均值。
年代=生存(coxMdl)
年代
coxMdl
t
coxMdl。基线
例子
年代=生存(coxMdl,X)当预测因子有值时,估计生存函数X.在这种情况下,年代每一行对应的列是X.
年代=生存(coxMdl,X)
X
年代=生存(coxMdl,X,分层)对给定的分层变量的生存函数进行估计分层.你必须有一行分层的每一行X.
年代=生存(coxMdl,X,分层)
分层
请注意
当你训练coxMdl使用分层变量和通过预测变量X,生存还需要你通过分层变量。
生存
年代=生存(___“时间”,T)有时计算生存T使用前面语法中的任何输入参数组合。
年代=生存(___“时间”,T)
T
[年代,吹捧) =生存(___)还有《泰晤士报》吹捧每一个生存估计都是在这里计算的。
[年代,吹捧) =生存(___)
吹捧
全部折叠
执行一个Cox比例风险回归灯泡数据集,其中包含灯泡的模拟寿命。灯泡数据的第一列包含两种不同类型灯泡的生存期(以小时为单位)。第二列包含一个二元变量,该二元变量指示灯泡是荧光灯还是白炽灯;0表示灯泡为荧光灯,1表示灯泡为白炽灯。第三列为截尾信息,其中0表示观察到灯泡失效,1表示观察到截尾。
灯泡
拟合灯泡寿命的Cox比例危险模型,考虑了审查。预测变量是灯泡的类型。
负载灯泡coxMdl = fitcox(灯泡(:,2),灯泡(:1),...“审查”灯泡(:3));
计算基线生存函数作为时间的函数t,意思是灯泡在一段时间后失效的概率t.默认情况下,基线是为预测器的平均值计算的,在本例中就是这样意思是(灯泡(:,2))=0.5.返回时间吹捧求出存活函数。
意思是(灯泡(:,2))
0.5
[年代,兜售]=生存(coxMdl);
用时间的阶梯图来绘制生存。(《纽约时报》吹捧也在coxMdl.Hazard (: 1).)
coxMdl.Hazard (: 1)
持有在;楼梯(吹捧,年代,“b -”)包含“时间\ t”ylabel“时间过后的失败概率”
用荧光灯和白炽灯的生存功能覆盖情节。
s_fluorescent =生存(coxMdl, 0);s_incandescent =生存(coxMdl, 1);楼梯(s_fluorescent吹捧,的r -)楼梯(s_incandescent吹捧,“k -”)传说(“基线”,“荧光”,“白炽灯”)举行从
要在不创建生存数据的情况下创建情节,请使用plotSurvival.
plotSurvival
加载coxModel数据。(这个模拟数据是在示例中生成的Cox比例风险模型对象.)模型的命名coxMdl是否有三个分层水平(1、2和3)和一个预测因子X有三个分类值(1,1 /20和1/100)。
coxModel
负载coxModel
求生存函数X = 1在三个层次上。
X = 1
c1 =分类(1);X = (c1, c1, c1);分层= [1,2,3];s =生存(coxMdl X,分层);
绘制三个生存函数。首先,找出三个层次的时间。
t1 = find(coxmdd . hazard (:,3) == 1);t1 = coxMdl.Hazard (t1, 1);t2 = find(coxmdd . hazard (:,3) == 2);t2 = coxMdl.Hazard (t2, 1);t3 = find(coxmdd . hazard (:,3) == 3);t3 = coxMdl.Hazard (t3, 1);
绘制三个关卡的生存情况。查看从1到30的图。
情节(t1, {1}, t2, {2}, t3, {3}) xlim([1, 30])传说(“分层级别1”,“分层级别2”,“分层三级”,“位置”,“东北”)包含(“t”) ylabel (“超过t存活的概率”)
或者,通过指定时间论点。
时间
t = linspace (30300);圣=生存(coxMdl X,分层,“时间”t);图绘制(圣{1}t, t,圣{2}t, st{3})传说(“分层级别1”,“分层级别2”,“分层三级”,“位置”,“东北”)包含(“t”) ylabel (“超过t存活的概率”)
CoxModel
拟合的Cox比例风险模型,指定为CoxModel对象。创建coxMdl使用fitcox.
fitcox
0
模型的预测器,指定为用于训练的相同类型的预测器数组coxMdl.每一行的X表示一组预测器。
数据类型:双|表格|分类
双
表格
分类
分层水平,指定为用于培训的变量或相同类型的变量coxMdl.中指定相同的行数分层就像在X.
数据类型:单|双|逻辑|字符|字符串|表格|细胞|分类
单
逻辑
字符
字符串
细胞
生存估计时间,指定为真实向量。生存对指定的时间进行排序,并在必要时将其转换为列向量。结果值由训练数据中的时间线性插值。
例子:0:40
0:40
数据类型:双
生存估计,作为数字列向量或数字列向量的单元格数组返回。
对于非分层模型,年代是估计概率的排序数值列向量。
对于分层模型,年代是每个层的估计概率的排序数值列向量的单元数组。
生存的每一行返回生存估计的列X.
mdl.Hazard (: 1)
生存估计时间,返回如下。
对于非分层模型,吹捧为训练集中次数排序的数值列向量。
对于分层模型,吹捧为每个层次的训练集中训练次数的排序数值列向量的单元数组。
的coxMdl.Hazard (: 1)向量包含分层模型和非分层模型的时间。对于分层模型,不同分层级别的时间由0项分隔。
数据类型:双|细胞
CoxModel|fitcox|hazardratio|plotSurvival
hazardratio
您有这个示例的修改版本。您想打开这个示例与您的编辑吗?
你点击一个链接对应于这个MATLAB命令:
通过在MATLAB命令窗口中输入命令来运行命令。Web浏览器不支持MATLAB命令。金宝app
选择一个网站,在那里获得翻译的内容,并看到当地的活动和优惠。根据您的位置,我们建议您选择:.
你也可以从以下列表中选择一个网站:
选择中国网站(中文或英文)以获得最佳网站性能。其他MathWorks国家站点没有针对您所在位置的访问进行优化。
与当地办事处联系