求问如何把0-1面板变量变成生存分析变量？

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

想要计算1的持续期，如何转换呢？
似乎要进行KM分析需要一个截面数据。
一般假设是各持续段是不相关的。

数据如下：

1. * Example generated by -dataex-. To install: ssc install dataex
   
2. clear
   
3. input int year float down
   
4. 2001 0
   
5. 2002 1
   
6. 2003 0
   
7. 2004 0
   
8. 2005 1
   
9. 2006 1
   
10. 2007 1
    
11. 2008 0
    
12. 2009 1
    
13. 2010 0
    
14. 2011 0
    
15. 2012 1
    
16. 2013 1
    
17. 2014 0
    
18. 2015 1
    
19. 2016 1
    
20. 2017 0
    
21. 2018 1
    
22. end

复制代码

扫码加我 拉你入群

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

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

分享0 收藏2 回帖

关键词：生存分析 generated generate install example

按照逻辑，最终数据格式应该如下：


* Example generated by -dataex-. To install: ssc install dataex
clear
input int year float down byte downmax
2001 0 .
2002 1 1
2003 0 .
2004 0 .
2005 1 .
2006 1 .
2007 1 3
2008 0 .
2009 1 1
2010 0 .
2011 0 .
2012 1 .
2013 1 2
2014 0 .
2015 1 .
2016 1 2
2017 0 .
2018 1 1
end

按照上述数据，样本期内有6个持续段，分别持续了2002（1年）、2005-2007（3年）、2009（1年）、2012-2013（2年）、2015-2016（2年）、2018（1年）。

你要的是这种效果吗？

1. * Example generated by -dataex-. To install: ssc install dataex
   
2. clear
   
3. input int year float down
   
4. 2001 0
   
5. 2002 1
   
6. 2003 0
   
7. 2004 0
   
8. 2005 1
   
9. 2006 1
   
10. 2007 1
    
11. 2008 0
    
12. 2009 1
    
13. 2010 0
    
14. 2011 0
    
15. 2012 1
    
16. 2013 1
    
17. 2014 0
    
18. 2015 1
    
19. 2016 1
    
20. 2017 0
    
21. 2018 1
    
22. end
    
23. 
    
24. spell down,cond(down==1)
    
25. bys  _spell  :egen wanted=max(_seq*_end)
    
26. sort year down
    
27. 
    
28. . list
    
29. 
    
30.      +---------------------------------------------+
    
31.      | year   down   _seq   _spell   _end   wanted |
    
32.      |---------------------------------------------|
    
33.   1. | 2001      0      0        0      0        0 |
    
34.   2. | 2002      1      1        1      1        1 |
    
35.   3. | 2003      0      0        0      0        0 |
    
36.   4. | 2004      0      0        0      0        0 |
    
37.   5. | 2005      1      1        2      0        3 |
    
38.      |---------------------------------------------|
    
39.   6. | 2006      1      2        2      0        3 |
    
40.   7. | 2007      1      3        2      1        3 |
    
41.   8. | 2008      0      0        0      0        0 |
    
42.   9. | 2009      1      1        3      1        1 |
    
43. 10. | 2010      0      0        0      0        0 |
    
44.      |---------------------------------------------|
    
45. 11. | 2011      0      0        0      0        0 |
    
46. 12. | 2012      1      1        4      0        2 |
    
47. 13. | 2013      1      2        4      1        2 |
    
48. 14. | 2014      0      0        0      0        0 |
    
49. 15. | 2015      1      1        5      0        2 |
    
50.      |---------------------------------------------|
    
51. 16. | 2016      1      2        5      1        2 |
    
52. 17. | 2017      0      0        0      0        0 |
    
53. 18. | 2018      1      1        6      1        1 |
    
54.      +---------------------------------------------+

复制代码

zdlspace 发表于 2021-1-15 16:28
你要的是这种效果吗？

谢谢，只需要在您的基础上加上一行代码就是最终想要的了！万分感谢！

你要的那个更简单一点：

1. *或者
   
2. clear
   
3. input int year float down
   
4. 2001 0
   
5. 2002 1
   
6. 2003 0
   
7. 2004 0
   
8. 2005 1
   
9. 2006 1
   
10. 2007 1
    
11. 2008 0
    
12. 2009 1
    
13. 2010 0
    
14. 2011 0
    
15. 2012 1
    
16. 2013 1
    
17. 2014 0
    
18. 2015 1
    
19. 2016 1
    
20. 2017 0
    
21. 2018 1
    
22. end
    
23. 
    
24. spell down,cond(down==1)
    
25. gen wanted2=_seq*_end
    
26. mvdecode wanted2,mv(0)      //或者replace wanted2=. if wanted2==0,我只是为了多用mvencode才用的，看上去高大上一点，哈哈

复制代码

zdlspace 发表于 2021-1-15 16:40
你要的那个更简单一点：

感谢！太赞了！

我是工作狂2009 发表于 2021-1-15 16:41
感谢！太赞了！

一般很多人要求更高一点，需要第一种效果，比较麻烦一点。不过你预期的效果好做多了

zdlspace 发表于 2021-1-15 16:43
一般很多人要求更高一点，需要第一种效果，比较麻烦一点。不过你预期的效果好做多了

第一次接触这个spell外部命令，真是强大的工具！

我是工作狂2009 发表于 2021-1-15 16:47
第一次接触这个spell外部命令，真是强大的工具！

多逛逛论坛，我们的论坛大神黄河泉老师经常推荐此命令