Heterogenous Compound Symmetry

CathyVI — Wed, 22 Mar 2023 18:34:39 GMT

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.

Re: Heterogenous Compound Symmetry

pink_poodle — Wed, 22 Mar 2023 20:13:07 GMT

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

Re: Heterogenous Compound Symmetry

CathyVI — Fri, 24 Mar 2023 18:39:12 GMT

@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.

Re: Heterogenous Compound Symmetry

pink_poodle — Mon, 27 Mar 2023 17:50:43 GMT

Hi @CathyVI ,

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

Re: Heterogenous Compound Symmetry

CathyVI — Mon, 27 Mar 2023 17:55:40 GMT

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

topic Re: Heterogenous Compound Symmetry in SAS Programming

Heterogenous Compound Symmetry

Re: Heterogenous Compound Symmetry

Re: Heterogenous Compound Symmetry

Re: Heterogenous Compound Symmetry

Re: Heterogenous Compound Symmetry