BookmarkSubscribeRSS Feed
☑ This topic is solved. Need further help from the community? Please sign in and ask a new question.
CathyVI
Pyrite | Level 9

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.
1 ACCEPTED SOLUTION

Accepted Solutions
CathyVI
Pyrite | Level 9

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

View solution in original post

4 REPLIES 4
pink_poodle
Barite | Level 11
I ran this code, and it shows an output but includes a warning: "Stopped because of infinite likelihood."
CathyVI
Pyrite | Level 9

@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
Barite | Level 11

Hi @CathyVI ,

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

pink_poodle_0-1679939417708.png

 

CathyVI
Pyrite | Level 9

@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: Save the Date

 SAS Innovate 2025 is scheduled for May 6-9 in Orlando, FL. Sign up to be first to learn about the agenda and registration!

Save the date!

How to Concatenate Values

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.

SAS Training: Just a Click Away

 Ready to level-up your skills? Choose your own adventure.

Browse our catalog!

Discussion stats
  • 4 replies
  • 826 views
  • 4 likes
  • 2 in conversation