生存分析预处理

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

对不起各位，原数据有误，已更新，非常抱歉！12小时分值变化的生存分析，v0代表起始分值，v1-12为每个小时的分值。规定：①分值降低且不再反弹为降低；②分值变为0且不再反弹为消失。数据类型如下：

obs	v0	v1	v2	v3	v4	v5	v6	v7	v8	v9	v10	v11	v12
1	10	9	4	5	2	0	2	1	0	0	0	0	0
2	4	9	4	1	2	0	2	1	0	1	2	3	3
3	10	3	4	2	1	0	2	1	0	0	0	0	0
4	4	9	4	5	1	0	2	1	0	0	4	4	4
5	5	2	0	2	1	0	0	0	0	0	0	0	0
……

第1条观测 降低为6小时（6小时为2，之后没有出现反弹），消失为8小时；
第2条观测 降低为11小时（11小时为3，之后没有出现反弹），消失为12小时（因在12结束时没有发生终点事件“消失”，所以定义为 截尾数据）；
第3条观测 降低为6小时，消失为8小时；
第4条观测 降低为12小时（因在12结束时没有发生终点事件“降低”，所以定义为截尾数据），消失为12小时（因在12结束时没有发生终点事件“消失”， 截尾数据）；
第5条观测 降低为3小时，消失为5小时；
……
观测数非常多，求指点程序怎么写能得到“降低”和“消失”小时数字，以及是否截尾，谢谢！！！

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏0 回帖

关键词：生存分析 预处理 截尾数据 数据类型 求指点

请把数据文件上传到论坛。

第1条观测的降低为什么是3小时？

pobel 发表于 2013-1-23 16:18
第1条观测的降低为什么是3小时？

因为5之后没有反弹，之前的有反弹，不算降低，游戏规则是这么定的。。。

ccgcs 发表于 2013-1-23 17:00
因为5之后没有反弹，之前的有反弹，不算降低，游戏规则是这么定的。。。

还是看不懂降低这个变量。
第一条观测v5到v6不算反弹？
第二条观测的降低为11小时，为什么第四条观测的降低为12小时，而不是10小时？
第五条观测的降低为什么不是3小时？

楼主，解析一下每条观测为什么降低是这么多
不理解降低是咋算的
不清楚，没办法编程。

pobel 发表于 2013-1-23 17:05
还是看不懂降低这个变量。
第一条观测v5到v6不算反弹？
第二条观测的降低为11小时，为什么第四条观测的 ...

对不起，原数据有误，已更新，谢谢！

楼主可以试试这个程序：

data have;
        input obs v0 v1 v2 v3 v4 v5 v6 v7 v8 v9 v10 v11 v12;
        cards;
1  10 9 4 5 2 0  2 1 0 0 0 0 0
2 4 9 4 1 2 0 2 1 0 1 2 3 3
3  10 3 4 2 1 0 2 1 0 0 0 0 0
4 4 9 4 5 1 0 2 1 0 0 4 4 4
5 5 2 0 2 1 0 0 0 0 0 0 0 0
;

data wanted;
    set have;
        censor=(v12 ne 0);
        array v(13) v0-v12;
        do i=13 to 2 by -1;
        if v(i-1) >=v(i) then remission=i-2;
                if v(i)=0 then resolve=i-1;
                if v(i-1)< v(i) then i=0;
        end;
        if v12>=V0 then do; remission=12; resolve=12; end;
        if censor=1 then resolve=12;
        drop i;
run;