Displaying 1 result from an estimated 1 matches for "bmdp8v".
Did you mean:
bmdp
2001 Aug 09
0
converting BMDP 8V mixed model treatment to R ?
...nd.
/finish.
I tried in R (some minor differences in variable naming):
latm <- read.table( "latm.dat" )
latm.aov <- aov( llat ~ direct * flicker *ground *retest +person
+ Error( person : ( direct * flicker *ground *retest )))
Indeed, summary(latm.aov) gives the same F tests than BMDP8V for the 1st part
of the evaluation. But I don't get variance components from this.
So it looks as if I have to change to mle from nmle package. I did not succeed
to translate this model to mle yet.
Trial-and-error fails mostly, because
the problem thrashes my system (FreeBSD, K6-233, 64MB). T...