Displaying 1 result from an estimated 1 matches for "cmbeale".
Did you mean:
beale
2007 May 25
0
Help with complex lme model fit
...the correct results. It seems there might be
some
identifiability issues that could possibly be dealt with by
specifying
starting parameters - but I can't see how to do this. I'm comparing
results from R to those got when using GenStat...
The raw data are available on the web at
http://cmbeale.freehostia.com/OutData.txt and can be read directly
into R using:
gpdat <- read.table("http://cmbeale.freehostia.com/OutData.txt",
header = T)
gpdat$X7 <- as.factor(gpdat$X7)
gpdat$X4 <- as.factor(gpdat$X4)
rand_mat <- as.matrix(gpdat[,11:26])
gpdat <- groupedData(Y1...