/*TYPE        Type of companion
0: newly pregnant female
1: virgin female
9: not applicable (when PARTNERS=0)*/


options ls=75;

proc format;
   value trtfmt
   1 = "control"
   2 = "1 preg."
   3 = "1 rec."
   4 = "8 preg."
   5 = "8 rec.";
   value typefmt
   0 = "pregnant"
   1 = "receptive"
   9 = "na";
run;

data flies;
   format treatment trtfmt.;
   format typecomp typefmt.;
   input id numcomps typecomp days thorax slpng treatment;
   tmp=((slpng/100)**(0.5));
   asleep=arsin(tmp);
   if (treatment>1) then companions="yes";
   else companions="no";
   control=(companions="no");
   cards;
01 8 0 35 0.64 22  4
02 8 0 37 0.68  9  4
03 8 0 49 0.68 49  4
04 8 0 46 0.72  1  4
05 8 0 63 0.72 23  4
06 8 0 39 0.76 83  4
07 8 0 46 0.76 23  4
08 8 0 56 0.76 15  4
09 8 0 63 0.76  9  4
10 8 0 65 0.76 81  4
11 8 0 56 0.80 12  4
12 8 0 65 0.80 15  4
13 8 0 70 0.80 37  4
14 8 0 63 0.84 24  4
15 8 0 65 0.84 26  4
16 8 0 70 0.84 17  4
17 8 0 77 0.84 14  4
18 8 0 81 0.84 14  4
19 8 0 86 0.84  6  4
20 8 0 70 0.88 25  4
21 8 0 70 0.88 18  4
22 8 0 77 0.92 26  4
23 8 0 77 0.92 24  4
24 8 0 81 0.92 29  4
25 8 0 77 0.94 27  4
01 0 9 40 0.64 18  1
02 0 9 37 0.70  6  1
03 0 9 44 0.72 19  1
04 0 9 47 0.72  7  1
05 0 9 47 0.72 16  1
06 0 9 47 0.76 13  1
07 0 9 68 0.78 35  1
08 0 9 47 0.80  2  1
09 0 9 54 0.84 35  1
10 0 9 61 0.84  6  1
11 0 9 71 0.84 15  1
12 0 9 75 0.84 14  1
13 0 9 89 0.84 18  1
14 0 9 58 0.88 50  1
15 0 9 59 0.88 25  1
16 0 9 62 0.88 10  1
17 0 9 79 0.88 33  1
18 0 9 96 0.88 43  1
19 0 9 58 0.92 35  1
20 0 9 62 0.92 17  1
21 0 9 70 0.92 27  1
22 0 9 72 0.92 22  1
23 0 9 75 0.92 16  1
24 0 9 96 0.92 20  1
25 0 9 75 0.94 37  1
01 1 0 46 0.64 23  2
02 1 0 42 0.68  4  2
03 1 0 65 0.72 20  2
04 1 0 46 0.76 42  2
05 1 0 58 0.76  9  2
06 1 0 42 0.80 32  2
07 1 0 48 0.80 66  2
08 1 0 58 0.80 28  2
09 1 0 50 0.82 10  2
10 1 0 80 0.82  4  2
11 1 0 63 0.84 12  2
12 1 0 65 0.84 17  2
13 1 0 70 0.84 12  2
14 1 0 70 0.84 23  2
15 1 0 72 0.84 40  2
16 1 0 97 0.84 18  2
17 1 0 46 0.88 10  2
18 1 0 56 0.88 38  2
19 1 0 70 0.88  7  2
20 1 0 70 0.88 23  2
21 1 0 72 0.88 36  2
22 1 0 76 0.88  9  2
23 1 0 90 0.88 21  2
24 1 0 76 0.92 62  2
25 1 0 92 0.92 36  2
01 1 1 21 0.68 23  3
02 1 1 40 0.68 62  3
03 1 1 44 0.72 28  3
04 1 1 54 0.76 18  3
05 1 1 36 0.78 10  3
06 1 1 40 0.80 28  3
07 1 1 56 0.80 22  3
08 1 1 60 0.80 29  3
09 1 1 48 0.84 15  3
10 1 1 53 0.84 73  3
11 1 1 60 0.84 10  3
12 1 1 60 0.84  5  3
13 1 1 65 0.84 13  3
14 1 1 68 0.84 27  3
15 1 1 60 0.88 20  3
16 1 1 81 0.88 21  3
17 1 1 81 0.88 12  3
18 1 1 48 0.90 49  3
19 1 1 48 0.90 17  3
20 1 1 56 0.90 22  3
21 1 1 68 0.90 71  3
22 1 1 75 0.90 17  3
23 1 1 81 0.90 10  3
24 1 1 48 0.92 24  3
25 1 1 68 0.92 18  3
01 8 1 16 0.64 34  5
02 8 1 19 0.64  6  5
03 8 1 19 0.68  4  5
04 8 1 32 0.72 22  5
05 8 1 33 0.72 28  5
06 8 1 33 0.74 31  5
07 8 1 30 0.76 16  5
08 8 1 42 0.76 27  5
09 8 1 42 0.76  8  5
10 8 1 33 0.78 32  5
11 8 1 26 0.80 20  5
12 8 1 30 0.80 35  5
13 8 1 40 0.82 12  5
14 8 1 54 0.82 14  5
15 8 1 34 0.84 17  5
16 8 1 34 0.84 29  5
17 8 1 47 0.84 31  5
18 8 1 47 0.84  6  5
19 8 1 42 0.88 30  5
20 8 1 47 0.88 27  5
21 8 1 54 0.88 40  5
22 8 1 54 0.88 19  5
23 8 1 56 0.88  8  5
24 8 1 60 0.88  8  5
25 8 1 44 0.92 15  5
;
run;

proc glm;
   title "nested analysis of fruitflies";
   class control numcomps typecomp;
   model days=         ;
   lsmeans          ;
run;