Heterogenous Compound Symmetry

[CathyVI]

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;

WARNING: Output 'FitStatistics' was not created. Make sure that the output object name, label, or path is spelled correctly.

Also, verify that the appropriate procedure options are used to produce the requested output object. For example, verify that the NOPRINT option is not used.

[pink_poodle]

I ran this code, and it shows an output but includes a warning: "Stopped because of infinite likelihood."

[CathyVI]

@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_poodle]

Hi @CathyVI ,

You are right, the code does not output AIC. The last table it outputs is this:

[CathyVI · Accepted solution]

@pink_poodleI found the solution. Just sort the data by ID  before running the proc mixed and  exclude the ID from the class statement.