Here is another case of ask a big group of statisticians and get a wide range of answers.
The two proposed methods seem to leave something to be desired. It would be nice to have a method that is free of sample size and doesn't require fitting a model with all the higher order interactions. Certainly different strokes for different folks and different types of data require different considerations. In clinical or experimental settings you might be able to elucidate good reasons from clinicians for or not for testing certain interactions, but there can be a lot of observational studies where you may not have as rich clinical thought.
May I suggest an initial screening of your data using trees, either classification or decision or regression tress (whatever you use and/or call them). They have been shown to be very versatile in both a separate analysis of the data and for model building. However in the ANCOVA setting with many variables the trees will help see which categorical/quantitative variables might interact in a regression model.
With the ABCD and XWZ model you said you were looking at then if a tree analysis finds many splits along ABC variables but only X splits when C splits then you may have an interaction worth putting into a regression model. However if there seems to be no method to where XWZ split or the split levels are around the same in different parts of the tree then you may not have any real interactive effects. It's not perfect but it is more exploratory in nature which lets see how tenable a equal slopes assumption is with the data you have.
Good luck!
Jason
-------------------------------------------
Jason Brinkley
East Carolina University
-------------------------------------------
Original Message:
Sent: 05-01-2014 16:47
From: George Milliken
Subject: Checking Model Assumptions for Equal Slopes in ANCOVA - Complicated Model
-------------------------------------------
George Milliken
-------------------------------------------
Beverly: Unfortunately there are no easy short cuts. What I do is put all of the effects in the model, A B A*B ... A*B*C*D and then include X and X interacting with all of the effects, same for Y and Z providing a model with 63 terms. Then I do a backward elimination process with the terms involving the covariates, focusing on the higher order interactions first. Not an easy process, but you can use proc glmselect to do some model building.