Page 283 The coefficients at the top of the page
get file 'c:cama4depress.sav'.
logistic regression cases /method = enter age income.
<some output omitted>
Page 283 Figure 12.1 Logistic function for the depression data set
NOTE: We were unable to reproduce this graph.
Page 285 Table 12.1 Classification of individuals by depression level and sex
crosstabs /tables = sex by cases.
Page 286 Odds ratios and coefficients
compute sex1 = sex - 1. logistic regression cases with sex1.
<some output omitted>
Page 287 Table of coefficients and standard errors
logistic regression cases /method = enter age income.
<some output omitted>
Page 288 These numbers are obtained from the output from page 287.
Page 290 Table at the top of the page
compute duminc = 1. if income >=10 duminc = 0. compute dumemp = 0. if employ = 2 dumemp = 1. if employ = 3 dumemp = 1. if employ = 7 dumemp = $sysmis. exe.
logistic cases with duminc dumemp.
<some output omitted>
Page 291 Table in the middle of the page
compute dincemp = duminc*dumemp. logistic cases with duminc dumemp dincemp.
<some output omitted>
Page 292 bottom of the page
NOTE: The likelihood ratio chi-square values needed are given in the output for the two models shown above: 16.83-8.6 = 8.23.
Page 298 middle of the page
compute age0 = 0. if age < 28 age0 = 1. compute age1 = 0. if age >=28 and age <= 42 age1 = 1. compute age2 = 0. if age >=43 and age <= 58 age2 = 1. compute age3 = 0. if age >=59 and age <= 89 age3 = 1. exe.
logistic cases with age1 age2 age3 income sex.
<some output omitted>
Page 299 Figure 12.2 Estimated coefficients for age quartiles by midpoint of the quartile
data list list /coef newage. begin data. 0 22.5 .0747259 35 -.5706276 50.5 -.8853103 74 end data.
formats coef (f3.1) newage (f2.0). GGRAPH /GRAPHDATASET NAME="GraphDataset" VARIABLES= coef newage /GRAPHSPEC SOURCE=INLINE . BEGIN GPL SOURCE: s=userSource( id( "GraphDataset" ) ) DATA: coef=col( source(s), name( "coef" ) ) DATA: newage=col( source(s), name( "newage" ) ) GUIDE: text.title( label( "Figure 12.2" ) ) GUIDE: axis( dim( 1 ), label( "New Age" ) ) GUIDE: axis( dim( 2 ), label( "Coeff" ), delta(.5) ) SCALE: linear(dim(2), min(-1), max(0)) ELEMENT: point( position( newage * coef ) ) ELEMENT: line( position( newage * coef ) ) END GPL.
Page 302 Figure 12.3 Delta beta measures to assess the influence of individual patterns on estimated coefficients
NOTE: We have skipped this graph for now.
Page 303 Table 12.2 Percent change in estimated parameters when including and excluding influential patterns
NOTE: We have skipped this graph for now.
Page 304 Table 12.3 Estimated probability of being a case (p-hat) for five influential observations
NOTE: We have skipped this graph for now.
Page 305 Figure 12.4 Delta chi-square measure to assess influence of pattern on overall fit with symbol size proportional to delta beta
NOTE: We have skipped this graph for now.
Page 307 Figure 12.5 Percentage of individuals correctly classified by logistic regression
NOTE: We have skipped this graph for now.
Page 307 Figure 12.6 ROC curve from logistic regression for the depression data set
NOTE: We have skipped this graph for now.