dm "output;clear;log;clear"; GOPTIONS reset=all goutmode=replace; OPTIONS PAGENO=1; OPTIONS LINESIZE=75 PAGESIZE=66 NODATE; ODS PDF File="C:\teaching\math593\ex7_1.pdf"; ODS Listing Close; GOptions NoPrompt Vsize=6 Hsize=6 Horigin=1.2 Vorigin=2.5 FText=SwissX FTitle=SwissX HText=1 HTitle=1; TITLE 'Ex. 7.1 - Bone Measurements on Turkeys'; DATA TURKEY; INPUT ID $ 1-5 S 6 T 7 SIDE $ 8 WGT 10-12 HUM 14-16 RAD 18-20 ULN 22-24 FEMUR 26-28 TIB 30-32 TIN 34-36 CAR 38-40 D3P 42-44 STL 46-48 STB 50-52 COR 54-56 PEL 58-60 MAX 62-64 MIN 66-68 SCA 70-72; IF S=1 THEN SEX='MALE '; ELSE SEX='FEMALE'; IF T=1 THEN TYPE='WILD '; ELSE TYPE='DOMESTIC'; CARDS; K766 11L 142 107 N399 11L 153 138 153 139 246 162 810 307 196 74 145 104 NEX1 11L 224 72 NEX2 11L 220 74 NEX3 11L 228 78 NEB4 11L 217 76 NEB5 11L 230 81 NEB6 11L 210 74 NEB7 11L 215 71 NE16 11L 204 69 B396 11 140 107 B457 11 143 151 108 132 B710 11 153 140 147 142 151 817 305 102 128 B790 11 156 137 151 146 155 814 305 111 137 B791 11 132 148 138 145 775 106 128 B795 11 151 134 151 144 789 292 116 126 B819 11 158 135 151 146 152 790 289 111 125 B081 11 135 149 149 789 111 123 B085 11 148 129 146 139 147 767 287 106 123 B089 11 157 140 154 140 159 818 301 116 136 B090 11 153 138 153 141 151 822 312 115 133 B091 11 156 138 156 145 150 835 310 118 133 B093 11 151 133 148 139 152 793 290 105 B097 11 153 135 150 144 158 772 276 102 123 B099 11 152 140 151 144 158 792 303 111 122 B101 11 151 814 304 B102 11 153 145 156 115 136 B103 11 149 133 149 139 151 107 132 B105 11 145 160 116 134 B106 11 147 130 144 136 145 765 289 108 131 B111 11 154 138 155 142 153 827 315 111 128 B114 11 154 138 154 138 155 802 287 111 132 B116 11 161 131 150 140 151 816 301 112 134 B117 11 142 159 146 156 824 298 113 136 L641 11L 150 133 226 141 216 68 111 177 133 99 122 L690 11L 140 126 137 123 218 141 770 290 150 64 95 145 128 87 118 L783 11L 154 144 237 157 219 77 113 186 161 107 135 L902 11L 126 145 790 300 L905 11L 240 169 154 255 163 241 76 122 207 141 101 143 K263 21L 172 54 110 80 K694 21L 57 119 83 K767 21L 128 51 103 85 NEB1 21L 188 59 NEB2 21L 184 54 B426 21 123 108 122 112 122 647 241 84 106 B455 21 126 109 123 117 123 682 260 86 111 B724 21 126 112 126 117 128 660 243 B821 21 125 110 122 118 125 657 251 86 112 B845 21 122 107 120 112 641 86 106 B893 21 120 109 120 115 121 657 256 86 103 B086 21 121 108 120 113 120 635 86 105 B112 21 125 110 123 115 123 656 248 90 111 B122 21 123 107 122 113 124 652 255 82 109 L782 21L 130 127 119 191 126 176 57 92 130 87 L790 21L 116 103 114 106 175 119 141 106 77 L637 12L 200 153 137 140 820 300 104 166 130 L674 12L 200 142 131 140 128 199 131 800 250 133 84 95 158 173 125 122 L678 12L 215 148 130 145 133 216 141 800 290 145 84 98 163 180 134 124 L679 12L 199 148 131 144 133 211 135 820 300 140 86 103 170 178 137 123 L682 12L 190 145 131 140 129 215 136 790 300 134 86 100 159 182 129 123 L684 12L 200 153 135 146 138 232 146 840 310 149 93 107 166 189 140 125 L695 12L 149 132 147 131 218 139 810 300 131 83 102 160 178 130 126 L735 12L 129 118 127 280 81 97 L746 12L 215 143 129 141 130 200 137 800 300 138 84 97 163 186 138 121 L747 12L 233 149 132 147 137 219 136 790 300 164 83 105 173 178 148 127 L748 12L 204 145 130 145 132 215 136 760 280 145 85 97 162 178 129 123 L749 12L 203 148 135 145 135 212 145 840 310 100 182 129 124 L750 12L 220 149 130 147 140 225 147 770 300 155 88 104 166 184 141 126 L753 12L 141 127 136 126 197 131 830 310 97 168 119 120 L755 12L 164 143 130 139 129 207 144 790 97 162 112 118 L758 12L 239 149 131 146 135 219 137 810 300 138 82 102 169 172 132 127 L762 12L 145 131 143 136 210 140 790 290 137 82 101 162 182 133 125 L763 12L 167 143 129 138 127 204 144 790 93 167 126 116 L765 12L 188 148 137 146 133 215 144 820 300 139 91 99 165 171 129 125 L770 12L 190 144 129 138 130 201 134 790 300 97 177 127 118 L774 12L 213 143 128 141 130 203 137 800 300 138 81 98 166 183 143 123 L776A12L 133 218 138 169 186 138 L776B12L 138 225 141 168 194 146 L779 12L 220 149 136 147 135 208 139 820 300 141 92 98 162 182 143 126 L784 12L 220 150 135 148 134 204 143 820 300 142 82 103 161 184 140 129 NED1 12L 145 84 NED2 12L 145 84 NEG1 12L 207 84 NEG2 12L 89 NEG3 12L 86 NEG4 12L 86 NEG5 12L 93 NEG6 12L 88 L638 22L 120 122 108 117 109 176 115 670 250 111 70 82 131 146 114 102 L671 22L 104 122 106 117 110 177 112 660 250 116 65 81 131 141 111 103 L672 22L 120 122 109 120 113 178 114 670 26 125 86 136 147 102 107 L673 22L 135 121 105 118 112 176 113 650 240 129 68 83 134 140 106 109 L675 22L 118 101 114 110 168 110 630 250 121 80 132 157 119 107 L676 22L 119 102 115 111 167 106 640 250 137 67 84 135 102 108 L677 22L 123 118 104 115 107 167 111 640 240 120 62 80 129 148 105 103 L680 22L 119 106 117 109 170 105 660 260 121 68 83 129 147 107 103 L683 22L 114 121 103 116 110 172 106 650 250 130 64 81 131 144 104 104 L685 22L 119 107 118 109 178 114 650 250 119 68 80 125 147 110 99 L691 22L 116 120 106 119 113 176 112 660 260 131 62 84 130 138 100 106 L692 22L 122 106 117 111 171 108 670 260 120 85 136 101 103 L694 22L 121 106 118 111 176 111 660 253 128 69 85 133 141 104 107 L696 22L 140 126 112 180 117 680 250 132 68 85 134 146 111 112 L697 22L 117 121 107 117 109 177 114 670 260 120 68 82 131 145 111 104 L740 22L 118 119 104 116 109 174 111 650 250 123 66 84 130 146 97 103 L743 22L 109 118 104 115 106 167 105 660 260 120 66 80 129 138 102 98 L751 22L 121 126 112 122 114 180 119 690 260 120 68 83 138 147 115 106 L754 22L 109 101 107 102 153 106 240 71 137 101 92 L756 22L 125 123 108 120 114 177 118 680 250 126 66 82 151 108 L757 22L 113 122 107 120 112 179 114 670 260 126 65 80 134 134 104 103 L759 22L 139 119 107 118 111 178 116 670 260 121 71 82 136 142 105 107 L760 22L 126 109 122 117 183 119 690 260 132 68 86 137 142 109 111 L761 22L 111 124 109 119 113 172 115 680 260 108 63 84 133 135 104 106 L764 22L 133 120 105 117 109 173 111 670 260 130 64 82 136 152 110 108 L766 22L 121 108 119 109 178 115 660 250 115 64 80 134 146 109 105 L767 22L 129 114 126 118 188 120 700 270 123 70 85 161 116 113 L768 22L 108 124 108 121 112 181 114 660 260 123 65 81 132 130 101 103 L769 22L 116 122 107 119 113 174 112 660 250 118 69 81 135 146 105 106 L771 22L 111 114 101 113 104 158 106 630 240 113 66 77 123 139 102 101 L772 22L 131 128 110 125 120 183 118 680 270 147 69 88 140 158 111 109 L773 22L 108 120 107 120 114 175 114 660 260 120 63 83 136 140 106 107 L775 22L 127 110 125 120 186 120 700 124 75 85 144 111 112 L786 22L 120 108 117 108 165 110 680 113 66 81 130 137 104 103 L792 22L 122 123 110 120 112 180 114 117 72 86 136 135 100 107 NEG1 22L 183 61 NEG2 22L 193 60 NEG3 22L 186 62 DOR1 12L 212 132 173 130 DOR2 12L 223 140 171 138 DOR3 12L 227 141 182 126 DOR4 12L 225 145 193 142 DOR5 12L 207 139 174 131 DOR6 12L 214 145 181 140 DOR7 12L 211 144 185 133 DOR8 12L 214 147 175 125 DOR9 12L 207 139 170 130 DO10 12L 204 136 179 130 DO11 22L 181 113 146 111 DO12 22L 181 124 152 104 DO13 22L 186 121 150 109 DO14 22L 178 115 135 101 DO15 22L 179 110 138 102 DO16 22L 182 116 131 103 DO17 22L 180 115 148 114 DO18 22L 177 114 141 109 DO19 22L 178 116 149 111 DO20 22L 179 113 142 112 DO21 22L 180 110 141 106 DO22 22L 182 113 146 108 DO23 22L 178 111 140 102 DO24 22L 176 115 148 112 DO25 22L 174 115 133 103 DO26 22L 170 107 152 106 DO27 22L 173 108 152 108 DO28 22L 181 118 147 111 DO29 22L 183 115 129 102 DO30 22L 177 113 135 101 DO31 22L 182 117 148 109 DO32 22L 180 112 146 112 RUN; DATA TURKEY2; SET TURKEY; IF SEX='MALE'; KEEP ID TYPE HUM RAD ULN FEMUR TIN CAR D3P COR SCA; PROC PRINT DATA=TURKEY2; RUN; PROC DISCRIM CROSSVALIDATE MAHALANOBIS; CLASS TYPE; VAR HUM--SCA; RUN; title2 'Stepwise Seletion'; PROC STEPDISC DATA=TURKEY2 METHOD=Stepwise SLE=.40 SLS=.15; CLASS TYPE; VAR HUM--SCA; RUN; /*********** PLOTTING ************/ ** For plotting TIN and COR are used since they have the highest F in STEP Process **; *** Following is for Contours and Linear DISC Func ********; Proc Means Data=TURKEY2 NoPrint; Var TIN COR; Output Out=Turkey3 Min=MinT MinC Max=MaxT MaxC; Run; Data PlotTurkey; If _N_=1 Then Set Turkey3; Inc1=(MaxT-MinT)/50; Inc2=(MaxC-MinC)/50; Do TIN = (MinT-Inc1) To (MaxT+Inc1) By Inc1; Do COR = (MinC-Inc2) To (MaxC+Inc2) By Inc2; Output; Keep TIN COR; End; End; Stop; Run; Proc DISCRIM Data=TURKEY2 Testdata=PlotTurkey TestOut=PlotP TestOutD=PlotD; Class TYPE; Var TIN COR; Run; ** PLOT #1: Scatter plot TIN vs COR by TYPE ******; PROC GPLOT DATA=TURKEY2; plot TIN*COR = TYPE/ HAxis=Axis1 VAxis=Axis2; Axis1 Order=(80 To 130 By 10.0); Axis2 Order=(125 To 165 By 10.0);; Symbol1 V=circle H=0.7 I=None C=RED; Symbol2 V=Star H=0.7 I=None C=BLACK; run; quit; ** PLOT #2 #3: contour of group means ******; Title2 'Plot of Estimated Densities'; Proc GContour Data=PlotD; Title3 "WILD"; Plot TIN*COR=wild / HAxis=Axis1 VAxis=Axis2; Axis1 Order=(80 To 130 By 10.0); Axis2 Order=(125 To 165 By 10.0); Run; quit; Proc GContour Data=PlotD; Title3 "DOMESTIC"; Plot TIN*COR=domestic/ HAxis=Axis1 VAxis=Axis2; Axis1 Order=(80 To 130 By 10.0); Axis2 Order=(125 To 165 By 10.0);; Run; quit; ** PLOT #4: Classfication Results ******; Title2 'Plot of Classification Results'; Proc GPlot Data=PlotP; Plot TIN*COR=_Into_; Symbol1 V=Dot H=0.7 I=None C=RED; Symbol2 V=square H=0.7 I=None C=BLUE; Run; quit; title2 'Different Cost function C(D|W)/C(W|D) = 4'; /******************************************************************************** To have a different costs for DA we may use adjusted priors for Priors statement. For example, C(D|W)/C(W|D) = 4 then use priors P(W)=C(D|W)/(C(D|W)+C(W|D)) = 0.8 and P(D)= 0.2 *********************************************************************************/ PROC DISCRIM DATA=TURKEY2 CROSSVALIDATE MAHALANOBIS; CLASS TYPE; VAR HUM--SCA; PRIORS 'WILD'=0.8 'DOMESTIC'=0.2; RUN; ODS Listing; ODS PDF Close;