Thanks for your input. Maybe I didn't make myself clear enough re
"age range". Age range refers to actual age in days at each
measurement. As all subjects were measured on the same day at each
measurement event, the difference between subject ages is the same
at each measurement event, so age at first measurement would be the
data in the "age range" variable. Do you still think this should be
inlcuded as a categorical variable?
________________________________________
From: [email protected] [[email protected]
] on behalf of Joseph Coveney [[email protected]]
Sent: Friday, 15 May 2009 11:48 AM
To: [email protected]
Subject: st: RE: Covariates in ANOVAs
Miss Gina Micke wrote:
I would like to know how to control for an independent continuous
variable when
using a repeated measures ANOVA. I have a 2x2 factorial design
experiment with
males and females in each of the treatment groups (T1 & T2) measured
at multiple
"events", each subject assigned an "id". The subjects are all
different ages so
I would like to be able to control for the effect of "age range".
If I didn't have "age range" as a covariate my model would be as
follows:
anova x T1 T2 sex T1*T2 T1*sex T2*sex T1*T2*sex/id|T1*event T2*event
T1*T2*event
sex*event sex*T1*event sex*T2*event T1*T2*sex*event, repeated(event)
partial
bse(id|T1*T2*sex)
Do I use the ANCOVA and cont(age range) command? Is this output the
effects of
my fixed factors taking into account the effect of age range?
--------------------------------------------------------------------------------
First, the command line text shown doesn't look at all like correct
syntax.
Second, by the name of it, a variable called "age range" would be
categorical
and so you would not be using the -continuous()- option for it in -
anova-.
Unless you have balanced data (equal representation of age ranges
within each
factorial cell, as well as equal sex-by-treatment group cell sizes
maintained
throughout the observation period), try something like the following
to start
with:
xi3 i.treatment_group*i.age_group*i.sex*i.observation_interval
xtmixed response _I* || patient_id: , nolrtest variance reml
You'll need to install the user-written command, -xi3-, from SSC if
you haven't
already.
Joseph Coveney
*
* For searches and help try:
* http://www.stata.com/help.cgi?search
* http://www.stata.com/support/statalist/faq
* http://www.ats.ucla.edu/stat/stata/
*
* For searches and help try:
* http://www.stata.com/help.cgi?search
* http://www.stata.com/support/statalist/faq
* http://www.ats.ucla.edu/stat/stata/
