请教连老师，谢谢！！

1229

收藏 2011-10-26

连老师，重新发帖请教您，下面是数据情况：
hhid家庭编号； line 家庭成员编号；wave 年份；gender 性别；age 年龄；inc家庭人均收入； pk 家庭贫困状况（0=非贫困；1=贫困，是通过命令：gen pk=0. replace pk=1 if inc<=2492.1835计算出来的）

hhid	line	gender	wave	age	inc	pk
211101001	2	2	1989	36.54	3986.22	0
211101001	1	1	1989	37.54	3986.22	0
211101001	3	2	1989	10.05	3986.22	0
211101001	2	2	1991	38.48	3457.497	0
211101001	1	1	1991	39.49	3457.497	0
211101001	3	2	1991	11.99	3457.497	0
211101001	3	2	1993	13.89	4609.963	0
211101001	2	2	1993	40.39	4609.963	0
211101001	1	1	1993	41.39	4609.963	0
211101001	3	2	2000	20.8	7785.401	0
211101001	2	2	2000	47.29	7785.401	0
211101001	1	1	2000	48.29	7785.401	0
211101002	3	2	1989	9.78	3680.05	0
211101002	2	2	1989	36.93	3680.05	0
211101002	1	1	1989	37	3680.05	0
211101002	2	2	1991	38.87	2904.244	0
211101002	3	2	1991	11.72	2904.244	0
211101002	1	1	1991	38.93	2904.244	0
211101002	2	2	1993	40.78	3687.57	0
211101002	3	2	1993	13.63	3687.57	0
211101002	1	1	1993	40.85	3687.57	0
211101003	3	1	1989	13.05	4231.679	0
211101003	2	2	1989	38.73	4231.679	0
211101003	1	1	1989	39.95	4231.679	0
211101003	2	2	1991	40.67	4634.417	0
211101003	3	1	1991	14.99	4634.417	0
211101003	1	1	1991	41.89	4634.417	0
211101003	1	1	1993	43.79	4923.602	0
211101003	2	2	1993	42.58	4923.602	0
211101003	3	1	1993	16.9	4923.602	0
211101003	1	1	2000	50.69	2210.787	1
211101003	3	1	2000	23.8	2210.787	1
211101003	2	2	2000	49.48	2210.787	1
211101003	2	2	2004	53.55	1780.151	1
211101003	102	2	2006	6	4522.926	0
211101003	101	2	2006	31.91	4522.926	0
211101003	2	2	2009	58.64	5400	0
211101004	3	2	1989	11.49	4786.485	0
211101004	1	2	1989	35.88	4786.485	0
211101004	2	1	1989	38.43	4786.485	0
211101004	3	2	1991	13.43	8035.13	0
211101004	2	1	1991	40.36	8035.13	0
211101004	1	2	1991	37.82	8035.13	0
211101004	1	2	1993	39.73	3935.297	0

因为看了相关文献，想用cox回归来进行生存分析，所以需要得出时间变量--------即每个家庭从第一次被认定为贫困的年份到脱离贫困的年份之间的时间差所形成的变量，以及家庭第一次被认定为贫困的那一年的年龄、性别情况，现在很困扰的是怎么通过命令求出第一次被认定为贫困家庭的年份和脱离贫困的年份。。。（不好意思，不知道是否说清楚了。。？）

谢谢！！

扫码加我拉你入群

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

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

全部回复

susanketty

2011-10-30 11:12:17

呵呵，求关注啊！！

扫码加我拉你入群

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

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

arlionn

2011-11-1 11:19:02

*-Duration of poverty

*-https://bbs.pinggu.org/thread-1217063-1-1.html

*-说明：无法以家庭为单位进行分析，而要以家庭中的每个成员为单位进行分析

  cd D:\stata11\ado\personal\temp
  *insheet using data_pov_duration.txt, clear

*------------------input data----------------------------------
  clear
  input hhid line gender wave age inc pk
211101001 2 2 1989 36.54 3986.22 0
211101001 1 1 1989 37.54 3986.22 0
211101001 3 2 1989 10.05 3986.22 0
211101001 2 2 1991 38.48 3457.497 0
211101001 1 1 1991 39.49 3457.497 0
211101001 3 2 1991 11.99 3457.497 0
211101001 3 2 1993 13.89 4609.963 0
211101001 2 2 1993 40.39 4609.963 0
211101001 1 1 1993 41.39 4609.963 0
211101001 3 2 2000 20.8 7785.401 0
211101001 2 2 2000 47.29 7785.401 0
211101001 1 1 2000 48.29 7785.401 0
211101002 3 2 1989 9.78 3680.05 0
211101002 2 2 1989 36.93 3680.05 0
211101002 1 1 1989 37 3680.05 0
211101002 2 2 1991 38.87 2904.244 0
211101002 3 2 1991 11.72 2904.244 0
211101002 1 1 1991 38.93 2904.244 0
211101002 2 2 1993 40.78 3687.57 0
211101002 3 2 1993 13.63 3687.57 0
211101002 1 1 1993 40.85 3687.57 0
211101003 3 1 1989 13.05 4231.679 0
211101003 2 2 1989 38.73 4231.679 0
211101003 1 1 1989 39.95 4231.679 0
211101003 2 2 1991 40.67 4634.417 0
211101003 3 1 1991 14.99 4634.417 0
211101003 1 1 1991 41.89 4634.417 0
211101003 1 1 1993 43.79 4923.602 0
211101003 2 2 1993 42.58 4923.602 0
211101003 3 1 1993 16.9 4923.602 0
211101003 1 1 2000 50.69 2210.787 1
211101003 3 1 2000 23.8 2210.787 1
211101003 2 2 2000 49.48 2210.787 1
211101003 2 2 2004 53.55 1780.151 1
211101003 102 2 2006 6 4522.926 0
211101003 101 2 2006 31.91 4522.926 0
211101003 2 2 2009 58.64 5400 0
211101004 3 2 1989 11.49 4786.485 0
211101004 1 2 1989 35.88 4786.485 0
211101004 2 1 1989 38.43 4786.485 0
211101004 3 2 1991 13.43 8035.13 0
211101004 2 1 1991 40.36 8035.13 0
211101004 1 2 1991 37.82 8035.13 0
211101004 1 2 1993 39.73 3935.297 0
211101004 2 1 1993 42.27 3935.297 0
211101004 3 2 1993 15.34 3935.297 0
211101005 2 2 1989 42.43 7995.566 0
211101005 3 2 1989 13.71 7995.566 0
211101005 1 1 1989 43.78 7995.566 0
211101005 2 2 1991 44.37 3890.304 0
211101005 1 1 1991 45.73 3890.304 0
211101005 3 2 1991 15.65 3890.304 0
211101005 3 2 1993 17.56 3945.643 0
211101005 1 1 1993 47.63 3945.643 0
211101005 2 2 1993 46.28 3945.643 0
211101006 2 2 1989 55.29 4883.226 0
211101006 1 1 1989 55.33 4883.226 0
211101006 1 1 1993 59.19 10082.58 0
211101006 2 2 1993 59.15 10082.58 0
211101007 4 1989 2754.467 0
211101007 1 1 1989 38.45 2754.467 0
211101007 2 2 1989 37.03 2754.467 0
211101007 3 2 1989 11.02 2754.467 0
211101007 3 2 1991 12.97 2494.61 0
211101007 2 2 1991 38.97 2494.61 0
211101007 4 1 1991 5.53 2494.61 0
211101007 1 1 1991 40.39 2494.61 0
211101007 4 1 1993 7.44 3225.757 0
211101007 3 2 1993 14.87 3225.757 0
211101007 2 2 1993 40.88 3225.757 0
211101007 1 1 1993 42.3 3225.757 0
211101008 2 2 1989 48.5 4149.261 0
211101008 3 1 1989 24.45 4149.261 0
211101008 1 1 1989 50.16 4149.261 0
211101008 1 1 1989 50.16 4149.261 0
211101008 2 2 1989 48.5 4149.261 0
211101008 2 2 1989 48.5 4149.261 0
211101008 2 2 1989 48.5 4149.261 0
211101008 4 2 1989 26.78 4149.261 0
211101008 1 1 1989 50.16 4149.261 0
211101008 3 1 1989 24.45 4149.261 0
211101008 2 2 1989 48.5 4149.261 0
211101008 2 2 1989 48.5 4149.261 0
211101008 1 1 1989 50.16 4149.261 0
211101008 1 1 1989 50.16 4149.261 0
211101008 1 1 1989 50.16 4149.261 0
211101008 1 1 1989 50.16 4149.261 0
211101008 4 2 1989 26.78 4149.261 0
211101008 5 1 1989 2.11 4149.261 0
211101008 1 1 1989 50.16 4149.261 0
211101008 3 1 1989 24.45 4149.261 0
211101008 2 2 1989 48.5 4149.261 0
211101008 1 1 1989 50.16 4149.261 0
211101008 1 1 1989 50.16 4149.261 0
211101008 2 2 1989 48.5 4149.261 0
211101008 2 2 1989 48.5 4149.261 0
211101008 2 2 1989 48.5 4149.261 0
211101008 2 2 1991 50.44 4094.428 0
211101008 2 2 1991 50.44 4094.428 0
211101008 1 1 1991 52.1 4094.428 0
211101008 3 1 1991 26.39 4094.428 0
211101008 1 1 1991 52.1 4094.428 0
211101008 2 2 1991 50.44 4094.428 0
211101008 2 2 1991 50.44 4094.428 0
211101008 3 1 1991 26.39 4094.428 0
211101008 1 1 1991 52.1 4094.428 0
211101008 21 2 1991 21.64 4094.428 0
211101008 2 2 1991 50.44 4094.428 0
211101008 2 2 1991 50.44 4094.428 0
211101008 3 1 1991 26.39 4094.428 0
211101008 2 2 1991 50.44 4094.428 0
211101008 1 1 1991 52.1 4094.428 0
211101008 3 1 1991 26.39 4094.428 0
211101008 3 1 1991 26.39 4094.428 0
211101008 1 1 1991 52.1 4094.428 0
211101008 3 1 1991 26.39 4094.428 0
211101008 21 2 1991 21.64 4094.428 0
211101008 3 1 1991 26.39 4094.428 0
211101008 2 2 1991 50.44 4094.428 0
211101008 2 2 1991 50.44 4094.428 0
211101008 2 2 1991 50.44 4094.428 0
211101008 4 2 1991 28.72 4094.428 0
211101008 2 2 1991 50.44 4094.428 0
211101008 3 1 1991 26.39 4094.428 0
211101008 5 1991 4094.428 0
211101008 3 1 1991 26.39 4094.428 0
211101008 4 2 1991 28.72 4094.428 0
211101008 1 1 1991 52.1 4094.428 0
211101008 2 2 1991 50.44 4094.428 0
211101008 2 2 1991 50.44 4094.428 0
211101008 3 1 1991 26.39 4094.428 0
211101008 2 2 1991 50.44 4094.428 0
211101008 3 1 1991 26.39 4094.428 0
211101008 1 1 1991 52.1 4094.428 0
211101008 1 1 1991 52.1 4094.428 0
211101008 1 1 1991 52.1 4094.428 0
211101008 1 1 1991 52.1 4094.428 0
211101008 3 1 1991 26.39 4094.428 0
211101008 1 1 1991 52.1 4094.428 0
211101008 1 1 1991 52.1 4094.428 0
211101008 1 1 1991 52.1 4094.428 0
211101008 1 1 1991 52.1 4094.428 0
211101008 1 1 1993 54.02 5254.783 0
211101008 1 1 1993 54.02 5254.783 0
211101008 5 1 1993 5.97 5254.783 0
211101008 2 2 1993 52.36 5254.783 0
211101008 2 2 1993 52.36 5254.783 0
211101008 1 1 1993 54.02 5254.783 0
211101008 1 1 1993 54.02 5254.783 0
211101008 2 2 1993 52.36 5254.783 0
211101008 1 1 1993 54.02 5254.783 0
211101008 2 2 1993 52.36 5254.783 0
211101008 1 1 1993 54.02 5254.783 0
211101008 3 1 1993 28.31 5254.783 0
211101008 3 1 1993 28.31 5254.783 0
211101008 1 1 1993 54.02 5254.783 0
211101008 2 2 1993 52.36 5254.783 0
211101008 2 2 1993 52.36 5254.783 0
211101008 2 2 1993 52.36 5254.783 0
211101008 4 2 1993 30.64 5254.783 0
211101008 2 2 1993 52.36 5254.783 0
211101008 2 2 1993 52.36 5254.783 0
211101008 21 2 1993 23.56 5254.783 0
211101008 1 1 1993 54.02 5254.783 0
211101008 1 1 1993 54.02 5254.783 0
211101008 1 1 1993 54.02 5254.783 0
211101008 3 1 1993 28.31 5254.783 0
211101008 3 1 1993 28.31 5254.783 0
211101008 2 2 1993 52.36 5254.783 0
211101008 3 1 1993 28.31 5254.783 0
211101008 3 1 1993 28.31 5254.783 0
211101008 1 1 1993 54.02 5254.783 0
211101008 1 1 1993 54.02 5254.783 0
211101008 2 2 1993 52.36 5254.783 0
211101008 2 2 1993 52.36 5254.783 0
211101008 2 2 1993 52.36 5254.783 0
211101008 4 2 1993 30.64 5254.783 0
211101008 1 1 1993 54.02 5254.783 0
211101008 3 1 1993 28.31 5254.783 0
211101008 3 1 1993 28.31 5254.783 0
211101008 3 1 1993 28.31 5254.783 0
211101008 3 1 1993 28.31 5254.783 0
211101008 3 1 1993 28.31 5254.783 0
211101008 2 2 1993 52.36 5254.783 0
211101008 21 2 1993 23.56 5254.783 0
211101008 1 1 1993 54.02 5254.783 0
211101008 2 2 1993 52.36 5254.783 0
211101008 1 1 1993 54.02 5254.783 0
211101008 1 1 2000 60.91 2309.891 1
211101008 3 1 2000 35.2 2309.891 1
211101008 2 2 2000 59.25 2309.891 1
211101008 21 2 2000 30.44 2309.891 1
211101008 5 1 2000 12.86 2309.891 1
211101008 3 1 2000 35.2 2309.891 1
211101008 3 1 2000 35.2 2309.891 1
211101008 61 1 2000 6.19 2309.891 1
211101008 2 2 2000 59.25 2309.891 1
211101008 21 2 2000 30.44 2309.891 1
211101008 1 1 2000 60.91 2309.891 1
211101008 3 1 2000 35.2 2309.891 1
211101008 3 1 2000 35.2 2309.891 1
211101008 2 2 2000 59.25 2309.891 1
211101008 1 1 2000 60.91 2309.891 1
211101008 4 2000 2309.891 1
211101008 3 1 2000 35.2 2309.891 1
211101008 1 1 2000 60.91 2309.891 1
211101008 2 2 2000 59.25 2309.891 1
211101008 1 1 2000 60.91 2309.891 1
211101008 2 2 2000 59.25 2309.891 1
211101008 2 2 2000 59.25 2309.891 1
211101008 2 2 2000 59.25 2309.891 1
211101008 1 1 2000 60.91 2309.891 1
211101008 1 1 2000 60.91 2309.891 1
211101008 3 1 2004 39.27 5357.143 0
211101008 2 2 2004 63.33 5357.143 0
211101008 21 2 2004 34.52 5357.143 0
211101008 2 2 2004 63.33 5357.143 0
211101008 3 1 2004 39.27 5357.143 0
211101008 2 2 2004 63.33 5357.143 0
211101008 2 2 2004 63.33 5357.143 0
211101008 3 1 2004 39.27 5357.143 0
211101008 81 2004 5357.143 0
211101008 61 1 2004 10.26 5357.143 0
211101008 1 1 2004 64.98 5357.143 0
211101008 2 2 2004 63.33 5357.143 0
211101008 21 2 2004 34.52 5357.143 0
211101008 1 1 2004 64.98 5357.143 0
211101008 1 1 2004 64.98 5357.143 0
211101008 5 2004 5357.143 0
211101008 1 1 2004 64.98 5357.143 0
211101008 2 2 2006 65.28 5084.217 0
211101008 1 1 2006 66.94 5084.217 0
211101008 81 2006 5084.217 0
211101008 61 1 2006 12.22 5084.217 0
211101008 2 2 2006 65.28 5084.217 0
211101008 3 1 2006 41.23 5084.217 0
211101008 1 1 2006 66.94 5084.217 0
211101008 3 1 2006 41.23 5084.217 0
211101008 101 2006 5084.217 0
211101008 5 1 2006 18.89 5084.217 0
211101008 2 2 2006 65.28 5084.217 0
211101008 102 2006 5084.217 0
211101008 1 1 2006 66.94 5084.217 0
211101008 21 2 2006 36.47 5084.217 0
211101008 4 2006 5084.217 0
211101008 1 1 2006 66.94 5084.217 0
211101008 2 2 2006 65.28 5084.217 0
211101008 1 1 2006 66.94 5084.217 0
211101008 21 2 2006 36.47 5084.217 0
211101008 2 2 2006 65.28 5084.217 0
211101008 21 2 2009 39.61 13756 0
211101008 2 2 2009 68.41 13756 0
211101008 102 2009 13756 0
211101008 101 2009 13756 0
211101008 2 2 2009 68.41 13756 0
211101008 3 2009 13756 0
211101008 1 2009 13756 0
211101008 2 2 2009 68.41 13756 0
211101008 4 2009 13756 0
211101008 2 2 2009 68.41 13756 0
211101008 1 2009 13756 0
211101008 1 2009 13756 0
211101008 5 1 2009 22.02 13756 0
211101008 2 2 2009 68.41 13756 0
211101008 61 1 2009 15.35 13756 0
end
*--------------------------------------------------------------


  duplicates drop hhid line wave, force // 删除重复样本
  sort hhid line wave

  bysort hhid line: gen  pk0 = wave if pk[_n-1]==0&pk[_n]==1
  bysort hhid line: egen pk0yr = min(pk0)  // 进入贫困时间
  bysort hhid line: gen pkT = wave if pk[_n-1]==1&pk[_n]==0
  bysort hhid line: egen pkTyr = max(pkT)  // 脱贫时间
  gen pk_duration = pkTyr-pk0yr          // 持续时间
  replace pk_duration=0 if pk_duration==.  // 贫困状态未发生变化的观察值

  order hhid-pk  pk0 pkT pk0yr pkTyr pk_dura

  browse

扫码加我拉你入群

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

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

arlionn

2011-11-1 11:20:25

你可以逐条执行命令，查看数据的变化，以便确认我写的命令是否有误。

扫码加我拉你入群

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

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

susanketty

2011-11-1 12:19:30

感谢连老师！！试过了，没什么问题！
还想请教您：
1、如果像这样调查时间有间隔的，不是每年连续调查的，会不会不能使用生存分析，或是对结果有一定的影响？因为贫困时间变量，计算上是没有问题，但是从准确意义上是不对的，因为时间有间断，间断年份的情况是不清楚的。
2、您命令中对于家庭成员分析，是不是因为数据里有关于家庭成员的变量，如果这些变量不是必须，我直接对家庭整体分析是否也是可以的，只不过将您命令中有关于家庭成员的换为家庭的？
3、如果下一步运用生存分析对于家庭进行分析，命令是否这样书写：
stset pk_duration,failure(pk) id (hhid)
感谢。。。。。。

扫码加我拉你入群

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

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

arlionn

2011-11-1 14:26:26

我对生存分析知之甚少，不敢妄言。
你需要搜索一些采用了类似数据研究类似问题的论文，看看他们是怎么处理的。

扫码加我拉你入群

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

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

点击查看更多内容…

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

分享

扫码加好友，拉您进群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群