[求助]在使用logistic回归中遇到的一个问题

4598

收藏 2009-06-07

我在使用logistic回归中遇到一个问题：（以下为我的应用过程及问题）gen fardeau=. (686 missing values generated) . replace fardeau=1 if totzarit&lt;=20 (320 real changes made) . replace fardeau=2 if totzarit&gt;=21 &amp; totzarit&lt;=40 (235 real changes made) . replace fardeau=3 if totzarit&gt;=41 &amp; totzarit&lt;=60 (65 real changes made) . replace fardeau=4 if totzarit&gt;=61 &amp; totzarit&lt;=88 (18 real changes made) . tab fardeau,missing &nbsp;&nbsp; fardeau |&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;Freq.&nbsp;&nbsp;&nbsp;&nbsp; Percent&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;Cum. ------------+----------------------------------- &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 1 |&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;320&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 46.65&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 46.65 &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 2 |&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;235&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 34.26&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 80.90 &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 3 |&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 65&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;9.48&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 90.38 &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 4 |&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 18&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;2.62&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 93.00 &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; . |&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 48&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;7.00&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;100.00 ------------+----------------------------------- &nbsp;&nbsp;&nbsp;&nbsp; Total |&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;686&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;100.00 xi:logit stable i.fardeau i.fardeau&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; _Ifardeau_1-4&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; (naturally coded; _Ifardeau_1 omitted) note: _Ifardeau_4 != 0 predicts failure perfectly &nbsp;&nbsp;&nbsp;&nbsp; _Ifardeau_4 dropped and 1 obs not used（我想找到一个命令，可以阻止_Ifardeau_4被删除，而且象前3项一样得出数据分析。不知哪位高手可以指点我这个菜鸟） Iteration 0:&nbsp;&nbsp; log likelihood = -95.043737 Iteration 1:&nbsp;&nbsp; log likelihood =&nbsp;&nbsp;-94.13092 Iteration 2:&nbsp;&nbsp; log likelihood = -94.070877 Iteration 3:&nbsp;&nbsp; log likelihood = -94.070021 Iteration 4:&nbsp;&nbsp; log likelihood = -94.070021 Logistic regression&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; Number of obs&nbsp;&nbsp; =&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;170 &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; LR chi2(2)&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;=&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 1.95 &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; Prob &gt; chi2&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp;&nbsp;&nbsp;&nbsp; 0.3777 Log likelihood = -94.070021&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; Pseudo R2&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp;&nbsp;&nbsp;&nbsp; 0.0102 ------------------------------------------------------------------------------ &nbsp;&nbsp;&nbsp;&nbsp; stable |&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;Coef.&nbsp;&nbsp; Std. Err.&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;z&nbsp;&nbsp;&nbsp;&nbsp;P&gt;|z|&nbsp;&nbsp;&nbsp;&nbsp; [95% Conf. Interval] -------------+---------------------------------------------------------------- _Ifardeau_2 |&nbsp;&nbsp;-.1040009&nbsp;&nbsp; .3706763&nbsp;&nbsp;&nbsp;&nbsp;-0.28&nbsp;&nbsp; 0.779&nbsp;&nbsp;&nbsp;&nbsp;-.8305131&nbsp;&nbsp;&nbsp;&nbsp;.6225113 _Ifardeau_3 |&nbsp;&nbsp;-1.287354&nbsp;&nbsp; 1.074473&nbsp;&nbsp;&nbsp;&nbsp;-1.20&nbsp;&nbsp; 0.231&nbsp;&nbsp;&nbsp;&nbsp;-3.393283&nbsp;&nbsp;&nbsp;&nbsp;.8185743 &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;_cons |&nbsp;&nbsp;-1.015231&nbsp;&nbsp; .2334368&nbsp;&nbsp;&nbsp;&nbsp;-4.35&nbsp;&nbsp; 0.000&nbsp;&nbsp;&nbsp;&nbsp;-1.472758&nbsp;&nbsp; -.5577029 ------------------------------------------------------------------------------

[此贴子已经被作者于2009-6-7 21:14:49编辑过]

扫码加我拉你入群

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

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

全部回复

eblog

2009-6-7 21:28:00

把虚拟变量的相关知识看看

扫码加我拉你入群

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

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

nolen56

2009-6-7 22:09:00

回复：（eblog）把虚拟变量的相关知识看看

以下是引用eblog在2009-6-7 21:28:00的发言：

把虚拟变量的相关知识看看

请问楼上的高手，可否说得更具体些，我应该看虚拟变量相关的那些内容？本人确实是能力有限，特在此求教。

扫码加我拉你入群

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

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

sungmoo

2009-6-8 07:54:00

可不可以先把你的stable与totzarit的数据帖出来？

扫码加我拉你入群

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

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

sungmoo

2009-6-8 08:57:00

以下是引用nolen56在2009-6-7 21:13:00的发言：

gen fardeau=.
(686 missing values generated)

. replace fardeau=1 if totzarit<=20
(320 real changes made)

. replace fardeau=2 if totzarit>=21 & totzarit<=40
(235 real changes made)

. replace fardeau=3 if totzarit>=41 & totzarit<=60
(65 real changes made)

. replace fardeau=4 if totzarit>=61 & totzarit<=88
(18 real changes made)

*这几步可以一步实现：

recode totzarit (min/20=1)(21/40=2)(41/60=3)(61/88=4)(else=.),g(fardeau)

[此贴子已经被作者于2009-6-8 8:59:02编辑过]

扫码加我拉你入群

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

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

nolen56

2009-6-8 18:00:00

回复：（sungmoo）可不可以先把你的stable与totzarit的...

Gen vmms=totmms4a - totmms

sum vmms
Variable Obs Mean Std. Dev. Min Max
vmms 184 -7.141304 5.745193 -26 9

list numcp totmms4a vmms if vmms==. & totmms4a!=.
Numcp Totmms Totmms4a Vmms
16141 . 16 .

3)
sum vmms,det
gen vmms4=.
replace vmms4=1 if vmms<=-11
replace vmms4=2 if vmms>-11 & vmms<=-7
replace vmms4=3 if vmms>-7 & vmms<=-3
replace vmms4=4 if vmms>-3 & vmms<.

tab vmms4
vmms4 Freq. Percent Cum.
1 51 27.72 27.72
2 48 26.09 53.80
3 40 21.74 75.54
4 45 24.46 100.00
Totle 184 100.00

4)
gen stable=0
replace stable=1 if vmms4==4
replace stable=. if vmms==.

tab stable, miss
stable Freq. Percent Cum.
0 139 20.26 20.26
1 45 6.56 26.82
. 502 73.18 100.00
Total 686 100.00

sum totzarit

Variable | Obs Mean Std. Dev. Min Max
-------------+--------------------------------------------------------
totzarit | 638 22.60815 15.8566 0 82

. sum totzarit,det

zarit total t0
-------------------------------------------------------------
Percentiles Smallest
1% 0 0
5% 2 0
10% 4 0 Obs 638
25% 10 0 Sum of Wgt. 638

50% 20 Mean 22.60815
Largest Std. Dev. 15.8566
75% 33 71
90% 44 71 Variance 251.4318
95% 55 72 Skewness .8159247
99% 68 82 Kurtosis 3.345145

还请sungmoo多指教

扫码加我拉你入群

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

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

点击查看更多内容…

扫码加我拉你入群

扫码加我拉你入群

回复：（eblog）把虚拟变量的相关知识看看

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

回复：（sungmoo）原始数据？

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加好友，拉您进群


	1	2	3	4	Total

0	72	46	10	2	130
1	26	15	1	0	42

Total	98	61	11	2	172

扫码加我 拉你入群

扫码加我 拉你入群

回复：（eblog）把虚拟变量的相关知识看看

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

回复：（sungmoo）原始数据？

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

分享

扫码加好友，拉您进群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群