BookmarkSubscribeRSS Feed
🔒 This topic is solved and locked. Need further help from the community? Please sign in and ask a new question.
plf515
Lapis Lazuli | Level 10

Using SAS 9.4.

 

16 animals were randomly assigned to two treatments (control and treatment). Each animal had its level on the dependent variable measured 24 times over the course of the day. For each animal, these numbers vary wildly, from 0 into the thousands, with no strong correlations.

 

Interest is in whether the two treatments are different.

 

I compared three models:

a) T-test on total activity level

proc ttest data = wide;
 class group;
 var total;
run;

 

b) A multilevel model with unstructured covariances

title "MLM for day 1, unstructured covariance, only group entered";
proc mixed data = long;
 class group mouse_id;
 model act = group/solution;
 repeated/type = un subject = mouse_id;
run;

c) A multilevel model with AR(1) covariances

 

title "MLM for day 1, AR1 covariance, only group entered";
proc mixed data = long;
 class group mouse_id;
 model act = group/solution;
 repeated/type = ar(1) subject = mouse_id;
run;

c) Had slightly better AIC and BIC than b).

 

The effect size of group was roughly equal across all three models (after accounting for the fact that the t-test was on 24 values added up) but the standard errors for b) and c) were wildly different. b) Was significant, a) and c) were not.

 

In the past, I've not seen such huge changes in the standard errors.

 

Any thoughts on why this might occur and what to do?

 

NOTE: It's the same data, just one in long format and one in wide format.

1 ACCEPTED SOLUTION

Accepted Solutions
PaigeMiller
Diamond | Level 26

There's not much we can say since we don't have your data, and you have changed the data sets from one analysis to another anyway (not even sure how you can compare different data sets in this way).

 

But ... if you fit a model with the wrong covariance structure (which doesn't match the data), I can understand why the model doesn't fit well and the standard errors are large.

--
Paige Miller

View solution in original post

1 REPLY 1
PaigeMiller
Diamond | Level 26

There's not much we can say since we don't have your data, and you have changed the data sets from one analysis to another anyway (not even sure how you can compare different data sets in this way).

 

But ... if you fit a model with the wrong covariance structure (which doesn't match the data), I can understand why the model doesn't fit well and the standard errors are large.

--
Paige Miller

sas-innovate-2024.png

Join us for SAS Innovate April 16-19 at the Aria in Las Vegas. Bring the team and save big with our group pricing for a limited time only.

Pre-conference courses and tutorials are filling up fast and are always a sellout. Register today to reserve your seat.

 

Register now!

What is ANOVA?

ANOVA, or Analysis Of Variance, is used to compare the averages or means of two or more populations to better understand how they differ. Watch this tutorial for more.

Find more tutorials on the SAS Users YouTube channel.

Discussion stats
  • 1 reply
  • 1151 views
  • 1 like
  • 2 in conversation