* For 590AT/593CA * Logistic regression * * SAS commands that create data set of the high school and beyond data; data hsb; input id sex race ses sctyp hsp locus concpt mot car rdg wrtg math sci civ; type='case'; achieve =( rdg+wrtg+math+sci+civ)/5; if hsp=2 then academic=1; else academic=0; ncases=1; label id ='ID number ' sex ='Gender: 1=male, 2=female' ses ='Socio-Economic Status' sctyp ='School Type: 1=public, 2=private' hsp = 'High School Program' locus ='Locus of Control' concpt ='Self Concept' mot ='Motivation' car ='Career Choice' rdg ='Reading T-Score' wrtg ='Writing T-Score' math ='Math T-Score' sci ='Science T-Score' civ ='Civics T-Score'; datalines; 1 2 1 1 1 3 0.29 0.88 0.67 10 33.6 43.7 40.2 39.0 40.6 2 1 1 1 1 1 -0.42 0.03 0.33 2 46.9 35.9 41.9 36.3 45.6 3 2 1 1 1 1 0.71 0.03 0.67 9 41.6 59.3 41.9 44.4 45.6 4 2 1 2 1 3 0.06 0.03 0.00 15 38.9 41.1 32.7 41.7 40.6 5 2 1 2 1 3 0.22 -0.28 0.00 1 36.3 48.9 39.5 41.7 45.6 6 1 1 2 1 1 0.46 0.03 0.00 11 49.5 46.3 46.2 41.7 35.6 7 1 1 1 1 2 0.44 -0.47 0.33 10 62.7 64.5 48.0 63.4 55.6 8 2 1 1 1 2 0.68 0.25 1.00 9 44.2 51.5 36.9 49.8 55.6 9 1 1 2 1 2 0.06 0.56 0.33 9 46.9 41.1 45.3 47.1 55.6 10 2 1 1 1 2 0.05 0.15 1.00 11 44.2 49.5 40.5 39.0 50.6 11 2 1 1 1 2 0.25 0.34 1.00 9 49.5 61.9 42.9 41.7 50.6 12 1 1 1 1 1 0.06 -0.59 0.00 8 44.2 41.1 59.8 44.4 35.6 13 1 1 1 1 2 -1.10 0.03 0.00 7 46.9 40.4 42.8 45.0 30.6 14 2 1 2 1 1 0.52 -0.59 0.67 4 41.6 59.3 38.6 36.3 50.6 15 1 1 2 1 3 -1.28 1.19 1.00 9 41.6 41.1 45.9 44.4 40.6 16 2 1 1 1 2 0.06 0.03 1.00 10 46.9 41.1 46.2 39.6 40.6 17 2 1 3 1 2 0.26 0.03 1.00 10 60.1 64.5 55.7 63.4 60.5 18 1 1 2 1 3 0.21 0.94 0.00 12 38.9 48.9 41.8 58.5 45.6 19 1 1 2 1 2 0.06 0.03 1.00 6 57.4 54.1 58.9 47.1 50.6 20 1 1 3 1 2 0.28 0.03 1.00 10 73.3 64.5 68.0 63.4 70.5 21 1 1 2 1 1 -0.44 -0.84 0.33 2 46.9 46.3 45.9 37.9 45.6 22 2 1 1 1 2 -0.43 0.13 0.33 7 38.9 44.3 51.6 43.9 48.1 23 2 1 1 1 3 -0.22 0.32 0.33 1 35.2 38.5 39.9 34.7 40.6 24 2 1 1 1 3 0.73 -1.42 0.33 15 38.9 41.1 41.0 36.3 45.6 25 1 1 2 1 3 0.26 0.03 0.33 11 48.4 48.9 52.2 43.9 50.6 26 2 1 2 1 2 1.12 -0.74 0.67 15 31.0 41.1 36.0 36.9 45.6 27 2 1 3 1 1 0.46 0.59 1.00 16 52.1 62.5 53.6 56.3 65.5 28 2 1 2 1 1 -0.60 -0.47 0.67 9 46.9 54.1 49.0 52.6 60.5 29 2 1 2 1 2 -0.90 0.03 0.67 4 36.3 44.3 36.1 33.6 40.6 30 2 1 2 1 2 0.31 1.19 1.00 9 46.9 44.3 41.9 39.0 50.6 31 1 1 2 1 2 -0.21 -1.38 0.00 9 34.2 46.3 44.5 39.0 35.6 32 2 1 2 1 2 0.69 0.34 1.00 9 52.1 56.7 53.4 60.7 45.6 33 2 1 2 1 1 -0.19 -1.67 0.00 13 42.1 39.8 42.2 38.5 40.6 34 1 1 2 1 3 -1.78 0.56 0.33 8 37.3 43.7 45.4 39.0 45.6 35 1 1 2 1 3 0.30 0.56 1.00 16 44.2 33.3 37.7 30.9 41.9 36 1 1 2 1 2 1.36 0.44 0.67 9 70.7 58.0 65.4 63.4 60.5 37 2 1 2 1 3 -1.05 -1.65 0.00 9 46.9 46.3 38.6 47.1 60.5 38 1 1 2 1 3 0.29 0.39 0.33 7 41.6 33.3 37.6 28.2 45.6 39 1 1 2 1 2 -0.83 0.13 1.00 9 41.6 35.9 36.7 36.3 55.6 40 1 1 3 1 2 0.21 -0.28 0.67 12 46.9 41.1 53.6 41.7 55.6 41 1 1 3 1 2 -0.03 0.63 0.67 9 52.1 54.1 48.2 55.3 60.5 42 1 1 2 1 2 0.03 0.28 0.33 8 52.1 54.1 52.0 55.3 50.6 43 1 1 3 1 3 -0.83 1.19 0.67 9 38.9 38.5 44.4 26.0 41.9 44 2 1 1 1 3 0.22 -0.76 1.00 15 31.0 41.1 49.2 33.6 35.6 45 1 1 3 1 2 1.16 0.32 1.00 9 60.1 54.1 48.6 58.0 30.6 46 1 1 1 1 3 -0.22 -1.34 0.67 16 46.9 31.3 43.6 36.3 35.6 47 1 1 1 1 1 -0.18 0.15 0.00 15 38.9 28.1 36.0 36.3 35.6 48 2 1 3 1 3 0.02 0.03 0.33 4 41.6 35.9 37.7 41.7 45.6 49 2 1 2 1 2 -1.06 0.03 0.67 6 46.9 56.7 48.0 44.4 40.6 50 1 1 2 1 1 -1.33 0.03 1.00 9 41.6 48.9 33.7 44.4 45.6 51 1 1 1 1 3 -1.55 0.03 0.00 9 36.3 41.1 43.5 33.6 40.6 52 1 1 2 1 3 -1.58 0.56 0.67 2 49.5 33.3 48.7 44.4 35.6 53 2 1 3 1 3 0.47 -0.09 0.33 9 33.6 33.9 38.8 39.6 30.6 54 1 1 3 1 3 0.46 0.03 0.67 9 54.8 54.7 56.9 58.0 40.6 55 2 1 1 1 1 -1.78 0.56 1.00 9 28.3 46.3 42.8 44.4 50.6 56 1 1 1 1 2 -0.84 -1.09 0.33 2 41.6 41.7 45.1 36.3 40.6 57 2 1 2 1 1 -0.89 -0.18 0.67 10 44.2 43.7 44.5 41.7 50.6 58 1 1 3 1 2 -0.68 -0.26 0.00 2 60.1 51.5 56.8 60.7 60.5 59 2 1 1 1 3 -0.60 0.32 1.00 1 38.9 56.7 41.2 33.6 55.6 60 1 1 2 1 3 -0.59 1.19 0.67 6 36.3 36.5 35.1 33.6 35.6 61 1 1 2 1 1 -0.04 0.65 1.00 16 44.2 44.3 61.1 49.8 45.6 62 2 1 1 1 1 0.52 -0.28 0.67 9 38.9 41.7 33.7 30.9 50.6 63 2 1 1 1 1 -0.34 0.59 1.00 7 38.9 33.9 35.1 44.4 40.6 64 1 1 2 1 3 -0.44 0.03 0.33 9 41.6 38.5 38.7 56.3 45.6 65 1 1 2 1 3 -0.21 0.03 1.00 6 41.6 43.7 51.1 49.8 40.6 66 2 2 3 1 2 0.32 -1.17 0.33 9 52.1 54.1 58.1 47.1 60.5 67 2 2 1 1 2 -0.71 -2.29 0.33 9 65.4 64.5 63.7 58.0 70.5 68 1 2 1 1 1 0.45 0.03 0.67 2 41.6 25.5 44.3 36.3 35.6 69 1 2 1 1 1 0.71 0.03 1.00 10 46.9 61.9 46.2 60.7 45.6 70 1 2 2 1 2 0.06 0.44 0.67 9 52.1 61.9 66.3 47.1 45.6 71 2 2 2 1 3 0.06 0.28 1.00 1 44.2 54.1 47.1 58.0 50.6 72 1 2 3 1 2 0.91 0.34 1.00 10 52.1 54.1 58.1 55.8 50.6 73 2 2 2 1 1 -1.09 0.03 0.33 2 46.9 43.7 41.9 41.7 35.6 74 2 2 3 1 2 0.02 0.03 0.67 10 65.4 56.7 62.3 58.5 55.6 75 1 2 3 2 2 -1.05 0.28 0.67 9 62.7 64.5 66.3 58.0 60.5 76 2 2 3 1 2 0.11 -1.42 1.00 16 60.1 59.3 62.1 60.7 50.6 77 1 2 2 1 3 0.89 0.32 0.67 10 46.9 44.3 48.7 53.1 55.6 78 1 2 1 1 3 -0.45 -0.64 0.33 9 65.4 56.7 55.2 52.6 55.6 79 1 2 2 1 2 0.46 -0.30 1.00 10 53.2 60.6 61.2 56.9 55.6 80 2 2 2 1 2 0.46 -1.17 0.00 9 68.0 62.5 65.4 60.7 60.5 81 2 2 3 1 2 0.28 0.94 1.00 10 54.8 56.0 43.5 49.8 50.6 82 2 2 2 1 1 0.06 0.90 0.67 9 38.9 52.8 54.4 49.8 40.6 83 1 2 2 1 2 0.06 0.94 1.00 7 52.1 51.5 57.9 60.7 35.6 84 1 2 3 1 2 -0.17 0.31 1.00 10 65.4 67.1 75.5 71.5 70.5 85 1 2 1 1 1 -0.38 0.03 0.00 2 52.1 43.7 49.4 55.3 40.6 86 2 2 3 1 2 0.23 -0.60 1.00 10 57.4 67.1 68.0 59.1 65.5 87 1 2 2 1 2 1.36 0.59 1.00 10 62.7 67.1 69.6 68.8 55.6 88 2 2 3 1 2 -0.39 1.19 1.00 9 40.5 59.3 41.9 33.6 50.6 89 1 2 2 1 1 -1.28 -1.05 0.33 9 57.4 43.7 57.8 60.7 45.6 90 2 2 1 1 3 -0.65 -1.17 0.00 8 36.3 30.7 45.9 36.3 40.6 91 2 2 2 2 1 -0.40 -0.47 1.00 10 54.8 59.3 51.8 41.7 55.6 92 2 2 3 1 3 0.06 -0.29 0.67 1 49.5 56.7 68.1 49.8 55.6 93 2 2 2 1 2 0.55 0.90 1.00 9 62.7 61.9 59.6 60.7 65.5 94 2 2 3 1 3 -0.35 -0.89 0.00 1 49.5 67.1 66.4 66.1 55.6 95 2 2 2 1 2 0.30 -1.09 0.00 9 51.1 56.7 55.3 47.7 56.8 96 2 2 1 1 3 -0.93 0.65 1.00 5 47.4 45.6 49.4 42.3 45.6 97 2 2 1 1 2 1.16 -0.72 0.67 10 56.9 64.5 72.2 54.2 55.6 98 1 2 3 1 2 0.00 0.15 0.33 10 73.3 64.5 75.5 60.7 50.6 99 1 2 2 1 2 0.52 0.34 1.00 9 52.1 54.1 51.4 49.8 50.6 100 2 1 3 2 2 0.42 -0.47 0.33 10 73.3 61.2 57.4 54.7 65.5 101 1 1 2 2 1 -0.67 -0.47 0.00 8 44.2 48.9 43.1 55.3 45.6 102 2 1 3 2 2 0.45 0.65 1.00 9 57.4 56.7 46.9 52.6 55.6 103 2 1 1 2 1 0.06 0.32 0.00 15 60.1 54.1 49.9 49.8 50.6 104 1 1 2 2 2 0.26 0.88 1.00 9 65.4 59.3 62.3 60.7 65.5 105 2 1 3 2 2 0.46 0.34 1.00 9 57.4 54.1 59.6 60.7 70.5 106 2 3 1 1 1 -0.89 -1.23 1.00 10 44.2 48.9 43.8 35.2 50.6 107 2 3 2 1 3 0.57 -2.62 0.00 11 38.9 28.1 38.4 44.4 35.6 108 2 3 3 1 3 0.26 0.94 0.00 17 54.8 38.5 46.8 36.3 45.6 109 2 3 1 1 3 0.31 1.19 1.00 1 40.5 46.9 40.4 39.0 50.6 110 1 3 1 1 3 0.93 -1.63 0.67 8 38.9 41.1 40.3 34.1 45.6 111 1 3 3 1 2 0.68 0.94 1.00 9 62.7 52.1 66.1 68.8 50.6 112 1 3 1 1 2 0.03 0.32 1.00 6 44.7 56.7 50.3 30.9 55.6 113 1 3 3 1 2 -0.03 0.87 1.00 7 53.2 46.3 43.0 41.7 50.6 114 1 3 2 1 2 0.24 -0.35 1.00 8 30.5 35.9 36.9 33.6 40.6 115 2 3 3 1 2 1.14 1.19 1.00 9 65.9 67.1 67.1 60.7 65.5 116 2 3 2 1 2 0.27 0.34 1.00 7 49.5 48.9 48.0 41.7 60.5 117 2 3 1 1 3 -1.34 1.19 1.00 9 44.2 41.1 44.2 44.4 40.6 118 1 3 1 1 1 0.45 -0.91 0.33 2 41.6 41.1 42.8 49.8 40.6 119 2 3 1 1 2 0.28 1.19 1.00 1 56.4 59.3 38.5 55.3 56.8 120 2 3 1 1 3 -1.33 -0.26 0.67 15 38.4 37.2 33.4 30.9 36.9 121 1 3 2 1 2 0.49 0.94 1.00 9 49.5 39.8 44.5 55.3 55.6 122 2 3 1 1 2 0.10 0.32 0.67 9 40.0 40.4 42.4 48.8 60.5 123 2 3 2 1 1 1.11 -1.05 0.33 9 70.7 67.1 63.0 63.4 65.5 124 2 3 2 1 3 0.22 -0.76 0.67 16 45.8 52.1 54.9 44.4 55.6 125 2 3 3 1 1 0.75 -0.72 0.67 14 42.6 46.3 41.2 28.2 45.6 126 1 3 2 1 2 0.89 0.15 0.67 10 68.0 59.3 56.6 50.4 60.5 127 2 3 1 1 2 -0.49 0.03 0.00 8 46.9 37.2 42.8 41.7 45.6 128 1 3 2 1 2 0.49 0.44 0.67 10 53.7 43.7 51.1 52.0 43.1 129 1 3 2 1 3 -0.38 0.34 0.67 7 38.9 28.1 35.3 39.0 40.6 130 2 3 1 1 3 -0.60 0.90 0.67 1 46.9 61.9 45.4 33.6 45.6 131 2 3 2 1 2 1.36 0.65 0.67 10 59.0 56.7 52.9 48.8 56.8 132 2 3 1 1 2 -0.24 -0.16 0.33 16 62.7 61.9 52.4 47.1 60.5 133 2 3 1 1 3 -0.23 0.44 0.33 14 33.6 35.2 40.9 29.3 25.7 134 2 3 2 1 2 0.51 -1.29 0.33 1 39.4 41.7 46.9 35.2 50.6 135 1 3 2 1 3 -0.60 -1.43 1.00 10 38.9 51.5 38.6 39.0 35.6 136 2 3 1 1 2 0.06 0.94 0.67 9 45.3 54.7 44.3 33.6 40.6 137 1 3 2 1 2 -0.37 -1.90 0.67 6 54.8 36.5 37.7 49.8 60.5 138 1 3 2 1 2 0.46 0.94 0.33 5 36.3 30.7 44.8 39.6 40.6 139 2 3 1 1 2 0.44 0.65 1.00 8 46.9 46.3 48.5 32.5 40.6 140 2 3 1 1 3 0.10 1.19 1.00 15 38.9 35.9 46.8 39.0 40.6 141 1 3 2 1 3 -2.23 1.19 1.00 8 36.3 38.5 39.3 39.0 45.6 142 1 3 2 1 2 1.13 0.87 1.00 10 57.4 54.7 51.7 49.8 50.6 143 1 3 2 1 2 0.26 1.19 1.00 9 52.1 52.8 37.6 52.6 50.6 144 2 3 1 1 1 -0.04 -0.29 0.67 10 41.6 38.5 40.2 33.6 45.6 145 1 3 3 1 3 0.69 0.03 0.67 9 49.5 39.8 38.6 49.3 46.8 146 2 3 1 1 1 -1.76 -0.93 0.67 6 36.3 38.5 36.4 36.3 50.6 147 1 3 3 1 3 -0.66 -0.08 1.00 12 44.7 33.3 33.7 39.0 35.6 148 1 3 2 1 1 -0.44 -0.31 0.67 9 49.5 59.3 42.1 52.6 60.5 149 1 3 1 1 2 1.16 1.19 0.33 9 55.3 46.9 49.1 54.2 51.8 150 2 3 2 1 2 0.91 -0.28 1.00 7 60.1 67.1 56.2 37.4 50.6 151 2 3 3 1 1 0.75 -0.26 0.67 15 41.6 35.9 42.0 31.4 39.4 152 1 3 2 1 2 -1.14 -0.06 1.00 6 36.3 46.3 36.9 41.7 35.6 153 2 3 2 1 1 0.46 0.03 1.00 14 52.1 56.7 62.8 47.1 60.5 154 2 3 1 1 2 0.45 0.03 1.00 1 49.5 46.3 52.9 52.6 65.5 155 2 3 1 1 3 -0.54 -0.59 1.00 4 38.9 35.9 38.4 33.6 35.6 156 1 3 3 1 3 -0.42 -0.60 1.00 14 54.8 56.7 64.7 58.0 60.5 157 1 3 2 1 3 -0.38 0.56 1.00 2 33.6 36.5 45.7 39.0 30.6 158 2 3 1 1 1 1.16 0.59 1.00 10 36.3 45.0 44.3 39.0 40.6 159 1 3 1 1 1 -0.68 0.32 0.67 9 36.3 44.3 34.4 47.1 45.6 160 2 3 1 2 1 0.70 0.87 1.00 1 46.9 54.1 46.4 49.8 55.6 161 2 3 2 2 2 0.48 0.63 1.00 10 54.8 48.2 58.7 51.5 45.6 162 2 3 3 2 2 -0.40 -1.80 0.33 1 50.0 45.0 42.9 41.7 45.6 163 2 3 2 2 2 0.66 0.00 1.00 10 52.1 48.9 48.5 44.4 60.5 164 1 4 3 1 2 0.06 0.03 1.00 9 62.7 64.5 71.3 55.3 70.5 165 1 4 1 1 3 0.46 -1.17 0.33 5 46.9 30.7 40.3 58.0 35.6 166 1 4 2 1 3 -0.40 0.88 0.33 3 54.8 45.0 45.9 58.0 50.6 167 2 4 3 1 2 0.06 0.03 0.67 10 54.8 54.1 47.1 49.8 50.6 168 2 4 2 1 1 -0.69 -1.13 0.67 1 46.9 54.1 45.4 44.4 45.6 169 2 4 2 1 2 1.16 1.19 1.00 10 70.7 64.5 72.2 66.1 55.6 170 1 4 1 1 1 0.30 -0.59 1.00 16 54.8 46.3 45.5 58.0 60.5 171 2 4 2 1 2 0.08 0.59 1.00 9 57.4 56.7 53.9 60.7 60.5 172 1 4 2 1 3 -1.30 0.03 1.00 7 54.8 41.1 40.2 44.4 40.6 173 2 4 1 1 1 -0.61 0.03 0.33 11 44.2 54.1 40.3 52.6 40.6 174 1 4 2 1 2 -0.16 0.65 1.00 7 68.0 61.9 69.7 71.5 60.5 175 2 4 2 1 2 1.16 0.65 1.00 1 65.4 67.1 63.1 55.3 70.5 176 2 4 3 1 1 0.91 0.65 0.33 11 54.8 48.9 41.9 52.6 55.6 177 1 4 2 1 2 0.68 0.03 0.67 10 71.2 63.2 60.2 65.5 65.5 178 1 4 2 1 2 0.06 -0.29 0.67 3 57.4 64.5 51.1 63.4 60.5 179 2 4 2 1 2 1.16 0.03 1.00 6 70.7 67.1 61.2 60.7 70.5 180 2 4 2 1 3 -0.59 0.03 0.33 14 33.6 54.1 41.0 41.7 50.6 181 2 4 3 1 2 1.13 0.87 1.00 10 76.0 62.5 60.0 67.2 65.5 182 1 4 3 1 2 0.42 0.56 0.67 17 70.7 56.0 62.0 67.7 63.0 183 2 4 1 1 2 1.16 0.63 0.67 13 58.5 54.1 46.3 48.2 53.1 184 1 4 3 1 1 1.16 -0.81 0.67 9 65.4 64.5 48.0 63.4 65.5 185 2 4 2 1 1 0.45 0.56 1.00 5 52.1 43.7 41.9 47.1 45.6 186 1 4 2 1 3 -0.16 -0.53 0.67 5 62.7 61.9 67.0 66.1 45.6 187 2 4 1 1 1 0.68 -0.60 0.33 14 52.1 64.5 59.7 55.8 50.6 188 2 4 3 1 2 0.06 0.03 1.00 10 41.6 54.1 41.2 41.7 40.6 189 2 4 3 1 2 1.11 0.90 0.33 6 62.7 64.5 61.4 58.0 60.5 190 2 4 3 1 3 -0.41 0.03 0.67 15 49.5 51.5 45.1 58.0 35.6 191 1 4 3 1 2 0.49 0.94 1.00 16 54.8 56.7 63.8 60.7 65.5 192 2 4 2 1 2 0.67 0.03 0.67 9 52.1 56.7 51.1 55.3 55.6 193 2 4 2 1 2 0.51 0.03 1.00 14 54.8 54.1 66.4 41.7 55.6 194 1 4 3 1 1 -0.84 0.34 0.00 10 41.6 46.3 59.6 63.4 65.5 195 1 4 3 1 2 0.27 0.03 1.00 10 60.1 51.5 48.1 52.6 55.6 196 2 4 2 1 3 0.06 0.03 1.00 1 68.0 61.9 55.5 49.8 50.6 197 2 4 1 1 3 -1.15 -1.38 1.00 17 49.5 56.7 46.1 53.6 40.6 198 2 4 2 1 2 0.66 -0.47 0.67 9 54.8 59.3 50.6 60.7 50.6 199 1 4 1 1 3 -0.82 0.63 1.00 12 36.8 36.5 41.5 33.1 31.9 200 1 4 3 1 2 -0.27 0.88 1.00 13 52.1 64.5 60.6 60.7 45.6 201 2 4 2 1 3 0.31 0.63 1.00 1 38.9 56.7 46.3 58.0 55.6 202 1 4 2 1 2 1.13 -0.60 0.67 16 73.3 67.1 71.3 63.4 65.5 203 2 4 2 1 2 -0.44 -0.78 1.00 10 58.0 48.9 43.8 44.4 65.5 204 2 4 2 1 2 0.44 0.03 1.00 9 55.3 51.5 48.0 58.0 65.5 205 2 4 1 1 3 -0.90 -2.54 0.00 4 44.2 43.7 40.2 39.6 30.6 206 2 4 1 1 3 -0.18 0.28 0.33 4 44.2 51.5 46.1 47.1 50.6 207 2 4 2 1 1 1.11 0.90 0.33 4 55.3 50.2 41.7 58.5 55.6 208 1 4 1 1 1 0.05 -0.60 0.00 8 57.4 51.5 40.6 46.6 56.8 209 2 4 2 1 1 -0.63 0.44 0.33 4 41.6 43.7 46.8 36.3 50.6 210 2 4 2 1 1 0.06 0.03 1.00 9 54.8 59.3 46.2 47.1 65.5 211 2 4 2 1 1 0.71 -0.29 0.67 9 57.4 61.9 55.5 58.0 65.5 212 1 4 2 1 1 0.26 0.03 1.00 2 57.4 51.5 55.3 60.7 60.5 213 2 4 1 1 1 0.04 -0.47 0.67 1 33.6 59.3 44.7 47.1 45.6 214 2 4 2 1 3 -1.81 -0.64 1.00 15 41.6 54.1 46.9 47.1 45.6 215 2 4 2 1 2 -0.40 -0.22 0.67 9 52.1 38.5 36.1 44.4 45.6 216 2 4 3 1 2 0.68 0.03 1.00 13 65.4 67.1 64.5 60.7 70.5 217 2 4 2 1 2 0.00 -0.42 0.67 11 49.5 51.5 52.9 39.0 55.6 218 2 4 2 1 3 0.46 0.03 1.00 1 57.4 67.1 47.7 55.8 50.6 219 2 4 2 1 3 0.10 0.03 0.33 6 49.5 56.7 48.0 47.1 60.5 220 2 4 2 1 2 0.66 0.34 1.00 4 56.9 49.5 49.5 50.9 58.1 221 1 4 1 1 1 -0.36 -1.67 1.00 3 44.2 43.7 56.4 58.0 60.5 222 1 4 3 1 1 0.93 0.34 1.00 10 46.9 54.1 54.6 55.3 50.6 223 1 4 2 1 3 0.48 0.03 1.00 10 60.1 46.3 51.4 52.6 60.5 224 2 4 2 1 1 -1.15 -0.47 0.67 4 44.2 56.7 52.2 44.4 55.6 225 1 4 2 1 2 0.89 0.59 0.67 6 60.6 56.7 70.5 58.0 65.5 226 1 4 3 1 2 0.49 1.19 1.00 10 46.9 51.5 50.6 49.8 55.6 227 1 4 2 1 2 0.44 0.65 1.00 3 52.1 54.1 54.6 41.7 55.6 228 2 4 2 1 2 0.00 -0.82 0.67 10 54.8 67.1 63.8 49.8 50.6 229 2 4 1 1 2 0.91 0.03 1.00 1 60.6 59.3 49.4 44.4 65.5 230 2 4 3 1 1 -0.20 1.19 1.00 15 41.6 46.9 41.7 47.1 60.5 231 2 4 1 1 1 0.00 -1.09 0.67 9 62.7 51.5 45.9 47.1 45.6 232 2 4 2 1 2 0.20 -0.47 0.33 11 38.9 54.1 53.7 52.6 40.6 233 2 4 2 1 3 0.06 0.03 0.67 15 62.7 61.9 58.8 60.7 60.5 234 2 4 3 1 2 0.96 0.87 1.00 9 57.4 64.5 48.0 47.1 50.6 235 2 4 2 1 2 0.20 -0.55 0.67 9 60.1 61.9 48.7 49.8 50.6 236 1 4 2 1 2 0.08 0.94 0.67 10 73.3 61.9 60.5 71.5 55.6 237 1 4 2 1 2 1.36 -0.80 0.67 13 52.1 59.3 59.6 59.1 40.6 238 1 4 3 1 2 0.46 0.94 1.00 10 65.4 61.9 67.9 66.1 65.5 239 2 4 1 1 3 -0.84 -0.18 0.00 1 33.6 46.3 38.4 36.3 45.6 240 1 4 2 1 2 0.46 0.03 0.67 9 62.7 54.1 54.8 55.3 65.5 241 1 4 1 1 2 0.48 0.69 0.67 16 62.7 43.0 58.9 65.0 44.3 242 1 4 3 1 2 0.75 1.19 1.00 10 68.0 61.9 52.2 60.7 65.5 243 2 4 2 1 2 -0.40 0.03 0.00 1 54.8 49.5 55.3 47.1 51.8 244 2 4 3 1 2 0.71 0.34 0.67 9 68.0 61.9 64.5 68.8 60.5 245 1 4 2 1 1 0.46 0.34 0.67 9 62.2 56.7 51.9 54.7 58.1 246 2 4 2 1 3 0.43 -0.47 1.00 1 42.6 50.8 44.6 39.0 55.6 247 2 4 2 1 3 0.46 -0.47 0.00 9 49.5 61.9 41.4 55.3 30.6 248 1 4 2 1 1 -0.49 0.03 0.33 9 65.4 59.3 56.8 66.1 60.5 249 2 4 2 1 2 -1.56 0.03 0.33 15 62.7 54.1 53.0 55.3 45.6 250 1 4 2 1 1 -0.17 -0.84 0.67 15 62.7 56.7 54.3 58.0 50.6 251 1 4 3 1 1 0.66 0.34 1.00 7 60.1 51.5 53.0 63.4 55.6 252 2 4 1 1 2 0.48 0.32 1.00 14 49.5 54.1 38.7 49.8 50.6 253 1 4 2 1 1 -0.14 -0.24 0.67 9 54.8 38.5 57.0 52.6 45.6 254 2 4 3 1 2 0.48 -0.47 0.33 9 52.1 61.9 55.5 60.7 60.5 255 1 4 2 1 2 -1.13 -0.55 0.33 12 44.2 46.9 45.5 39.0 50.6 256 1 4 3 1 1 0.06 0.03 0.00 10 44.2 33.3 53.9 58.0 30.6 257 1 4 2 1 1 0.20 -0.47 0.33 3 38.9 38.5 53.0 52.6 50.6 258 2 4 2 1 3 -0.64 -0.51 0.67 4 49.5 56.7 47.7 44.4 60.5 259 2 4 2 1 1 0.23 0.44 1.00 2 49.5 51.5 45.5 49.8 55.6 260 2 4 2 1 2 0.10 1.19 1.00 9 52.1 56.7 62.3 66.1 70.5 261 1 4 2 1 3 0.28 0.03 1.00 9 54.8 59.3 49.5 63.4 56.8 262 2 4 3 1 2 0.52 0.59 1.00 9 68.0 59.9 64.0 68.8 65.5 263 1 4 2 1 2 0.46 0.03 0.67 8 54.8 64.5 58.6 55.3 60.5 264 2 4 1 1 1 -0.80 0.15 0.33 15 41.6 41.1 39.5 47.1 60.5 265 2 4 1 1 3 -1.33 -0.60 0.33 1 34.7 35.2 40.2 50.9 33.1 266 2 4 1 1 1 -1.33 0.34 0.00 4 44.2 41.1 36.9 41.7 45.6 267 1 4 3 1 2 0.68 0.34 1.00 16 57.4 56.7 62.8 58.0 50.6 268 2 4 3 1 2 0.46 -0.55 0.33 10 41.6 54.1 50.3 49.8 51.8 269 1 4 2 1 3 -0.44 -0.22 0.00 2 44.2 33.3 40.9 47.1 45.6 270 2 4 2 1 2 0.23 -0.76 1.00 9 65.4 64.5 59.6 66.1 60.5 271 2 4 3 1 3 0.29 0.03 1.00 4 49.5 48.9 56.2 47.1 45.6 272 2 4 1 1 1 0.46 0.03 0.67 16 52.1 56.7 53.0 47.1 50.6 273 2 4 2 1 1 0.46 0.03 0.33 13 44.7 38.5 45.9 44.4 35.6 274 2 4 3 1 1 0.33 -0.26 1.00 4 52.1 67.1 57.0 63.4 60.5 275 2 4 3 1 1 -0.44 -0.47 1.00 9 54.8 65.1 66.1 49.8 60.5 276 1 4 2 1 3 -1.05 -0.34 0.67 1 44.2 38.5 41.4 42.3 35.6 277 2 4 3 1 2 1.11 0.34 1.00 10 73.3 67.1 62.3 58.0 65.5 278 1 4 1 1 2 0.46 -0.18 1.00 10 52.1 46.3 48.3 58.0 60.5 279 1 4 1 1 3 0.08 0.94 0.33 2 49.5 52.8 50.6 48.8 66.8 280 1 4 3 1 2 0.31 0.03 1.00 10 54.8 48.9 61.2 60.7 55.6 281 1 4 3 1 3 -0.18 0.03 0.67 2 44.2 48.9 56.3 41.7 60.5 282 2 4 2 1 1 0.71 -0.05 0.67 1 52.7 47.6 52.2 45.0 50.6 283 1 4 3 1 2 0.02 0.13 0.33 9 73.3 59.9 70.5 60.7 70.5 284 1 4 3 1 2 1.36 1.19 1.00 9 65.4 48.9 66.3 58.0 35.6 285 2 4 2 1 2 0.00 0.03 0.67 1 61.1 61.2 62.2 57.4 53.1 286 2 4 3 1 2 0.52 0.34 1.00 10 65.4 54.1 61.4 58.0 55.6 287 2 4 3 1 2 0.66 0.34 0.67 9 54.8 56.7 61.9 63.4 55.6 288 1 4 1 1 1 0.25 0.65 0.00 3 49.5 48.9 50.4 63.4 55.6 289 1 4 3 1 2 -0.38 0.03 1.00 9 60.1 54.1 58.0 58.0 60.5 290 2 4 3 1 1 0.46 -1.11 0.33 9 65.4 67.1 67.0 66.1 70.5 291 2 4 2 1 1 0.23 0.03 1.00 13 65.4 56.7 65.4 58.0 70.5 292 2 4 1 1 3 -1.09 -0.90 0.67 9 57.4 59.9 50.5 52.6 36.9 293 2 4 2 1 2 0.71 -0.55 1.00 9 60.1 59.3 64.3 60.7 65.5 294 1 4 1 1 2 -0.82 -0.76 0.00 1 57.4 43.7 59.6 52.6 50.6 295 2 4 3 1 1 1.36 0.03 0.67 14 46.9 59.3 56.3 66.1 60.5 296 2 4 3 1 3 -0.14 -0.86 0.00 9 38.9 46.3 50.1 36.3 45.6 297 2 4 3 1 2 0.46 1.19 0.67 10 65.9 64.5 61.8 60.7 60.5 298 2 4 3 1 2 0.06 0.32 1.00 1 62.7 64.5 70.5 68.8 70.5 299 1 4 3 1 3 -0.51 -1.98 0.00 17 60.1 57.3 58.8 58.0 50.6 300 1 4 3 1 2 0.68 0.03 0.33 9 73.3 64.5 64.0 60.7 65.5 301 2 4 3 1 2 0.75 1.19 1.00 9 60.1 61.9 67.1 49.8 55.6 302 1 4 2 1 2 0.69 0.05 0.67 10 60.1 59.3 57.9 40.6 65.5 303 2 4 1 1 1 -0.84 -0.24 1.00 4 54.8 64.5 44.5 52.6 50.6 304 1 4 3 1 2 0.01 0.32 0.00 16 52.1 41.1 50.6 52.6 55.6 305 2 4 2 1 1 -0.16 -0.45 0.67 4 49.5 54.1 61.9 52.6 50.6 306 2 4 3 1 2 -0.60 0.28 1.00 9 33.6 41.1 42.9 41.7 50.6 307 1 4 3 1 2 0.46 0.34 1.00 9 76.0 52.1 64.1 63.9 60.5 308 1 4 1 1 2 0.00 0.65 0.67 7 46.9 33.3 50.4 47.1 55.6 309 2 4 3 1 1 0.32 0.90 0.67 10 33.6 51.5 41.9 49.8 50.6 310 1 4 3 1 2 0.52 0.65 0.33 2 54.3 62.5 56.6 54.7 45.6 311 1 4 3 1 3 0.46 0.63 0.67 2 46.9 52.8 49.3 53.1 35.6 312 2 4 2 1 2 0.66 0.03 1.00 1 52.1 56.7 41.9 52.6 60.5 313 2 4 2 1 2 0.46 -0.47 0.67 10 50.0 41.1 45.1 44.4 55.6 314 1 4 2 1 1 0.26 0.03 1.00 2 49.5 35.9 47.8 52.6 35.6 315 2 4 3 1 2 0.91 -0.47 0.67 10 57.4 59.3 57.2 59.1 65.5 316 2 4 2 1 3 1.36 0.87 1.00 4 36.3 43.7 37.2 41.7 40.6 317 2 4 2 1 2 0.02 -1.05 0.33 9 60.1 64.5 61.2 60.7 50.6 318 2 4 1 1 1 -0.44 -0.47 0.33 9 38.9 46.3 41.9 41.7 45.6 319 1 4 1 1 3 -0.44 -0.55 0.33 4 46.9 38.5 47.1 41.7 25.7 320 1 4 1 1 3 0.22 0.03 0.67 16 60.1 41.7 52.0 52.6 50.6 321 1 4 1 1 3 -0.91 -0.59 0.67 8 44.2 39.1 46.3 36.3 40.6 322 1 4 2 1 1 -0.84 -0.57 0.33 2 33.6 33.3 41.0 36.3 35.6 323 1 4 2 1 3 0.66 0.34 0.33 2 33.6 38.5 42.9 47.1 45.6 324 2 4 3 1 2 0.94 0.03 1.00 10 68.0 56.7 59.6 58.0 25.7 325 2 4 2 1 1 -0.65 0.03 0.33 9 41.6 39.1 42.3 41.7 40.6 326 1 4 1 1 1 1.13 0.03 0.33 5 45.8 41.7 43.1 53.6 36.9 327 2 4 1 1 3 1.36 0.63 1.00 1 68.0 64.5 58.9 63.4 45.6 328 2 4 2 1 3 0.22 0.03 0.33 1 52.1 54.7 49.5 53.6 60.5 329 1 4 3 1 2 -0.43 0.03 1.00 9 52.1 44.3 53.1 58.0 45.6 330 1 4 1 1 1 0.28 -0.06 0.33 2 44.2 30.7 35.3 47.1 45.6 331 2 4 1 1 1 -1.99 0.03 0.00 9 39.4 54.1 38.7 47.1 35.6 332 1 4 2 1 2 -0.40 0.03 0.00 10 44.2 33.3 47.8 44.4 45.6 333 1 4 2 1 1 -0.03 -0.16 0.00 6 33.6 43.7 42.6 47.1 35.6 334 2 4 2 1 2 0.08 -0.24 0.33 10 52.1 59.3 48.0 55.3 60.5 335 2 4 2 1 1 0.46 0.34 0.67 9 62.7 61.9 52.9 44.4 40.6 336 1 4 3 1 1 0.06 0.03 0.67 11 41.6 39.1 56.3 45.0 50.6 337 1 4 2 1 2 -0.39 -0.28 0.67 9 44.2 51.5 51.1 63.4 60.5 338 2 4 1 1 3 0.26 -0.26 0.00 6 49.5 56.7 43.8 47.1 40.6 339 1 4 2 1 2 0.46 0.03 1.00 7 70.7 56.7 51.3 71.5 55.6 340 1 4 3 1 2 0.47 0.01 1.00 10 68.0 64.5 61.5 55.3 60.5 341 2 4 3 1 2 0.71 0.03 1.00 14 68.0 61.9 72.2 66.1 56.8 342 2 4 2 1 2 -0.16 -0.47 0.67 9 36.3 61.9 52.2 39.0 45.6 343 1 4 1 1 1 -0.37 -0.47 0.67 9 41.6 48.9 42.6 49.8 55.6 344 2 4 1 1 1 0.27 -2.52 0.33 8 38.9 32.0 35.3 46.6 51.8 345 1 4 2 1 2 0.71 0.34 0.67 9 52.1 41.1 51.3 60.7 65.5 346 2 4 2 1 2 0.46 0.03 0.67 9 57.4 59.3 53.9 49.8 55.6 347 2 4 3 1 2 0.28 0.03 0.67 9 62.7 67.1 65.6 61.8 65.5 348 2 4 2 1 1 0.46 0.03 1.00 10 65.4 51.5 61.4 60.7 50.6 349 1 4 3 1 3 -1.50 0.03 0.67 9 33.6 48.9 38.6 42.3 55.6 350 1 4 3 1 2 -0.41 0.28 0.33 16 46.9 35.9 42.6 58.0 55.6 351 2 4 2 1 2 0.32 0.90 0.67 1 52.1 59.3 58.1 47.1 45.6 352 2 4 2 1 1 0.32 0.65 1.00 9 54.8 61.9 58.1 58.0 60.5 353 1 4 2 1 2 -0.21 -1.13 0.00 2 38.9 41.1 43.6 55.3 45.6 354 2 4 2 1 2 -0.60 -1.18 0.67 1 54.8 59.3 68.0 49.3 65.5 355 2 4 1 1 1 -0.40 -0.89 0.33 8 41.6 56.7 45.0 50.4 43.1 356 2 4 1 1 3 -0.89 -1.29 0.33 9 60.1 64.5 41.0 52.6 50.6 357 2 4 2 1 3 0.06 -0.28 1.00 1 65.9 64.5 54.6 60.7 60.5 358 2 4 3 1 2 0.45 0.03 1.00 4 62.7 51.5 54.4 49.8 50.6 359 1 4 2 1 2 0.22 0.32 1.00 14 49.5 56.7 62.8 49.8 70.5 360 2 4 3 1 2 0.71 0.28 1.00 11 60.1 64.5 63.0 68.8 60.5 361 2 4 2 1 3 0.91 0.34 1.00 1 68.0 59.3 53.0 63.4 60.5 362 1 4 2 1 1 0.00 0.65 1.00 6 54.8 64.5 51.1 49.8 55.6 363 1 4 1 1 3 -0.38 0.03 1.00 8 38.9 41.7 37.7 44.4 45.6 364 2 4 2 1 2 -0.16 -2.58 1.00 10 52.1 59.3 58.1 52.6 65.5 365 2 4 1 1 1 -1.10 -0.28 0.67 4 38.9 41.1 45.9 47.1 45.6 366 2 4 2 1 1 0.93 0.03 1.00 10 65.9 67.1 70.5 52.6 65.5 367 1 4 3 1 1 0.70 -0.16 0.33 9 68.0 59.3 55.7 63.4 65.5 368 1 4 2 1 3 -0.45 -0.26 0.00 2 38.9 35.9 44.2 55.3 55.6 369 2 4 2 1 3 0.02 -2.03 0.67 1 38.9 59.3 39.5 39.0 50.6 370 2 4 1 1 3 -0.19 -2.58 0.00 1 41.6 54.1 41.0 41.7 40.6 371 1 4 2 1 1 0.91 -0.47 0.67 8 44.2 41.1 50.5 58.0 40.6 372 1 4 2 1 1 0.09 -1.42 0.33 11 52.1 61.9 59.6 58.0 55.6 373 2 4 1 1 2 0.00 0.65 1.00 16 68.0 64.5 58.3 58.5 55.6 374 1 4 1 1 2 0.70 -0.16 1.00 10 49.5 61.9 50.9 46.1 66.8 375 1 4 1 1 1 -0.60 0.56 0.33 2 36.3 43.7 44.3 47.1 30.6 376 1 4 2 1 1 -1.33 0.03 0.67 2 41.6 30.7 56.9 47.1 50.6 377 2 4 2 1 2 -0.16 0.03 1.00 8 57.4 57.3 64.0 58.0 65.5 378 1 4 3 1 2 0.46 0.65 1.00 9 68.0 54.1 63.0 63.4 60.5 379 1 4 3 1 2 -0.33 0.38 0.67 6 62.7 59.3 56.5 55.3 55.6 380 2 4 3 1 1 0.25 0.03 0.67 4 49.5 51.5 55.5 44.4 55.6 381 2 4 2 1 2 0.22 -0.57 1.00 1 61.1 67.1 64.5 50.4 55.6 382 1 4 3 1 2 0.22 1.19 0.33 13 38.9 33.3 37.7 47.1 40.6 383 2 4 2 1 1 0.27 0.03 0.67 1 46.9 56.7 63.8 55.3 55.6 384 1 4 2 1 3 -2.23 1.19 0.33 5 44.2 38.5 41.2 49.8 45.6 385 2 4 1 1 1 -0.38 0.37 0.67 4 44.2 43.7 46.1 47.1 50.6 386 2 4 2 1 2 0.06 0.65 0.33 11 47.4 64.5 53.9 55.8 50.6 387 2 4 3 1 2 0.89 0.88 1.00 9 52.1 59.3 64.0 61.2 55.6 388 2 4 3 1 1 1.13 0.03 0.33 9 43.2 54.1 54.8 55.3 45.6 389 2 4 1 1 1 -1.02 -1.67 0.67 5 45.3 43.0 52.3 60.7 56.8 390 2 4 3 1 2 0.45 0.03 0.67 9 60.1 61.9 51.9 53.1 58.1 391 2 4 3 1 2 0.06 0.65 0.00 9 65.4 59.3 56.8 46.1 65.5 392 1 4 1 1 3 -1.74 0.61 0.33 17 45.8 34.6 37.9 39.0 45.6 393 1 4 2 1 3 -1.30 0.13 0.33 11 44.2 41.1 51.8 47.1 40.6 394 1 4 2 1 2 0.46 0.03 0.00 8 73.3 61.9 73.1 68.8 65.5 395 2 4 2 1 2 0.55 -0.24 0.67 13 46.9 56.7 54.6 47.1 55.6 396 2 4 2 1 3 -0.68 0.03 0.33 17 38.9 38.5 41.2 39.0 35.6 397 1 4 2 1 3 0.25 -0.26 0.33 2 49.5 30.7 40.3 33.6 30.6 398 2 4 1 1 3 0.96 0.94 0.67 9 46.9 50.2 46.5 47.1 51.8 399 2 4 1 1 1 -0.24 0.03 0.67 9 45.8 54.1 49.5 47.1 55.6 400 1 4 1 1 1 -0.40 0.65 0.67 5 44.2 43.7 38.6 33.6 45.6 401 2 4 1 1 2 -0.40 0.03 1.00 15 44.2 54.1 59.3 58.0 50.6 402 2 4 2 1 2 0.93 0.65 1.00 6 62.7 59.3 71.3 68.8 60.5 403 2 4 1 1 2 0.00 0.34 1.00 15 62.7 59.9 65.4 53.6 65.5 404 2 4 2 1 1 -0.38 -0.47 0.67 4 62.7 43.7 44.7 52.6 41.9 405 1 4 2 1 2 0.20 0.56 0.67 12 60.1 52.8 51.5 58.0 50.6 406 1 4 2 1 2 0.73 0.03 1.00 16 65.4 59.3 69.6 63.4 50.6 407 2 4 2 1 3 0.06 0.03 1.00 4 43.2 47.6 45.6 47.1 45.6 408 1 4 3 1 2 -0.45 -0.60 0.67 10 52.1 59.9 49.5 48.8 60.5 409 2 4 3 1 2 1.13 0.03 0.67 9 62.7 64.5 64.5 52.6 60.5 410 2 4 1 1 2 -0.18 -1.09 0.67 15 62.7 56.7 58.8 71.5 60.5 411 1 4 3 1 2 0.32 0.28 0.67 16 62.7 63.2 68.9 74.2 70.5 412 2 4 2 1 3 0.66 -0.60 0.00 15 42.6 56.7 40.3 49.8 50.6 413 2 4 2 1 1 0.51 0.03 0.33 6 60.1 51.5 53.9 63.4 50.6 414 2 4 2 1 1 0.08 0.03 1.00 13 54.8 67.1 47.4 49.8 65.5 415 2 4 2 1 2 0.46 -0.86 1.00 14 68.0 59.3 60.5 55.3 70.5 416 1 4 1 1 2 1.11 -0.09 0.00 9 44.2 48.9 48.0 49.8 55.6 417 2 4 2 1 1 0.07 0.34 1.00 9 65.4 59.3 51.2 60.7 55.6 418 1 4 2 1 3 -0.89 0.56 0.33 2 44.2 41.1 40.3 49.8 25.7 419 1 4 2 1 3 -0.19 0.03 0.33 2 54.8 51.5 42.8 60.7 50.6 420 1 4 1 1 2 -0.27 0.94 0.00 17 44.7 54.1 50.4 45.0 50.6 421 1 4 3 1 3 -0.40 -1.09 0.67 7 62.7 43.7 47.1 52.6 55.6 422 2 4 2 1 2 0.00 0.03 0.67 9 62.7 64.5 38.3 55.8 55.6 423 1 4 2 1 2 0.46 -0.47 0.67 10 68.0 64.5 70.6 68.8 59.3 424 2 4 2 1 3 0.75 0.34 1.00 9 47.4 41.7 51.6 38.5 50.6 425 2 4 1 1 1 1.36 0.65 1.00 15 46.9 51.5 57.0 49.8 40.6 426 1 4 1 1 2 0.68 0.87 1.00 10 54.8 56.7 56.9 58.0 70.5 427 1 4 2 1 3 0.05 -0.88 0.67 2 62.7 63.2 74.6 71.5 65.5 428 2 4 2 1 3 -0.84 0.03 0.00 9 36.3 59.3 47.7 39.0 55.6 429 1 4 3 1 3 1.13 0.56 1.00 6 60.1 48.9 50.4 55.3 65.5 430 1 4 3 1 1 1.13 0.63 0.67 10 60.1 64.5 57.9 60.7 65.5 431 2 4 2 1 1 0.55 -0.60 1.00 1 54.8 54.1 51.3 41.7 45.6 432 1 4 3 1 1 0.06 0.03 1.00 9 44.2 54.1 51.4 39.0 45.6 433 2 4 2 1 3 -0.64 0.03 1.00 4 41.6 46.3 38.4 36.3 45.6 434 2 4 3 1 2 0.48 0.03 1.00 4 65.4 61.9 60.4 63.4 60.5 435 2 4 1 1 1 -0.19 -0.85 0.00 15 54.8 60.6 55.5 58.0 45.6 436 1 4 3 1 2 0.20 -2.54 0.00 16 54.8 61.9 63.6 63.4 65.5 437 2 4 2 1 3 -0.60 -0.26 0.67 4 49.5 41.1 38.4 52.6 35.6 438 1 4 2 1 2 -1.33 0.03 1.00 16 36.3 28.1 53.7 33.6 25.7 439 2 4 1 1 2 0.68 0.88 1.00 9 46.9 61.9 53.0 52.6 60.5 440 2 4 3 1 2 1.36 0.94 1.00 4 52.1 48.9 51.3 41.7 45.6 441 2 4 2 1 1 0.89 0.65 0.67 9 68.0 48.9 56.3 63.4 55.6 442 2 4 2 1 3 -0.14 0.13 0.67 9 41.6 56.7 51.3 47.1 60.5 443 1 4 2 1 3 -0.43 -1.22 0.00 5 41.6 46.9 55.5 44.4 40.6 444 1 4 2 1 2 0.04 0.03 0.67 9 65.4 51.5 61.2 68.8 60.5 445 1 4 1 1 1 0.44 -0.16 0.67 9 62.7 48.9 48.8 66.1 45.6 446 1 4 1 1 2 0.65 -0.30 1.00 10 60.1 61.9 74.6 58.0 60.5 447 1 4 2 1 3 -1.33 0.65 0.00 2 40.5 38.5 47.1 44.4 30.6 448 1 4 2 1 3 0.68 1.19 0.33 14 41.6 41.1 57.2 71.5 30.6 449 1 4 2 1 3 0.10 0.03 0.00 3 62.7 48.9 42.1 66.1 70.5 450 2 4 3 1 2 0.50 0.03 0.67 17 68.0 59.3 58.8 66.1 65.5 451 2 4 2 1 3 0.02 -0.32 1.00 4 46.9 46.3 52.0 47.7 45.6 452 2 4 2 1 1 0.06 0.03 1.00 17 41.6 46.3 46.2 39.0 45.6 453 2 4 2 1 1 -0.83 0.65 0.33 4 44.2 59.3 55.5 36.3 55.6 454 2 4 3 1 2 0.91 0.15 0.00 9 54.8 56.7 55.5 58.0 60.5 455 2 4 1 1 1 0.06 0.59 1.00 1 52.1 51.5 45.4 39.0 50.6 456 1 4 2 1 3 -1.57 0.65 0.33 6 52.1 56.7 45.2 55.3 50.6 457 1 4 2 1 3 -1.33 0.03 0.33 5 38.9 30.7 40.4 39.0 50.6 458 1 4 2 1 2 0.46 0.65 1.00 16 49.5 48.9 60.5 55.3 55.6 459 2 4 2 1 1 0.03 0.03 1.00 11 41.6 43.0 44.5 39.0 50.6 460 1 4 3 1 2 -0.93 -0.22 0.67 10 65.4 64.5 65.5 60.7 65.5 461 1 4 1 1 2 0.24 -0.43 0.33 16 70.7 43.7 58.8 66.1 60.5 462 1 4 2 1 1 0.33 0.28 0.67 9 52.1 46.3 48.1 55.3 40.6 463 1 4 2 1 1 -0.11 0.25 1.00 9 44.2 44.3 45.6 39.0 50.6 464 2 4 2 1 2 0.94 -0.30 1.00 10 60.1 67.1 52.4 55.3 50.6 465 1 4 3 1 1 0.75 -0.49 0.33 16 41.6 41.1 41.2 44.4 35.6 466 2 4 2 1 2 0.30 -0.60 1.00 9 49.5 59.3 53.1 60.7 60.5 467 1 4 2 1 1 1.13 0.03 1.00 6 38.9 54.1 47.5 49.8 40.6 468 1 4 3 1 1 -0.16 0.03 0.33 11 44.2 59.3 49.5 55.3 45.6 469 1 4 2 1 1 0.46 0.65 1.00 16 68.0 59.3 57.9 74.2 65.5 470 2 4 3 1 2 1.36 0.94 1.00 10 57.4 64.5 49.9 55.3 55.6 471 2 4 3 1 2 0.46 0.32 0.67 16 49.5 54.1 47.8 60.7 50.6 472 2 4 2 1 1 -0.40 -0.26 0.33 15 52.1 54.1 55.3 52.6 50.6 473 2 4 2 1 3 0.06 0.28 1.00 13 60.1 59.3 56.2 60.7 60.5 474 2 4 1 1 3 -0.61 0.36 0.00 4 44.2 48.9 46.1 47.1 55.6 475 1 4 3 1 2 0.71 0.03 1.00 10 54.8 61.2 53.7 48.8 60.5 476 2 4 3 1 1 0.23 -0.26 0.00 9 52.1 61.9 56.3 58.0 55.6 477 2 4 2 1 2 0.46 -0.28 1.00 9 57.4 56.7 56.4 63.4 60.5 478 2 4 2 1 2 0.67 -0.59 0.67 9 54.8 64.5 55.3 49.8 60.5 479 2 4 2 1 2 0.46 0.03 1.00 10 60.1 64.5 56.3 63.4 60.5 480 1 4 1 1 3 -0.88 -0.29 0.33 8 44.2 54.1 41.0 49.8 45.6 481 2 4 1 1 2 0.89 0.34 0.33 10 57.4 61.9 72.2 60.7 60.5 482 2 4 2 1 3 0.06 -1.07 0.00 9 38.9 48.9 45.2 44.4 50.6 483 2 4 1 1 3 0.06 -0.60 0.67 9 62.7 56.0 57.3 60.7 56.8 484 2 4 2 1 3 -0.22 -2.27 0.00 1 57.4 51.5 40.3 60.7 55.6 485 2 4 2 1 2 0.71 0.34 0.00 2 62.7 64.5 57.4 60.7 70.5 486 2 4 2 1 2 0.53 0.81 0.67 9 54.8 59.3 61.4 47.1 55.6 487 2 4 1 1 2 -0.47 -0.47 0.00 1 52.1 57.3 63.7 58.0 55.6 488 1 4 3 1 3 0.43 -0.05 0.67 12 54.8 45.6 51.8 60.7 45.6 489 2 4 2 1 2 -0.44 -1.13 1.00 9 54.8 61.9 69.6 60.7 55.6 490 1 4 2 1 3 -1.28 0.34 0.33 9 31.0 35.9 46.1 39.0 45.6 491 2 4 2 1 1 -0.14 0.56 1.00 4 54.8 61.9 54.6 47.1 45.6 492 2 4 1 1 2 -0.14 -1.05 1.00 9 70.7 65.1 66.4 63.4 60.5 493 1 4 1 1 3 -0.86 0.28 1.00 3 36.3 48.9 54.4 60.7 35.6 494 2 4 1 1 3 0.30 0.03 0.67 1 36.8 59.3 40.7 49.8 35.6 495 1 4 2 1 1 0.68 0.03 1.00 10 65.4 64.5 51.3 66.1 65.5 496 2 4 2 1 2 0.23 0.03 1.00 9 52.1 59.3 52.9 60.7 50.6 497 2 4 2 1 3 0.00 -0.18 0.00 9 62.7 59.3 55.5 45.5 46.8 498 1 4 2 1 3 0.26 0.03 1.00 8 62.7 59.3 53.2 63.4 55.6 499 1 4 2 1 1 -1.58 -0.26 0.33 7 62.7 48.9 35.3 66.1 40.6 500 2 4 1 1 3 0.06 0.28 1.00 3 65.9 63.2 50.3 60.1 55.6 501 1 4 2 1 2 0.91 0.59 1.00 10 65.4 67.1 67.1 66.1 61.8 502 1 4 2 1 2 0.28 0.32 1.00 10 52.1 54.1 56.5 55.3 50.6 503 1 4 3 1 2 -0.19 -0.73 0.67 10 73.3 60.6 64.7 66.1 65.5 504 1 4 2 1 2 1.36 0.63 1.00 10 62.7 61.9 68.6 68.8 55.6 505 1 4 1 1 1 0.71 1.19 0.33 3 54.8 59.3 62.5 68.8 45.6 506 2 4 2 1 2 0.25 -0.47 0.33 10 76.0 67.1 72.2 66.1 65.5 507 1 4 2 1 3 -0.40 -0.76 0.67 7 52.1 35.9 50.3 47.1 50.6 508 1 4 3 1 2 0.68 -0.47 0.33 7 46.9 61.9 60.5 68.8 65.5 509 1 4 2 1 3 -0.41 0.34 0.67 16 52.1 38.5 57.9 52.6 50.6 510 1 4 2 1 3 -0.86 1.19 0.33 8 33.6 28.1 31.8 39.6 35.6 511 1 4 2 1 2 -0.44 1.19 1.00 2 60.1 54.1 59.6 55.3 65.5 512 1 4 1 1 3 -1.10 0.03 1.00 8 38.9 38.5 42.8 41.7 40.6 513 1 4 3 1 2 0.06 0.32 0.67 12 49.5 43.7 55.5 68.8 65.5 514 1 4 2 1 2 0.06 0.03 0.00 11 46.9 51.5 57.2 52.6 60.5 515 1 4 2 1 2 0.68 -0.26 1.00 10 62.7 61.9 56.2 47.1 45.6 516 1 4 2 1 2 0.66 0.03 0.33 6 65.4 59.3 67.1 66.1 65.5 517 2 4 1 1 1 0.26 0.87 1.00 3 49.5 61.9 61.4 63.4 50.6 518 1 4 2 1 1 0.06 0.03 1.00 2 54.8 41.1 42.8 47.1 40.6 519 1 4 2 2 2 0.32 1.19 1.00 9 57.4 61.9 64.0 55.8 55.6 520 1 4 2 2 2 0.06 0.03 0.33 9 68.0 59.3 71.3 66.1 55.6 521 1 4 3 2 2 0.46 0.34 1.00 16 52.1 46.3 50.4 52.6 50.6 522 1 4 2 2 2 0.00 0.65 1.00 7 52.1 61.9 62.1 58.0 60.5 523 2 4 3 2 1 0.68 0.32 1.00 15 36.3 56.7 41.9 49.8 40.6 524 2 4 2 2 2 -0.66 -1.07 0.67 1 49.5 61.9 60.4 47.1 50.6 525 2 4 3 2 2 0.48 -0.82 0.33 9 62.7 64.5 58.2 61.8 61.8 526 1 4 2 2 2 0.71 0.03 0.33 11 46.9 46.9 41.0 41.7 50.6 527 2 4 2 2 2 1.36 -1.18 1.00 9 62.7 61.9 63.1 60.7 60.5 528 2 4 3 2 2 -0.19 1.19 1.00 9 60.1 61.9 54.6 60.7 65.5 529 1 4 2 2 2 0.55 0.63 1.00 16 54.8 59.3 62.3 58.0 50.6 530 2 4 3 2 2 0.26 0.94 1.00 10 60.1 59.3 51.8 53.1 55.6 531 2 4 3 2 2 0.68 0.06 0.67 1 49.5 51.5 41.2 41.7 45.6 532 1 4 2 2 3 0.00 0.34 0.33 2 46.9 56.7 57.2 58.0 45.6 533 2 4 1 2 2 -0.17 0.03 0.67 15 44.2 59.3 45.7 55.3 50.6 534 1 4 2 2 2 0.06 0.56 1.00 10 49.5 56.7 46.2 55.3 50.6 535 2 4 2 2 2 0.52 0.34 1.00 10 46.9 64.5 59.6 49.8 55.6 536 2 4 3 2 2 0.91 -1.67 0.33 9 57.4 54.1 54.8 60.7 60.5 537 1 4 2 2 1 -0.40 0.65 0.00 10 52.1 46.3 51.3 58.0 55.6 538 2 4 3 2 2 0.46 -0.22 1.00 10 70.7 64.5 68.7 58.0 70.5 539 2 4 2 2 2 -0.40 -1.34 1.00 8 57.4 61.9 55.5 49.8 55.6 540 1 4 2 2 2 0.71 1.19 0.67 2 54.8 48.9 52.4 58.0 60.5 541 1 4 2 2 2 -0.14 0.34 0.33 9 49.5 46.3 44.8 58.0 60.5 542 2 4 2 2 2 0.52 0.03 1.00 14 46.9 51.5 53.7 41.7 40.6 543 2 4 3 2 2 -0.47 0.28 1.00 10 65.4 61.9 47.1 60.7 60.5 544 2 4 2 2 2 0.96 0.03 0.67 11 44.2 51.5 43.2 44.4 50.6 545 2 4 2 2 2 0.91 -0.47 1.00 9 46.9 64.5 55.8 61.8 60.5 546 2 4 3 2 2 0.02 -0.14 0.00 9 60.1 64.5 67.0 58.5 66.8 547 1 4 2 2 2 0.47 0.03 0.67 13 62.7 59.3 48.8 55.3 70.5 548 2 4 2 2 2 0.47 0.34 0.67 10 33.6 43.0 41.0 49.8 35.6 549 2 4 3 2 2 0.46 0.03 1.00 9 41.6 64.5 47.1 53.1 60.5 550 2 4 2 2 3 1.13 0.88 0.33 1 49.5 51.5 52.7 55.3 55.6 551 1 4 2 2 2 0.96 0.63 1.00 12 65.4 64.5 70.3 66.1 65.5 552 2 4 2 2 2 0.50 -0.22 1.00 9 54.8 61.9 61.3 60.7 50.6 553 1 4 2 2 2 0.71 0.28 0.67 16 62.7 56.7 54.7 58.0 40.6 554 2 4 1 2 2 -0.14 -0.60 0.33 15 44.2 48.9 54.3 49.8 45.6 555 1 4 3 2 2 0.06 1.19 1.00 11 65.4 64.5 67.9 63.4 55.6 556 2 4 2 2 2 0.87 0.87 0.67 10 57.4 61.9 62.8 55.3 40.6 557 2 4 2 2 2 1.16 1.19 1.00 16 60.1 59.3 55.3 58.0 45.6 558 2 4 2 2 3 0.46 0.56 0.67 9 44.2 54.1 52.2 44.4 45.6 559 2 4 2 2 1 0.27 0.03 0.67 4 57.4 41.1 57.0 55.3 51.8 560 2 4 1 2 2 -0.19 -1.73 1.00 1 44.7 47.6 50.3 39.0 50.6 561 2 4 3 2 2 0.68 0.03 1.00 6 44.2 48.9 46.2 47.7 45.6 562 2 4 3 2 2 0.23 0.03 1.00 1 62.7 61.9 55.5 55.3 60.5 563 2 4 3 2 2 0.91 0.03 1.00 9 57.4 59.3 55.5 55.3 50.6 564 1 4 3 2 2 0.32 -0.88 0.33 9 57.4 54.1 42.2 44.4 55.6 565 1 4 2 2 2 0.71 0.03 0.33 5 46.9 59.3 63.0 52.6 45.6 566 1 4 3 2 2 1.11 0.34 1.00 10 65.4 64.5 64.1 66.1 70.5 567 2 4 2 2 2 0.26 -0.57 1.00 10 41.6 58.6 54.6 55.3 60.5 568 2 4 2 2 2 0.00 0.34 0.33 9 46.9 59.3 53.7 58.0 45.6 569 1 4 1 2 2 0.20 -0.47 0.67 2 56.4 49.5 51.9 47.1 48.1 570 2 4 2 2 2 -0.84 -2.60 0.33 15 44.2 64.5 53.0 55.3 60.5 571 2 4 3 2 2 1.16 0.09 1.00 9 46.9 52.1 42.5 47.7 60.5 572 1 4 2 2 3 0.23 0.15 0.33 2 49.5 51.5 48.0 52.6 50.6 573 2 4 1 2 1 -0.23 0.69 0.67 15 46.9 61.9 48.0 39.0 50.6 574 1 4 3 2 2 0.03 0.56 0.67 2 65.4 67.1 63.0 66.1 70.5 575 2 4 2 2 1 0.68 0.03 0.00 10 44.2 35.9 43.6 47.1 40.6 576 2 4 2 2 3 0.43 -1.09 0.67 14 63.3 64.5 52.1 58.0 65.5 577 2 4 2 2 2 0.47 0.34 0.67 16 44.2 48.9 48.0 39.0 50.6 578 2 4 2 2 3 -0.59 0.28 1.00 14 46.9 51.5 48.5 49.8 50.6 579 2 4 1 2 3 -0.38 0.03 1.00 9 45.3 38.5 43.5 40.1 45.6 580 2 4 3 2 2 0.73 0.03 1.00 9 62.7 62.5 68.9 60.7 60.5 581 1 4 2 2 3 -0.66 -0.47 1.00 9 52.7 41.1 40.2 58.0 60.5 582 1 4 3 2 2 0.10 -0.16 0.33 9 49.5 59.3 51.0 47.1 45.6 583 1 4 2 2 1 -1.33 0.03 0.67 16 57.4 56.7 59.8 58.0 55.6 584 1 4 3 2 2 0.22 1.19 1.00 11 57.4 54.1 48.8 55.3 40.6 585 1 4 3 2 2 0.10 0.56 0.67 9 54.8 61.9 60.4 55.3 65.5 586 1 4 3 2 2 0.71 -0.06 0.67 9 44.2 37.8 49.2 39.0 45.6 587 1 4 3 2 2 0.68 -0.53 1.00 11 48.4 48.2 60.5 60.1 49.3 588 1 4 2 2 1 -0.93 -0.80 0.33 9 52.1 51.5 49.5 55.3 45.6 589 1 4 3 2 2 -0.60 0.34 0.00 2 49.5 41.7 50.3 36.3 60.5 590 1 4 2 2 2 0.06 0.03 0.67 16 46.9 56.7 56.1 60.7 55.6 591 2 4 3 2 2 1.14 0.63 1.00 15 55.3 53.4 44.5 53.6 48.1 592 2 4 3 2 2 0.06 0.03 0.67 14 46.9 60.6 51.3 63.4 30.6 593 1 4 2 2 2 -1.09 -0.26 0.33 11 44.2 41.1 45.1 47.1 45.6 594 1 4 3 2 2 0.02 0.28 0.67 9 44.2 51.5 52.9 52.6 45.6 595 1 4 3 2 2 0.00 0.03 0.33 10 52.1 59.3 50.4 60.7 60.5 596 1 4 3 2 2 0.27 -1.05 0.33 9 60.1 54.1 56.3 55.3 55.6 597 1 4 2 2 2 0.55 0.34 1.00 16 62.7 61.9 72.9 63.4 55.6 598 1 4 2 2 2 -0.61 0.34 0.67 16 68.0 54.1 74.6 66.1 65.5 599 2 4 2 2 2 0.23 0.94 1.00 10 57.4 67.1 57.9 60.7 55.6 600 1 4 2 2 2 0.02 0.03 0.33 9 62.7 54.1 64.7 58.0 50.6 ; run; run; /* Model I: */ proc genmod data=hsb; model academic/ncases = math / link=logit dist=binomial obstats type3; output out=preds pred=fitted; run; proc sort data=preds; by fitted; goptions reset=(axis, legend, pattern, symbol, title, footnote) norotate hpos=0 vpos=0 htext=2.5 ftext=swiss ctext= target= gaccess= gsfmode= ; goptions device=WIN ctext=blue graphrc interpol=join; axis1 color=black width=2.0 label=('Math Achievement Score') ; axis2 color=black width=2.0 label=(angle=90 'Probability Attend Academic') order=0 to 1 by .1 ; symbol1 value=none color=blue i=join; proc gplot data=preds; plot fitted*math / haxis=axis1 vaxis=axis2 frame ; run; /* Model II: */ proc genmod data=hsb; class SES; model academic/ncases = math ses/ link=logit dist=binomial obstats type3; output out=preds pred=fitted; run; proc sort data=preds; by fitted; goptions reset=(axis, legend, pattern, symbol, title, footnote) norotate hpos=0 vpos=0 htext=2.5 ftext=swiss ctext= target= gaccess= gsfmode= ; goptions device=WIN ctext=blue graphrc interpol=join; axis1 color=black width=2.0 label=('Math Achievement Score') ; axis2 color=black width=2.0 label=(angle=90 'Probability Attend Academic') order=0 to 1 by .1 ; symbol1 value=none color=red i=join; symbol2 value=none color=blue i=join; symbol3 value=none color=black i=join; legend1 position=(top center outside) frame cshadow=pink label=('SES:') value=('Low' 'Middle' 'High') ; proc gplot data=preds; plot fitted*math=ses / haxis=axis1 vaxis=axis2 frame legend=legend1 ; run; /* Model III: */ proc genmod data=hsb; model academic/ncases = math ses/ link=logit dist=binomial obstats type3; output out=preds pred=fitted; run; proc sort data=preds; by fitted; goptions reset=(axis, legend, pattern, symbol, title, footnote) norotate hpos=0 vpos=0 htext=2.5 ftext=swiss ctext= target= gaccess= gsfmode= ; goptions device=WIN ctext=blue graphrc interpol=join; axis1 color=black width=2.0 label=('Math Achievement Score') ; axis2 color=black width=2.0 label=(angle=90 'Probability Attend Academic') order=0 to 1 by .1 ; symbol1 value=none color=red i=join; symbol2 value=none color=blue i=join; symbol3 value=none color=black i=join; legend1 position=(top center outside) frame cshadow=pink label=('SES:') value=('Low' 'Middle' 'High') ; proc gplot data=preds; plot fitted*math=ses / haxis=axis1 vaxis=axis2 frame legend=legend1 ; run; /* Model IV: */ proc genmod data=hsb; class sctyp; model academic/ncases = math ses sctyp / link=logit dist=binomial obstats type3; output out=preds pred=fitted; run; data pred2; set preds; if ses=1 and sctyp=1 then grp=1; else if ses=2 and sctyp=1 then grp=2; else if ses=3 and sctyp=1 then grp=3; else if ses=1 and sctyp=2 then grp=4; else if ses=2 and sctyp=2 then grp=5; else if ses=3 and sctyp=3 then grp=6; proc sort data=pred2; by math grp; goptions reset=(axis, legend, pattern, symbol, title, footnote) norotate hpos=0 vpos=0 htext=2.5 ftext=swiss ctext= target= gaccess= gsfmode= ; goptions device=WIN ctext=blue graphrc interpol=join; axis1 color=black width=2.0 label=('Math Achievement Score') ; axis2 color=black width=2.0 label=(angle=90 'Probability Attend Academic') order=0 to 1 by .1 ; symbol1 value=none color=red i=join; symbol2 value=none color=blue i=join; symbol3 value=none color=black i=join; symbol4 value=none color=red i=join; symbol5 value=none color=blue i=join; symbol6 value=none color=black i=join; legend1 position=(top center outside) frame cshadow=pink label=none value=( 'private Low' 'private Middle' 'private High' 'public Low' 'public Middle' 'public High') ; proc gplot data=pred2; plot fitted*math=grp / haxis=axis1 vaxis=axis2 frame legend=legend1 ; note move=(25,30) 'Private'; note move=(125,41) 'Public'; run; /* Model V: */ proc genmod data=hsb; class sctyp; model academic/ncases = math ses sctyp ses*sctyp / link=logit dist=binomial obstats type3; output out=preds pred=fitted stdreschi=adjusted; run; proc sort data=preds; by math grp; goptions reset=(axis, legend, pattern, symbol, title, footnote) norotate hpos=0 vpos=0 htext=2.5 ftext=swiss ctext= target= gaccess= gsfmode= ; goptions device=WIN ctext=blue graphrc interpol=join; axis1 color=black width=2.0 label=('Math Achievement Score') ; axis2 color=black width=2.0 label=(angle=90 'Probability Attend Academic') order=0 to 1 by .1 ; axis3 color=black width=2.0 label=(angle=90 'Adjusted Residuals') ; symbol1 value=none color=red i=join; symbol2 value=none color=blue i=join; symbol3 value=none color=black i=join; symbol4 value=none color=green i=join; symbol5 value=none color=orange i=join; symbol6 value=none color=brown i=join; symbol7 value=circle color=blue i=none; proc gplot data=preds; plot fitted*math=grp / haxis=axis1 vaxis=axis2 frame nolegend ; note move=(25,49) 'high private'; note move=(22,35) 'middle private'; note move=(28,26) 'high public'; note move=(65.5,30.5) 'middle'; note move=(59.5,29) 'public'; note move=(57,20) 'low public & private'; data Quants; do i=1 to 600 by 1; p = i/601; z=quantile('normal',p); output; end; run; proc sort data=preds; by adjusted; data QQplot; merge Quants preds; run; axis1 color=black width=2.0 label=('Standard Normal Quantiles') ; proc gplot data=qqplot; plot adjusted*z=7 / haxis=axis1 vaxis=axis3 frame nolegend ; run; /* Model VI: */ proc genmod data=hsb; class sctyp; model academic/ncases = math ses sctyp ses*sctyp rdg / link=logit dist=binomial obstats type3; output out=preds pred=fitted stdreschi=adjusted; run; proc logistic data=hsb; class sctyp; model academic = math ses sctyp ses*sctyp rdg / lackfit influence; run; proc logistic data=hsb; class sctyp; model academic = math ses sctyp ses*sctyp rdg math*rdg/ lackfit influence; run; /* Model VII: */ proc genmod data=hsb; class sctyp; model academic/ncases = math ses sctyp ses*sctyp locus mot / link=logit dist=binomial obstats type3; output out=preds pred=fitted stdreschi=adjusted; run; /* Model Last one */ proc logistic data=hsb; class sctyp; model academic/ncases = math ses sctyp ses*sctyp math*ses math*sctyp math*sctyp*ses / lackfit; run; /* Model Last one */ proc logistic data=hsb; class sctyp; model academic/ncases = math ses sctyp ses*sctyp / lackfit; run; /* Model Last one */ proc logistic data=hsb; class sctyp; model academic/ncases = math ses sctyp ses*sctyp rdg sci wrtg civ / lackfit; run; /* Model V, Last one */ proc logistic data=hsb; class sctyp; model academic/ncases = math ses sctyp ses*sctyp sci / lackfit ctable outroc=roc pprob=0.53; output out=pred pred=fitted; run; data preds; set pred; if ses=1 and sctyp=1 then interaction="Low-Private"; else if ses=1 and sctyp=2 then interaction="Low-Public "; else if ses=2 and sctyp=1 then interaction="Mid-Private"; else if ses=2 and sctyp=2 then interaction="Mid-Public "; else if ses=3 and sctyp=1 then interaction="Hig-Private"; else if ses=3 and sctyp=2 then interaction="Hig-Public "; run; symbol1 c=blue i=join; symbol2 c=red i=join; symbol3 c=pink i=join; symbol4 c=cyan i=join; symbol5 c=green i=join; symbol6 c=grey i=join; axis1 label=(r=0 a=90) ; axis2 order=30 to 80 by 10; proc sort data=preds; by math; title "Easy Way to make plot"; proc gplot data= preds; plot fitted*math=interaction / vaxis=axis1 haxis=axis2; run; /* Model V to get ROC */ proc logistic data=hsb; class sctyp; model academic/ncases = math ses sctyp ses*sctyp / lackfit ctable outroc=roc ; run; goptions reset=(axis, legend, pattern, symbol, title, footnote) norotate hpos=0 vpos=0 htext=2.5 ftext=swiss ctext= target= gaccess= gsfmode= ; goptions device=WIN ctext=blue graphrc interpol=join; axis1 color=black width=2.0 label=('1 - Specificity (False Positives)') ; axis2 color=black width=2.0 label=(angle=90 'Sensitivity (True Positives)') order=0 to 1 by .1 ; symbol1 value=none color=red i=join; proc gplot data=roc; plot _SENSIT_ * _1MSPEC_ / haxis=axis1 vaxis=axis2 frame nolegend ; title 'c=.789'; run; /* Model VI */ proc logistic data=hsb; class sctyp ; model academic/ncases = math ses sctyp ses*sctyp math*ses math*sctyp math*ses*sctyp ; run; proc means data=hsb; var math ses; output out=mymeans mean=mmath mses ; run; data tmp; set mymeans; do id=1 to 600; output; end; run; data chsb; merge hsb tmp; cmath=math-mmath; cses = ses-mses; run; /* Model VI */ proc logistic data=chsb; class sctyp ; model academic/ncases = cmath cses sctyp cses*sctyp cmath*cses cmath*sctyp cmath*cses*sctyp ; run; /* Model VII */ proc logistic data=hsb; class sctyp ; model academic/ncases = math ses sctyp ses*sctyp rdg sci wrtg civ ; run; proc princomp data=hsb; var rdg wrtg math sci civ; run; /* Model VII with centered */ title "Model VII with centered "; proc logistic data=chsb; class sctyp ; model academic/ncases = cmath cses sctyp cses*sctyp crdg csci cwrtg cciv ; run;