Hello,
Am getting a warning from my 'Heterogenous Compound Symmetry' code and ods output will not generate the output (AIC, BIC etc.) I wanted. I have used the same code for other covariance structures and it worked perfectly.
I used the same code I posted below for another dataset with CSH and it worked. Here is a sample of my dataset. I would appreciate any help. Thanks
data have;
input Obs ID Group pre post1 post2 post3 time wk @@;
cards;
1 14789 2 10 10 15 15 1 10
2 14789 2 10 10 15 15 2 10
3 14789 2 10 10 15 15 3 15
4 14789 2 10 10 15 15 4 15
5 15170 2 16 20 20 20 1 16
6 15170 2 16 20 20 20 2 20
7 15170 2 16 20 20 20 3 20
8 15170 2 16 20 20 20 4 20
9 15611 1 10 20 18 22 1 10
10 15611 1 10 20 18 22 2 20
11 15611 1 10 20 18 22 3 18
12 15611 1 10 20 18 22 4 22
13 15697 1 10 10 10 11 1 10
14 15697 1 10 10 10 11 2 10
15 15697 1 10 10 10 11 3 10
16 15697 1 10 10 10 11 4 11
17 15844 1 21 26 26 27 1 21
18 15844 1 21 26 26 27 2 26
19 15844 1 21 26 26 27 3 26
20 15844 1 21 26 26 27 4 27
21 15846 2 26 29 28 24 1 26
22 15846 2 26 29 28 24 2 29
23 15846 2 26 29 28 24 3 28
24 15846 2 26 29 28 24 4 24
25 15998 2 10 10 10 12 1 10
26 15998 2 10 10 10 12 2 10
27 15998 2 10 10 10 12 3 10
28 15998 2 10 10 10 12 4 12
29 16021 1 11 13 12 15 1 11
30 16021 1 11 13 12 15 2 13
31 16021 1 11 13 12 15 3 12
32 16021 1 11 13 12 15 4 15
33 16961 2 14 20 18 18 1 14
34 16961 2 14 20 18 18 2 20
35 16961 2 14 20 18 18 3 18
36 16961 2 14 20 18 18 4 18
37 17129 2 14 16 18 15 1 14
38 17129 2 14 16 18 15 2 16
39 17129 2 14 16 18 15 3 18
40 17129 2 14 16 18 15 4 15
41 17571 2 17 21 23 . 1 17
42 17571 2 17 21 23 . 2 21
43 17571 2 17 21 23 . 3 23
44 17571 2 17 21 23 . 4 .
45 17578 1 21 40 37 43 1 21
46 17578 1 21 40 37 43 2 40
47 17578 1 21 40 37 43 3 37
48 17578 1 21 40 37 43 4 43
49 17939 1 16 13 14 14 1 16
;
run;
SAS code:
/* Heterogenous Compound Symmetry CSH */
proc mixed data=have;
class id group time;
model wk=id time time*group;
repeated time / type = csh subject = id r rcorr;
ods output FitStatistics=FitCSH(rename=(value=CSH))
Dimensions=ParmCSH(rename=(value=NumCSH));
run;
@pink_poodleI found the solution. Just sort the data by ID before running the proc mixed and exclude the ID from the class statement.
@pink_poodle Am assuming you didn't get the AIC as part of your output but if you did let me know. I did not get AIC and this is one of the output am looking for.
@pink_poodleI found the solution. Just sort the data by ID before running the proc mixed and exclude the ID from the class statement.
Are you ready for the spotlight? We're accepting content ideas for SAS Innovate 2025 to be held May 6-9 in Orlando, FL. The call is open until September 25. Read more here about why you should contribute and what is in it for you!
Learn how use the CAT functions in SAS to join values from multiple variables into a single value.
Find more tutorials on the SAS Users YouTube channel.