/* Combability Measurements on Human Hair. Garcia and Diaz.
J. Soc. Cosmet. Chem. 27,379-398 (1976)  */

options nodate nonumber ps=54 ls=76;

data one;
do swatch=1 to 16; 
do trt=0 to 1;
do run=1 to 5;
if swatch <= 8 then shampoo=1; else shampoo=2;
input combability @@;
output;
end; end; end;
cards;
223 155 160 193 173 183 133 190 153 173
228 333 310 233 143 173 173 115 198 150
103  90  58  65  50  80  75  68  70  58
123 143 170 100 148 115 125 120 125 148
115 158 110 108 138 145 113  98 138 140
160 128 140 135  88  73 123 100  75 115
 90 170 125 193 213 123 138 100 253  93
 60  50  40  50  50  38  55  35  38  53
138 113 178 143 150 255 110 195  93 213
173 148 165 175 108 118 200 145 155 108
170 165 215 268 190 133 155 145 240 230
120 130 120 150 140 150 130 110 185 105
100  70  95  90  95  65  60 105  90 100
 80  60  85  80  75  95  65  45  65  60
100 100 115  80 105 145 175 125 145 180
110  80  65  85 100  60  55  60  45  65
;
run;

data two;
set one;
file 'H:\public_html\data\combhair.dat';
put @1 shampoo 1.0 +2 trt 1.0 +2 swatch 2.0 +2 run 1.0 +2 combability 3.0;
run;

quit;

/* 

proc print;
run;


proc glm;
by shampoo;
class trt swatch;
model combability=trt swatch trt*swatch;
run;

quit;