Page 283. The coefficients at the top of the page.
use https://stats.idre.ucla.edu/stat/stata/examples/cama4/depress, clear
logit cases age income
Iteration 0: log likelihood = -134.06225
Iteration 1: log likelihood = -127.42024
Iteration 2: log likelihood = -127.01794
Iteration 3: log likelihood = -127.01305
Iteration 4: log likelihood = -127.01304
Logit estimates Number of obs = 294
LR chi2(2) = 14.10
Prob > chi2 = 0.0009
Log likelihood = -127.01304 Pseudo R2 = 0.0526
------------------------------------------------------------------------------
cases | Coef. Std. Err. z P>|z| [95% Conf. Interval]
-------------+----------------------------------------------------------------
age | -.0201671 .0088966 -2.27 0.023 -.0376041 -.0027301
income | -.0413479 .0140587 -2.94 0.003 -.0689025 -.0137933
_cons | .0279774 .4872007 0.06 0.954 -.9269184 .9828732
------------------------------------------------------------------------------
Table 12.1, page 285.
tab sex cases
| depressed is cesd >=16
|
sex | normal depressed | Total
-----------+----------------------+----------
male | 101 10 | 111
female | 143 40 | 183
-----------+----------------------+----------
Total | 244 50 | 294
Page 286. The odds ratios at the top of the page.
logit cases sex, or
Iteration 0: log likelihood = -134.06225
Iteration 1: log likelihood = -129.83832
Iteration 2: log likelihood = -129.69929
Iteration 3: log likelihood = -129.69883
Logit estimates Number of obs = 294
LR chi2(1) = 8.73
Prob > chi2 = 0.0031
Log likelihood = -129.69883 Pseudo R2 = 0.0325
------------------------------------------------------------------------------
cases | Odds Ratio Std. Err. z P>|z| [95% Conf. Interval]
-------------+----------------------------------------------------------------
sex | 2.825175 1.064196 2.76 0.006 1.350243 5.911243
------------------------------------------------------------------------------
Page 286. The estimates at the bottom of the page.
gen sex1 = sex - 1
logit cases sex1
Iteration 0: log likelihood = -134.06225
Iteration 1: log likelihood = -129.83832
Iteration 2: log likelihood = -129.69929
Iteration 3: log likelihood = -129.69883
Logit estimates Number of obs = 294
LR chi2(1) = 8.73
Prob > chi2 = 0.0031
Log likelihood = -129.69883 Pseudo R2 = 0.0325
------------------------------------------------------------------------------
cases | Coef. Std. Err. z P>|z| [95% Conf. Interval]
-------------+----------------------------------------------------------------
sex1 | 1.03857 .3766834 2.76 0.006 .3002844 1.776856
_cons | -2.312535 .3315077 -6.98 0.000 -2.962279 -1.662792
------------------------------------------------------------------------------
Page 287. Top of the page.
logit cases age income sex1
Iteration 0: log likelihood = -134.06225
Iteration 1: log likelihood = -124.45941
Iteration 2: log likelihood = -123.78013
Iteration 3: log likelihood = -123.76979
Iteration 4: log likelihood = -123.76979
Logit estimates Number of obs = 294
LR chi2(3) = 20.58
Prob > chi2 = 0.0001
Log likelihood = -123.76979 Pseudo R2 = 0.0768
------------------------------------------------------------------------------
cases | Coef. Std. Err. z P>|z| [95% Conf. Interval]
-------------+----------------------------------------------------------------
age | -.0209575 .0090401 -2.32 0.020 -.0386758 -.0032392
income | -.0365635 .0140897 -2.60 0.009 -.0641787 -.0089482
sex1 | .9294487 .3858256 2.41 0.016 .1732444 1.685653
_cons | -.6764559 .5788124 -1.17 0.243 -1.810907 .4579955
------------------------------------------------------------------------------
Page 287. Table of coefficients and standard errors at the bottom.
logit cases age income
Iteration 0: log likelihood = -134.06225
Iteration 1: log likelihood = -127.42024
Iteration 2: log likelihood = -127.01794
Iteration 3: log likelihood = -127.01305
Iteration 4: log likelihood = -127.01304
Logit estimates Number of obs = 294
LR chi2(2) = 14.10
Prob > chi2 = 0.0009
Log likelihood = -127.01304 Pseudo R2 = 0.0526
------------------------------------------------------------------------------
cases | Coef. Std. Err. z P>|z| [95% Conf. Interval]
-------------+----------------------------------------------------------------
age | -.0201671 .0088966 -2.27 0.023 -.0376041 -.0027301
income | -.0413479 .0140587 -2.94 0.003 -.0689025 -.0137933
_cons | .0279774 .4872007 0.06 0.954 -.9269184 .9828732
------------------------------------------------------------------------------
Page 288. These numbers are obtained from the output from page 287.
Page 290. Table at the top of the page.
gen duminc = 1
replace duminc = 0 if income >=10
gen dumemp = 0
replace dumemp = 1 if employ == 2 | employ == 3
replace dumemp = . if employ == 7
logit cases duminc dumemp
Iteration 0: log likelihood = -131.73021
Iteration 1: log likelihood = -127.6761
Iteration 2: log likelihood = -127.42821
Iteration 3: log likelihood = -127.42796
Logit estimates Number of obs = 290
LR chi2(2) = 8.60
Prob > chi2 = 0.0135
Log likelihood = -127.42796 Pseudo R2 = 0.0327
------------------------------------------------------------------------------
cases | Coef. Std. Err. z P>|z| [95% Conf. Interval]
-------------+----------------------------------------------------------------
duminc | .2722943 .3376895 0.81 0.420 -.3895651 .9341536
dumemp | 1.028483 .3487117 2.95 0.003 .3450204 1.711945
_cons | -1.934537 .2259083 -8.56 0.000 -2.377309 -1.491765
------------------------------------------------------------------------------
Page 291. Table in the middle of the page.
gen dincemp = duminc*dumemp
logit cases duminc dumemp dincemp
Iteration 0: log likelihood = -131.73021
Iteration 1: log likelihood = -125.18993
Iteration 2: log likelihood = -123.52872
Iteration 3: log likelihood = -123.31292
Iteration 4: log likelihood = -123.31287
Logit estimates Number of obs = 290
LR chi2(3) = 16.83
Prob > chi2 = 0.0008
Log likelihood = -123.31287 Pseudo R2 = 0.0639
------------------------------------------------------------------------------
cases | Coef. Std. Err. z P>|z| [95% Conf. Interval]
-------------+----------------------------------------------------------------
duminc | -.3756121 .4349381 -0.86 0.388 -1.228075 .4768508
dumemp | .317535 .4520206 0.70 0.482 -.5684091 1.203479
dincemp | 2.198143 .7888294 2.79 0.005 .6520661 3.74422
_cons | -1.734601 .2214037 -7.83 0.000 -2.168544 -1.300658
------------------------------------------------------------------------------
Page 292. Bottom of the page.
di 2.198143/.7888294
2.7865886
di (2.7865886)^2
7.765076
logit cases duminc dumemp
Iteration 0: log likelihood = -131.73021
Iteration 1: log likelihood = -127.6761
Iteration 2: log likelihood = -127.42821
Iteration 3: log likelihood = -127.42796
Logit estimates Number of obs = 290
LR chi2(2) = 8.60
Prob > chi2 = 0.0135
Log likelihood = -127.42796 Pseudo R2 = 0.0327
------------------------------------------------------------------------------
cases | Coef. Std. Err. z P>|z| [95% Conf. Interval]
-------------+----------------------------------------------------------------
duminc | .2722943 .3376895 0.81 0.420 -.3895651 .9341536
dumemp | 1.028483 .3487117 2.95 0.003 .3450204 1.711945
_cons | -1.934537 .2259083 -8.56 0.000 -2.377309 -1.491765
------------------------------------------------------------------------------
logit cases duminc dumemp dincemp
Iteration 0: log likelihood = -131.73021
Iteration 1: log likelihood = -125.18993
Iteration 2: log likelihood = -123.52872
Iteration 3: log likelihood = -123.31292
Iteration 4: log likelihood = -123.31287
Logit estimates Number of obs = 290
LR chi2(3) = 16.83
Prob > chi2 = 0.0008
Log likelihood = -123.31287 Pseudo R2 = 0.0639
------------------------------------------------------------------------------
cases | Coef. Std. Err. z P>|z| [95% Conf. Interval]
-------------+----------------------------------------------------------------
duminc | -.3756121 .4349381 -0.86 0.388 -1.228075 .4768508
dumemp | .317535 .4520206 0.70 0.482 -.5684091 1.203479
dincemp | 2.198143 .7888294 2.79 0.005 .6520661 3.74422
_cons | -1.734601 .2214037 -7.83 0.000 -2.168544 -1.300658
------------------------------------------------------------------------------
di 16.83-8.6
8.23
Page 298. Middle of the page.
gen age0 = 0
replace age0 = 1 if age < 28
gen age1 = 0
replace age1 = 1 if age >=28 & age <= 42
gen age2 = 0
replace age2 = 1 if age >=43 & age <= 58
gen age3 = 0
replace age3 = 1 if age >=59 & age <= 89
logit cases age1 age2 age3 income sex
Iteration 0: log likelihood = -134.06225
Iteration 1: log likelihood = -124.33006
Iteration 2: log likelihood = -123.58989
Iteration 3: log likelihood = -123.57661
Iteration 4: log likelihood = -123.57661
Logit estimates Number of obs = 294
LR chi2(5) = 20.97
Prob > chi2 = 0.0008
Log likelihood = -123.57661 Pseudo R2 = 0.0782
------------------------------------------------------------------------------
cases | Coef. Std. Err. z P>|z| [95% Conf. Interval]
-------------+----------------------------------------------------------------
age1 | .0747259 .4318019 0.17 0.863 -.7715902 .9210421
age2 | -.5706276 .4743977 -1.20 0.229 -1.50043 .3591747
age3 | -.8853103 .4563059 -1.94 0.052 -1.779654 .0090329
income | -.0379521 .014852 -2.56 0.011 -.0670616 -.0088426
sex | .9238035 .3864416 2.39 0.017 .1663919 1.681215
_cons | -2.159495 .7830424 -2.76 0.006 -3.69423 -.6247601
------------------------------------------------------------------------------
Figure 12.2, page 299.
gen female = (sex == 2)
logit cases age1 age2 age3 income female
matrix b = e(b)
gen b0 = 0 in 1
svmat b
list b0 - b6 in 1
* The point of the next two commands is to drop all of the
* unnecessary variables and observations. You may need to
* modify the variable names on the drop command to match
* your variables.
drop sex - age3
drop in 2/l
reshape long b, i(id)
list
gen newage = 22.5
replace newage = 35 if _j == 1
replace newage = 50.5 if _j == 2
replace newage = 74 if _j == 3
* You need to calculate the midpoint of each interval based on the
* table on page 298. The minimum of age is 18.
graph twoway scatter b newage in 1/4, msymbol(O) connect(L) ylabel(-1(.5)0, nogrid) ///
ytitle(Coefficient b) xlabel(20(10)80) xtitle(Age)
Figure 12.3, page 302.
use https://stats.idre.ucla.edu/stat/stata/examples/cama4/depress, clear logit cases sex income age predict p predict db, db graph twoway scatter db p
Table 12.2, page 303.
logit cases age income sex
Iteration 0: log likelihood = -134.06225
Iteration 1: log likelihood = -124.45941
Iteration 2: log likelihood = -123.78013
Iteration 3: log likelihood = -123.76979
Iteration 4: log likelihood = -123.76979
Logit estimates Number of obs = 294
LR chi2(3) = 20.58
Prob > chi2 = 0.0001
Log likelihood = -123.76979 Pseudo R2 = 0.0768
------------------------------------------------------------------------------
cases | Coef. Std. Err. z P>|z| [95% Conf. Interval]
-------------+----------------------------------------------------------------
age | -.0209575 .0090401 -2.32 0.020 -.0386758 -.0032392
income | -.0365635 .0140897 -2.60 0.009 -.0641787 -.0089482
sex | .9294487 .3858256 2.41 0.016 .1732444 1.685653
_cons | -1.605905 .8465372 -1.90 0.058 -3.265087 .0532779
------------------------------------------------------------------------------
list db in -4/l
+----------+
| db |
|----------|
291. | .0004829 |
292. | .0001442 |
293. | .0006498 |
294. | .0111896 |
+----------+
gen x = 0
replace x = 3 if db > .1637310
(5 real changes made)
replace x = 2 if db > .1789604
(4 real changes made)
replace x = 1 if db > .2084397
(3 real changes made)
logit cases age income sex if x ~= 1
Iteration 0: log likelihood = -131.91495
Iteration 1: log likelihood = -121.29438
Iteration 2: log likelihood = -120.39934
Iteration 3: log likelihood = -120.3796
Iteration 4: log likelihood = -120.37958
Logit estimates Number of obs = 291
LR chi2(3) = 23.07
Prob > chi2 = 0.0000
Log likelihood = -120.37958 Pseudo R2 = 0.0874
------------------------------------------------------------------------------
cases | Coef. Std. Err. z P>|z| [95% Conf. Interval]
-------------+----------------------------------------------------------------
age | -.0215037 .0091147 -2.36 0.018 -.0393683 -.0036391
income | -.0416866 .0149928 -2.78 0.005 -.071072 -.0123013
sex | 1.022995 .4007114 2.55 0.011 .2376152 1.808375
_cons | -1.692541 .8721181 -1.94 0.052 -3.401861 .016779
------------------------------------------------------------------------------
logit cases age income sex if x ~= 2
Iteration 0: log likelihood = -132.28233
Iteration 1: log likelihood = -121.82603
Iteration 2: log likelihood = -121.03222
Iteration 3: log likelihood = -121.01885
Iteration 4: log likelihood = -121.01885
Logit estimates Number of obs = 293
LR chi2(3) = 22.53
Prob > chi2 = 0.0001
Log likelihood = -121.01885 Pseudo R2 = 0.0851
------------------------------------------------------------------------------
cases | Coef. Std. Err. z P>|z| [95% Conf. Interval]
-------------+----------------------------------------------------------------
age | -.023414 .0092536 -2.53 0.011 -.0415507 -.0052774
income | -.0357778 .0141515 -2.53 0.011 -.0635142 -.0080415
sex | 1.050539 .4007857 2.62 0.009 .2650138 1.836065
_cons | -1.756961 .871175 -2.02 0.044 -3.464433 -.0494895
------------------------------------------------------------------------------
logit cases age income sex if x ~= 3
Iteration 0: log likelihood = -132.28233
Iteration 1: log likelihood = -121.54471
Iteration 2: log likelihood = -120.67498
Iteration 3: log likelihood = -120.65779
Iteration 4: log likelihood = -120.65778
Logit estimates Number of obs = 293
LR chi2(3) = 23.25
Prob > chi2 = 0.0000
Log likelihood = -120.65778 Pseudo R2 = 0.0879
------------------------------------------------------------------------------
cases | Coef. Std. Err. z P>|z| [95% Conf. Interval]
-------------+----------------------------------------------------------------
age | -.0228888 .0092002 -2.49 0.013 -.0409209 -.0048567
income | -.0389067 .0145405 -2.68 0.007 -.0674056 -.0104078
sex | 1.04192 .4008515 2.60 0.009 .2562651 1.827574
_cons | -1.71381 .8720871 -1.97 0.049 -3.423069 -.0045506
------------------------------------------------------------------------------
logit cases age income sex if x == 0
Iteration 0: log likelihood = -128.31745
Iteration 1: log likelihood = -115.44875
Iteration 2: log likelihood = -114.1166
Iteration 3: log likelihood = -114.074
Iteration 4: log likelihood = -114.07392
Logit estimates Number of obs = 289
LR chi2(3) = 28.49
Prob > chi2 = 0.0000
Log likelihood = -114.07392 Pseudo R2 = 0.1110
------------------------------------------------------------------------------
cases | Coef. Std. Err. z P>|z| [95% Conf. Interval]
-------------+----------------------------------------------------------------
age | -.0262663 .0095311 -2.76 0.006 -.0449469 -.0075856
income | -.0439743 .015642 -2.81 0.005 -.074632 -.0133166
sex | 1.30267 .4406845 2.96 0.003 .4389443 2.166396
_cons | -2.018118 .9403468 -2.15 0.032 -3.861163 -.1750718
------------------------------------------------------------------------------
Table 12.3, page 304.
list age income sex cases p if id == 288 | id==99 | id==143 | id == 232 | id == 68
+---------------------------------------+
| age income sex cases p |
|---------------------------------------|
68. | 40 45 1 1 .0406946 |
99. | 72 11 1 1 .069941 |
143. | 40 45 1 0 .0406946 |
232. | 40 45 1 0 .0406946 |
288. | 61 28 1 1 .0484001 |
+---------------------------------------+
Figure 12.4, page 305.
logit cases sex income age
Iteration 0: log likelihood = -134.06225
Iteration 1: log likelihood = -124.45941
Iteration 2: log likelihood = -123.78013
Iteration 3: log likelihood = -123.76979
Iteration 4: log likelihood = -123.76979
Logit estimates Number of obs = 294
LR chi2(3) = 20.58
Prob > chi2 = 0.0001
Log likelihood = -123.76979 Pseudo R2 = 0.0768
------------------------------------------------------------------------------
cases | Coef. Std. Err. z P>|z| [95% Conf. Interval]
-------------+----------------------------------------------------------------
sex | .9294487 .3858256 2.41 0.016 .1732444 1.685653
income | -.0365635 .0140897 -2.60 0.009 -.0641787 -.0089482
age | -.0209575 .0090401 -2.32 0.020 -.0386758 -.0032392
_cons | -1.605905 .8465372 -1.90 0.058 -3.265087 .0532779
------------------------------------------------------------------------------
predict dx2, dx2
graph twoway scatter dx2 p [w = db], msymbol(Oh) ylabel(0(5)25, nogrid) xlabel(0(.1).5)
Figure 12.5, page 307.
lsens
Figure 12.6, page 307.
lroc
