What's your error message? I see a syntax error in
random=list(Probe=pdBlocked(list(~1,pdCompSymm(~1
),End=~1,ProbeNo=pdCompSymm(~1))
I count 5 open parens and 3 close parens. I know that might have
been part of what you copied into this email and not what you gave R.
Dear R community,
I have trouble implementing a nested variance-covariance matrix in the
lme function.
The model has two fixed effects called End and logpgc, the response
variable is the logarithm to base 2 of Intensity ( log2(Intensity) )
and the random effects are called Probe and ProbeNo.
I completely forgot to supply the R code I tried:
vov1i2-read.table(VOV1_INHIBITED6-16-2006-13h35min33sec.txt,header=TRUE)
test.lme-lme(fixed=log2(Intensity)~End+logpgc,random=list(Probe=pdBlocked(list(~1,pdCompSymm(~1
),End=~1,ProbeNo=pdCompSymm(~1)),data=vov1i2)
It doesn't look right to me