%输入参数:
% *:一个n×2矩阵指的时间数据。两列的
%定义的间隔事件发生。如果在一行
%”时代”元素相等,那么这个对应于一个事件的时间。如果
%左列时间点是有限的和右列时间点是正无穷
%然后我们有权审查。如果时间点是负无穷
%的时间点是有限的我们审查。例如:
%
%
% 5正- >审查观察
%从8 - >左审查观察
% 8 11 - >间隔审查
% 9 9 - >完全(完全)
%
%注意,因为消极的时间值也是允许的,离开了审查
您需要指定%值左极限为负。需要澄清的是,正无穷
%无穷替换为“NA”在相应survfit R的函数。
%
%输出参数:
%时间:时间点生存的评估
% S:生存的估计
% SE:估计的标准误差
% L:相应的95%置信区间的下限
%的生存
% U:相应的95%置信区间的上限
%的生存
%
%参见示例文件。
%
%如果只有正确的审查存在那么通常的卡普兰迈耶是就业
%,没有必要使用这段代码从MATLAB已经建立
%的功能“ecdf”将更快。如果左/间隔审查
%的特恩布尔算法。这个算法也
%详细描述:
%
%”克莱因和Moeschberger(2003),生存分析(技术审查
%和截断数据),斯普林格出版社,纽约。”
%
%这个程序调用相应的R函数而已。
%有必要下载MATLAB R-link并遵循所有指令
%在描述:
%
%//www.tianjin-qmedu.com/matlabcentral/fileexchange/5051-matlab-r-link
%
%这意味着你还需要下载并安装Randfriends:
%
%http://rcom.univie.ac.at/download.html
%
%这适用于32位和64位的机器。
%还请运行演示包含在MATLAB R-link包,和
%确保一切工作之前这段代码。
%之后,您必须下载并安装包“生存”
%和matlab在r .这段代码实际上是提供的分析
%输出作为相应的生存提供的包(survfit)
%由特里Therneau教授提供。
%
%
%
%的引用:
%
%特恩布尔,b . w .双重生存函数的非参数估计
%审查数据。美国统计协会杂志》69 (1974):169 - 173。
%
%特恩布尔,b . w .经验分布函数和任意分组,审查
%和截断数据。英国皇家统计学会杂志38 (1976):290 - 295。
%
%克莱因和Moeschberger(2003),生存分析(技术审查
%和截断数据),斯普林格出版社,纽约。”
%
特里Therneau % R生存方案。参见:
%http://cran.r-project.org/web/packages/survival/survival.pdf
%
%