[R] syntex about a nested mixed linear model

[JongJoo Kim]

I am a novice R user, and have been in trouble to get the right mixed
model syntax in microarray analyses.

There are three factors: Dye(2 levels), Temperature(3 levels) and
Array(3 for each Temperature with a total of 9 arrays). I want to treat
array as random, and to regard array variation different between
Temperatures.

So the model I want to seek is:

Y = Dye + Temp + Dye + Temp*Dye + Array(Temp) + Dye*Array(Temp).

For F tests of Dye and Dye*Temp, then denominator factor is 
Dye*Array(Temp). For F test of Temp, then denominator factor is
Array(Temp).

I wrote this syntax to fit the above model by typing:

R2.0<-lme(dsnx.0 ~ temp+Dye+temp*Dye, random = list(
Dye*Array(Temp)=pdDiag(~-1+temp*Dye),  Arra=pdDiag(~-1+temp) ) ).

But there is syntax error message.
To obtain one solution, I typed

>int.dye<-interaction(Dye,temp:Arra,drop=T)
which created another factor, int.dye for 'Dye*Array(Temp)', and I
typed:

>R2.0<-lme(dsnx.0 ~ temp+Dye+temp*Dye, random = list(
int.dye=pdDiag(~-1+temp*Dye),  Arra=pdDiag(~-1+temp) ) )

It worked but the degree of freedom for int.dye is 12 instead of 6 <=
(2-1)*3*(3-1).

Could you solve the syntax to fit the right model for this nexted mixed
model?
Thanks so much for your guidance.
Jong
-- 
JongJoo Kim, Ph.D.

243 Shantz Building
Department of Animal Science
University of Arizona
Tucson, AZ 85721
Tel) 520-621-9757
Fax) 520-621-9435
kimjj at email.aizona.edu