options ls = 75 ps=55; data sasuser.drink; input Number /* Number in study */ Sex /* M = 1, F = 0 */ drinker /* regular drinker, >= 14 units/wk. Y = 1, N = 0 */ Age /* years */ smoker /* regular smoker Y = 1, N = 0 */ pure /* was drink "pure" (orange juice or lager)? Y = 1, N = 0 */ dr /* which drink? Lager = 1, orange juice = 0 */ concent /* a coded concentration: 0,1,2,3,4, according to the order of concentration in the paper. 0 = no added alcohol to either drink, but 4 is a different concentration for lager than for orange juice. Within drink, the numbers reflect an ordinal scale */ adj /* actual relative concentration, specific to each drink (e.g the %-ages for OJ have been re-scaled to 0,1,2,4,6, and for lager to 0,4,7,10,14, to reflect the relative magnitudes of the added alcohol. But the scales are specific to each drink - you'd have to multiply the OJ scale by some factor to express it on the same scale as the lager. Or just convert all back to % of drink volume (as explained in the paper) */ event /* which event? Medsoc ball = 0, Public Health Awayday = 1 */ spiked /* the complement of the variable "pure" (to get the probability that the drink has been spiked */ filter /* a filter for when I looked only at the OJ data */ ojadj ; /* the index used for the paper, added alcohol expressed as % weight by volume (whatever that means). */ lines; 104 0 1 19 0 1 1 0 0 0 0 0 0.00 100 1 1 19 1 1 1 0 0 0 0 0 0.00 95 1 1 19 0 0 1 0 0 0 1 0 0.00 92 0 0 20 0 1 1 0 0 0 0 0 0.00 76 0 1 19 0 0 1 0 0 0 1 0 0.00 50 0 1 19 0 1 1 0 0 0 0 0 0.00 34 0 1 19 0 1 1 0 0 0 0 0 0.00 12 1 1 19 0 1 1 0 0 0 0 0 0.00 9 1 1 22 0 1 1 0 0 0 0 0 0.00 89 1 0 20 0 1 1 1 4 0 0 0 3.20 88 1 1 20 0 1 1 1 4 0 0 0 3.20 63 0 1 21 0 1 1 1 4 0 0 0 3.20 46 1 1 20 0 1 1 1 4 0 0 0 3.20 42 0 1 19 0 1 1 1 4 0 0 0 3.20 30 1 1 19 0 1 1 1 4 0 0 0 3.20 27 1 1 19 0 1 1 1 4 0 0 0 3.20 11 1 1 20 0 1 1 1 4 0 0 0 3.20 10 1 0 29 0 0 1 1 4 0 1 0 3.20 4 1 0 40 0 1 1 1 4 0 0 0 3.20 110 1 1 19 0 1 1 2 7 0 0 0 5.50 106 1 1 19 0 0 1 2 7 0 1 0 5.50 97 0 1 20 0 0 1 2 7 0 1 0 5.50 79 1 1 21 0 0 1 2 7 0 1 0 5.50 77 0 1 20 0 1 1 2 7 0 0 0 5.50 57 1 1 20 0 0 1 2 7 0 1 0 5.50 56 0 1 21 0 0 1 2 7 0 1 0 5.50 43 0 1 22 0 1 1 2 7 0 0 0 5.50 8 1 1 31 0 1 1 2 7 0 0 0 5.50 83 0 1 19 1 1 1 3 10 0 0 0 7.90 53 0 1 18 0 0 1 3 10 0 1 0 7.90 52 1 1 19 0 0 1 3 10 0 1 0 7.90 45 0 1 20 0 0 1 3 10 0 1 0 7.90 68 1 1 20 0 0 1 3 10 0 1 0 7.90 38 1 1 22 0 0 1 3 10 0 1 0 7.90 35 1 1 19 0 0 1 3 10 0 1 0 7.90 25 1 1 19 0 0 1 3 10 0 1 0 7.90 23 1 1 21 0 0 1 3 10 0 1 0 7.90 14 1 1 19 0 0 1 3 10 0 1 0 7.90 113 0 1 19 0 0 1 4 14 0 1 0 11.10 105 0 0 19 0 0 1 4 14 0 1 0 11.10 84 1 1 20 0 0 1 4 14 0 1 0 11.10 69 1 0 20 0 0 1 4 14 0 1 0 11.10 48 1 1 19 0 0 1 4 14 0 1 0 11.10 31 1 1 19 0 0 1 4 14 0 1 0 11.10 29 1 1 20 1 0 1 4 14 0 1 0 11.10 26 1 1 21 0 0 1 4 14 0 1 0 11.10 2 1 1 20 0 0 1 4 14 0 1 0 11.10 81 1 1 20 0 1 0 0 0 0 0 1 0.00 78 1 1 28 0 1 0 0 0 0 0 1 0.00 67 0 1 20 0 1 0 0 0 0 0 1 0.00 64 1 1 20 0 1 0 0 0 0 0 1 0.00 62 1 1 22 0 1 0 0 0 0 0 1 0.00 58 0 1 20 0 1 0 0 0 0 0 1 0.00 44 0 1 19 0 1 0 0 0 0 0 1 0.00 40 0 1 20 0 1 0 0 0 0 0 1 0.00 36 1 1 23 0 1 0 0 0 0 0 1 0.00 32 0 1 18 0 1 0 0 0 0 0 1 0.00 24 1 1 19 0 1 0 0 0 0 0 1 0.00 18 1 1 20 0 1 0 0 0 0 0 1 0.00 17 1 1 18 0 0 0 0 0 0 1 1 0.00 13 0 1 20 0 1 0 0 0 0 0 1 0.00 1 0 1 20 0 1 0 0 0 0 0 1 0.00 91 0 0 20 0 1 0 1 1 0 0 1 2.00 90 0 1 20 0 0 0 1 1 0 1 1 2.00 82 1 1 21 0 1 0 1 1 0 0 1 2.00 74 1 1 20 0 1 0 1 1 0 0 1 2.00 70 0 1 21 0 1 0 1 1 0 0 1 2.00 65 0 1 19 0 1 0 1 1 0 0 1 2.00 59 0 1 20 0 1 0 1 1 0 0 1 2.00 39 0 0 20 0 1 0 1 1 0 0 1 2.00 22 1 1 21 0 1 0 1 1 0 0 1 2.00 15 1 1 20 1 1 0 1 1 0 0 1 2.00 7 0 0 25 0 0 0 1 1 0 1 1 2.00 3 0 1 20 0 1 0 1 1 0 0 1 2.00 107 0 1 18 0 1 0 2 2 0 0 1 4.00 114 1 1 19 0 1 0 2 2 0 0 1 4.00 112 0 1 19 0 0 0 2 2 0 1 1 4.00 98 1 1 20 0 1 0 2 2 0 0 1 4.00 94 0 1 20 0 1 0 2 2 0 0 1 4.00 93 0 0 21 0 1 0 2 2 0 0 1 4.00 85 0 1 19 0 1 0 2 2 0 0 1 4.00 73 0 1 19 1 0 0 2 2 0 1 1 4.00 61 1 0 21 0 0 0 2 2 0 1 1 4.00 55 0 0 21 0 1 0 2 2 0 0 1 4.00 49 0 0 19 0 0 0 2 2 0 1 1 4.00 41 0 1 20 0 1 0 2 2 0 0 1 4.00 33 0 1 99 0 1 0 2 2 0 0 1 4.00 21 1 0 19 0 0 0 2 2 0 1 1 4.00 103 0 0 18 0 0 0 3 4 0 1 1 7.90 102 0 1 19 0 0 0 3 4 0 1 1 7.90 101 1 1 19 0 0 0 3 4 0 1 1 7.90 96 1 1 19 0 0 0 3 4 0 1 1 7.90 72 0 1 21 0 0 0 3 4 0 1 1 7.90 71 0 1 20 0 0 0 3 4 0 1 1 7.90 66 0 1 20 0 0 0 3 4 0 1 1 7.90 60 1 1 19 0 0 0 3 4 0 1 1 7.90 47 0 1 19 0 0 0 3 4 0 1 1 7.90 28 1 1 21 0 0 0 3 4 0 1 1 7.90 20 1 1 18 0 0 0 3 4 0 1 1 7.90 16 1 1 20 0 0 0 3 4 0 1 1 7.90 6 0 1 24 0 0 0 3 4 0 1 1 7.90 5 1 0 53 0 0 0 3 4 0 1 1 7.90 111 0 0 18 0 0 0 4 6 0 1 1 11.90 109 0 1 19 0 0 0 4 6 0 1 1 11.90 108 0 1 19 0 0 0 4 6 0 1 1 11.90 99 1 1 20 0 0 0 4 6 0 1 1 11.90 87 0 1 20 0 0 0 4 6 0 1 1 11.90 86 0 1 21 0 0 0 4 6 0 1 1 11.90 80 0 1 20 0 0 0 4 6 0 1 1 11.90 75 0 1 21 1 0 0 4 6 0 1 1 11.90 54 0 1 21 0 0 0 4 6 0 1 1 11.90 51 0 0 20 0 0 0 4 6 0 1 1 11.90 37 0 1 19 0 0 0 4 6 0 1 1 11.90 19 1 0 20 0 0 0 4 6 0 1 1 11.90 115 0 1 41 0 0 0 1 1 1 1 1 2.00 116 1 0 44 0 1 0 2 2 1 0 1 4.00 117 0 0 55 0 0 0 2 2 1 1 1 4.00 118 0 0 50 0 0 0 3 4 1 1 1 7.90 119 0 0 32 0 1 0 0 0 1 0 1 0.00 120 0 0 38 0 1 0 0 0 1 0 1 0.00 121 0 0 39 0 1 0 2 2 1 0 1 4.00 122 0 0 39 0 1 0 0 0 1 0 1 0.00 123 0 1 42 1 0 0 3 4 1 1 1 7.90 124 1 0 46 0 1 0 0 0 1 0 1 0.00 125 1 1 35 0 1 0 3 4 1 0 1 7.90 126 1 0 38 0 0 0 2 2 1 1 1 4.00 127 1 0 37 0 1 0 2 2 1 0 1 4.00 128 0 0 48 0 1 0 0 0 1 0 1 0.00 129 0 0 53 0 1 0 2 2 1 0 1 4.00 130 0 0 38 0 1 0 2 2 1 0 1 4.00 131 0 1 40 0 0 0 3 4 1 1 1 7.90 132 0 0 27 0 1 0 1 1 1 0 1 2.00 133 0 1 38 0 0 0 3 4 1 1 1 7.90 134 0 1 45 0 0 0 2 2 1 1 1 4.00 135 0 1 31 0 0 0 3 4 1 1 1 7.90 136 0 0 29 0 1 0 1 1 1 0 1 2.00 137 0 0 34 0 0 0 2 2 1 1 1 4.00 138 1 1 23 0 0 0 3 4 1 1 1 7.90 139 0 1 39 0 0 0 3 4 1 1 1 7.90 140 0 0 37 0 1 0 0 0 1 0 1 0.00 141 1 1 38 0 0 0 1 1 1 1 1 2.00 142 1 1 26 1 0 0 2 2 1 1 1 4.00 143 0 0 42 0 1 0 2 2 1 0 1 4.00 144 0 1 25 0 1 0 0 0 1 0 1 0.00 145 1 0 30 0 0 0 3 4 1 1 1 7.90 146 0 0 41 0 1 0 0 0 1 0 1 0.00 147 1 1 33 0 0 0 3 4 1 1 1 7.90 148 1 1 45 0 1 0 1 1 1 0 1 2.00 149 0 0 37 1 0 0 1 1 1 1 1 2.00 150 0 0 21 1 0 0 0 0 1 1 1 0.00 151 1 1 47 0 1 0 1 1 1 0 1 2.00 152 0 1 30 0 0 0 2 2 1 1 1 4.00 153 1 0 28 0 0 0 3 4 1 1 1 7.90 154 1 1 23 1 0 0 1 1 1 1 1 2.00 155 0 1 40 0 0 0 2 2 1 1 1 4.00 ; run; proc means; run;