BookmarkSubscribeRSS Feed
deleted_user
Not applicable
Hi there,

A fellow researcher and I were discussing stats and he indicated that when asked to review papers for publication (we are both livestock researchers) he likes to see ANOVA tables. He likes to do this especially when several effects are in the P = 0.05-0.1 range to ensure that the error terms are not inflated and a Type 2 error (failing to reject the null hypothesis) is not being committed.

I analyze all of my experiments using PROC MIXED (where appropriate) and was wondering whether there is some way to get it to include the ANOVA tables in the output, in the event that I am ever asked to produce them to a reviewer. It seems to me that PROC GLM does this automatically and was wondering if it was as simple as just adding a line of code to my programs.

If anyone could offer a suggestion it would be much appreciated.
4 REPLIES 4
statsplank
Calcite | Level 5
Hi Matt O,

You might have noticed a "Type 3 Tests of Fixed Effects" table in the output produced by PROC MIXED. This table contains Type III F-statistics and p-values for the variables (fixed effects) in your model and can be considered an analog to an ANOVA table that PROC GLM produces. There are no sum of squares or mean sum of squares there because PROC MIXED uses maximum-likelihood estimation methods.
deleted_user
Not applicable
Hi Matt O,

If in the PROC MIXED statement you add METHOD=TYPE3, for example, you'll get expected mean squares and a Type 3 ANOVA table. This, however, does not work if you have subject effects, or a repeated statement. You can also use METHOD=TYPE1 or METHOD=TYPE2 depending on the type of sum of squares you want. This should match PROC GLM.

Hope this helps.

j.
Dale
Pyrite | Level 9
Yes, adding the METHOD=TYPE3 option to the PROC MIXED statement is the key to getting an ANOVA table from the MIXED procedure. It should be noted that doing so results in moment estimates of the parameters of the model. The moment estimators and the REML estimators may be identical for some balanced designs. However, for unbalanced designs, moment and REML estimators are usually not identical. Be forewarned that if you redo some analyses using moment estimators which were already performed using REML estimation, then some of those marginally significant effects could show different p-values.

D.

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
  • 4 replies
  • 4396 views
  • 0 likes
  • 3 in conversation