. use d:\统计程序安装\stata数据\nhefs
. svyset psu2[pw=swgt2],strata(strata2)
pweight: swgt2
VCE: linearized
Single unit: missing
Strata 1: strata2
SU 1: psu2
FPC 1: <zero>
. stset age_lung_cancer[pw=swgt2],fail(lung_cancer)
failure event: lung_cancer != 0 & lung_cancer < .
obs. time interval: (0, age_lung_cancer]
exit on or before: failure
weight: [pweight=swgt2]
14407 total obs.
5126 event time missing (age_lung_cancer>=.) PROBABLE ERROR
------------------------------------------------------------------------------
9281 obs. remaining, representing
83 failures in single record/single failure data
599691 total analysis time at risk, at risk from t = 0
earliest observed entry t = 0
last observed exit t = 97
. svy:stcox former_smoker smoker male urban1 rural
(running stcox on estimation sample)
Survey: Cox regression
Number of strata = 35 Number of obs = 9149
Number of PSUs = 105 Population size = 1.513e+08
Design df = 70
F( 5, 66) = 14.07
Prob > F = 0.0000
------------------------------------------------------------------------------
| Linearized
_t | Haz. Ratio Std. Err. t P>|t| [95% Conf. Interval]
-------------+----------------------------------------------------------------
former_smo~r | 2.788113 .6205102 4.61 0.000 1.788705 4.345923
smoker | 7.849483 2.593249 6.24 0.000 4.061457 15.17051
male | 1.187611 .3445315 0.59 0.555 .6658757 2.118142
urban1 | .8035074 .3285144 -0.54 0.594 .3555123 1.816039
rural | 1.581674 .5281859 1.37 0.174 .8125799 3.078702
这是stata中自带例子的程序及软件输出的结果,我觉得和我所模拟的一样呀,但是不知道怎么就是最后拟合cox模型时出不来结果,请大家帮帮忙,不吝赐教,谢谢。