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.
SAS Innovate 2025 is scheduled for May 6-9 in Orlando, FL. Sign up to be first to learn about the agenda and registration!
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.
Ready to level-up your skills? Choose your own adventure.