A variable is considered to be nested when it is unique within that
combination. i.e. if your measuring Apparatus was different at each sample.
Hence, as Lorenz suggested it is reasonable to treat Apparatus as fixed,
specially if you are considering making comparisons between the 2 apparatus.
I hope that this helps
Francisco
From: [EMAIL PROTECTED]
To: [EMAIL PROTECTED]
Subject: RE: [R] grouping for lme with nested repeated measurements
Date: Mon, 25 Oct 2004 08:11:21 +0200
> 1.The simple problem ist that i have different Samples , from
> which i make repeated measurements (each sample is measured 6
> times) and i repeat this experiment over several Days, so i
> get the lme grouping term "random=~1|Days/Sample".
I would rather specify this as: random= ~ 1 | Sample/Days, but I
am not quite sure how this affects the model.
> 2. Now i am measuring with 2 different measuring Apparatus
> the same Sample each 6 times, to see how big the difference
> from the appratus is.
> Because Apparatus is on the level of repeated measurements i
> cant write Days/Sample/Apparatus. the lme function offers a
> list() feature to design the grouping, but i didnt understand
> this, if it is the solution to the problem.
Why would you want to include the Apparatus in the random effect?
I assume that you are interested in differences and thus, this
is a fixed effect:
lme (fixed= response ~ apparatus, data= XX, random= ~ 1 | Sample/Days)
Lorenz
-
Lorenz Gygax, [EMAIL PROTECTED]
Centre for proper housing of ruminants and pigs
Swiss Federal Veterinary Office
______________________________________________
[EMAIL PROTECTED] mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide!
http://www.R-project.org/posting-guide.html
______________________________________________
[EMAIL PROTECTED] mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
