SurvivalModelFit
SurvivalModelFit[{e1,e2,…}]
创建一个对应于事件时间 ei 的生存模型.
更多信息和选项
- SurvivalModelFit 用于生存、可靠性和耐久性分析. 给出非参数生存分布以及直接来自于不完全数据的置信测量.
- SurvivalModelFit 返回一个符号式 SurvivalModel 对象表示它所构建的生存模式. 可以从 model["property"] 获得模型的属性和诊断.
- 事件 ei 的形式遵循 EventData 中使用的形式.
- 通过指定 model[h] 中的形式 h 获取生存模式的不同函数形式. 可以用以下形式:
-
"CDF" 累积分布函数 "CHF" 累积风险函数 "SF" 生存函数 - 在特殊时间 t 来自 SurvivalModelFit 的拟合函数 h 的值可以从 model[h][t] 中求得. model[t] 等同于 model["SF"][t].
- 指定 Normal[model] 给出估计生存函数的纯函数形式.
- 可用模型属性的列表可用 model["Properties"] 获取.
- 模型拟合函数 h 的属性可用 model["property",h] 获取. 表达式 model["property"] 等同于 model["property","SF"].
- 与拟合函数相关的某些属性可用 model["property"] 获取:
-
"EventTable" 总结拟合函数的格式化表格 "EventTableEntries" 事件表格的项 "FittedFunction" 拟合函数的纯函数形式 "FullEventTable" 扩展的事件表格 "FullEventTableEntries" 完全事件表格的项 "MeanSurvival" 平均生存时间 "MedianSurvival" 中位数生存时间 "StandardErrorFunction" 拟合函数的标准误差的函数形式 "StandardErrors" 每个时间点的标准误差列表 "SurvivalDistribution" 拟合中使用的 DataDistribution 对象 - 拟合函数的置信区间是由对应于由 model["EstimationPoints"] 给出的每个时间点的成对值 {{h1,min,h1,max},…} 列表给出.
- 以纯函数形式给出置信带,等同于在 model["EstimationPoints"] 中所有时间点计算的置信区间.
- 与置信区间和置信带相关的属性包括:
-
"PointwiseBands" 逐点置信带 "PointwiseIntervals" 逐点置信区间 "EqualPrecisionBands" 同步置信带 "EqualPrecisionIntervals" 同步置信间隔 "HallWellnerBands" 同步 Hall–Wellner置信带 "HallWellnerIntervals" 同步 Hall–Wellner置信间隔 - 与事件数据相关的属性:
-
"EstimationPoints" 估计中使用的时间点 "EventData" 估计中使用的 EventData 对象 "EventMatrix" 表明可能的事件位置的矩阵 "EventMatrixPlot" 事件矩阵自定义的矩阵图 "TruncationMatrix" 表明在哪单位可观察的矩阵 "TruncationMatrixPlot" 截断矩阵的自定义矩阵图 - 与观察的不同类型计数相关的属性:
-
"AtRisk" 每个时间点的高危个体数目 "EventCounts" 每个时间点的事件数目 "LeftCensoringCounts" 每个时间点的左删失事件数目 "ObservationCount" 试验单位的有效数目 "RightCensoringCounts" 每个时间点的右删失事件数目 - SurvivalModelFit 接受以下选项:
-
ConfidenceLevel 95/100 用于区间和带的级别 ConfidenceRange All 同步置信带的范围 ConfidenceTransform "LogLog" 使用的置信变换 Method Automatic 用于模型拟合的方法 - 当 ConfidenceLevel->p 时,概率 p 置信区间和置信带用于计算各种函数形式.
- ConfidenceRange->{tmin,tmax} 为 tmin 与 tmax 之间的拟合函数给出概率 p 同步置信区间和置信带.
- ConfidenceTransform 的可能设置包括 "Linear"、"LogLog"、"ArcSinSqrt"、"Log"、"Logit" 或一个纯函数 g.
- SurvivalModelFit 自动选择最适合数据的方法. Kaplan–Meier 估计器被用于右删失数据. 对于其他删失类型,使用自我一致性构建估计. 不同的方法可能只支持某些类型的删失或截断.
- 设置 Method->estimator 指定用于分布函数的估计器. 可能的设置包括:
-
"KaplanMeier" 右删失数据的乘积极限估计器 "NelsonAalen" 基于 Nelson–Aalen的累积风险估计器 "Turnbull" 区间删失数据的 Turnbull 算法 "Noncensored" 删失被忽略,使用区间中点 "SelfConsistency" 用于双删失数据的 Turnbull 算法 - 其他方法设置可以在范例中的选项部分找到. »
范例
打开所有单元关闭所有单元范围 (29)
基本用途 (4)
模型属性 (16)
在每个时间点的关于生存函数的 95% Hall–Wellner 区间:
对估计的 SurvivalFunction 的总结制表:
时间介于 25 和 35 之间,步长为 5 的具有事件历史数据的事件表格:
数据类型 (9)
删失 (7)
使用具有状态指示器的 EventData 表明删失观察:
选项 (10)
ConfidenceLevel (3)
ConfidenceTransform (3)
Method (3)
如果上次观察是删失的,超过那点的生存会是 Indeterminate:
子选项 "StandardErrorMethod" 可用于控制标准误差计算:
设置 Automatic 通过插入方法尝试计算标准误差:
应用 (8)
对一组26个精神病人进行跟踪调查直至死亡或者研究结束. 在研究结束后这些病人的寿命呈右删失分布. 比较病人存活的非参数和参数模型:
在 95% 每一点的置信带和 Gompertz–Makeham模型下,Kaplan–Meier估值:
下面数据包含70台电子发电机的连续运转时间(以小时为单位). 在给定最低限度的维修策略下,感兴趣的是各种运转时间下期望的平均失效数目. 有些发电机在观察器没有失效,并且是右删失的:
平均而言,在 87,500 个小时的运行时间内,失败了 0.337 次:
为缓解白血病人的痛苦,在临床试验中使用 6-巯基嘌呤或安慰剂. 比较两组的存活率:
具有 Tarone–Ware 加权的LogRankTest 与上面的发现是一致的:
191 名高中男孩被问及,他们第一次尝试大麻是什么年龄. 可能的回答包括,“我第一次尝试大麻是在年龄 ”,“我从没有抽过大麻”,或“我不记得何时第一次尝试大麻”:
一组38例急性淋巴细胞白血病的个人分别给予骨髓移植. 患者随访直至死亡,复发,或研究结束. 下面的数据可以被解释为无病存活时间:
感兴趣研究狒狒睡眠习惯的研究人员记录了有一半狒狒睡眠后从树上掉下来的所经过的时间. 在某些日子里,研究人员到达的时间足够早,可以精确记录这些时间. 在其他时候,他们在下降时间的中位数时间后到达,因此一天的观测被左删截了:
在7到9个小时时间段内,95% Hall–Wellner生存置信包(survival confidence envelope):
跟踪肾功能不全患者到第一导管出口处感染. 每个患者有两种类型的导管放置(手术或经皮). 最令人感兴趣的是导管的位置是最重要的,可以防止感染:
诊所招收了30个糖尿名患者并随访至死亡或研究结束. 并记录了患者在研究开始和结束时的年龄. 使用这些假设数据,比较60和65岁糖尿病患者的条件存活率和糖尿病存活率,忽略左截断:
属性和关系 (7)
对于右删失数据,默认情况下使用 Kaplan–Meier 估计器:
使用 CoxModelFit 估计协变量出现的生存概率:
使用 LogRankTest 比较多组的风险率:
使用 LogRankTest 比较风险率:
文本
Wolfram Research (2012),SurvivalModelFit,Wolfram 语言函数,https://reference.wolfram.com/language/ref/SurvivalModelFit.html.
CMS
Wolfram 语言. 2012. "SurvivalModelFit." Wolfram 语言与系统参考资料中心. Wolfram Research. https://reference.wolfram.com/language/ref/SurvivalModelFit.html.
APA
Wolfram 语言. (2012). SurvivalModelFit. Wolfram 语言与系统参考资料中心. 追溯自 https://reference.wolfram.com/language/ref/SurvivalModelFit.html 年