Well, I kind of feel like an idiot right now, and just deleted everything I had previously typed in response... While I was typing this up: "Just to clarify the data, each "pig" is an individual animal (the number is a unique identifier), and each pig is only in one treatment, it never received more than one treatment." And that's when I realized that one of the notes that I got in the log while trying your MIXED / Repeated code said this: "NOTE: An infinite likelihood is assumed in iteration 0 because of a nonpositive definite estimated R matrix for pig 30-4." Which is when I realized that TWO pigs were labeled as "30-4", so it registered two t0, t1, t2 for that animal. I fixed that, deleted the line "singular=1e-9", and then the problem with the "residual" went away. I now have output separated by treatment and by time. And my appologies for wasting your time/energy due to a simple typo, but I do appreciate the enormous help figuring out the Slice code, I don't even want to know how long it would have taken me to figure that out myself. As a biologist, I wish I'd taken more statistics courses and gotten more experience with the software, rather than mucking around with pigs all day! Data DSS;
Title 'DSS DSS-hLZ RBC';
Input pig $ Trt $ Time Resp;
Cards;
28-5 TRT1 0 7.08
28-8 TRT1 0 7.34
29-12 TRT1 0 8.09
29-14 TRT1 0 7.98
30-4a TRT1 0 7.34
30-4b TRT1 0 8.07
30-7 TRT1 0 8.41
28-6 TRT2 0 7
29-10 TRT2 0 6.9
30-6 TRT2 0 8.8
31-1 TRT2 0 7.08
28-1 TRT3 0 7.61
28-7 TRT3 0 7.57
30-5 TRT3 0 7.23
31-6 TRT3 0 6.6
28-3 CTRL 0 7.16
28-9 CTRL 0 7.54
30-3 CTRL 0 6.66
31-4 CTRL 0 8.28
31-5 CTRL 0 6.68
28-5 TRT1 1 7.11
28-8 TRT1 1 8.16
29-12 TRT1 1 8.25
29-14 TRT1 1 7.47
30-4a TRT1 1 7.86
30-4b TRT1 1 8.06
30-7 TRT1 1 7.89
28-6 TRT2 1 7.17
29-10 TRT2 1 7.54
30-6 TRT2 1 9.75
31-1 TRT2 1 7.66
28-1 TRT3 1 7.6
28-7 TRT3 1 7.26
30-5 TRT3 1 7.45
31-6 TRT3 1 6.4
28-3 CTRL 1 7.37
28-9 CTRL 1 7.84
30-3 CTRL 1 7
31-4 CTRL 1 8.08
31-5 CTRL 1 7.18
28-5 TRT1 2 7.09
28-8 TRT1 2 8.19
29-12 TRT1 2 8.28
29-14 TRT1 2 7.25
30-4a TRT1 2 6.67
30-4b TRT1 2 6.7
30-7 TRT1 2 7.73
28-6 TRT2 2 7.49
29-10 TRT2 2 7.05
30-6 TRT2 2 7.89
31-1 TRT2 2 6.47
28-1 TRT3 2 7.27
28-7 TRT3 2 7
30-5 TRT3 2 7.61
31-6 TRT3 2 6.38
28-3 CTRL 2 6.81
28-9 CTRL 2 7.54
30-3 CTRL 2 6.32
31-4 CTRL 2 7.06
31-5 CTRL 2 6.84
;
Proc GLIMMIX Data=DSS;
Class pig trt time;
Model Resp = trt time trt*time;
Random time / subject=pig type=ar(1) residual;
slice time*trt / sliceby=time sliceby=trt lines;
LSMEANS trt time/ diff lines;
Run;
... View more