Displaying 1 result from an estimated 1 matches for "m_sum".
Did you mean:
i_sum
2007 May 25
0
Help with complex lme model fit
...t 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 ~X1 + X2 + X3 + X4 + X5 + m_sum|.g, data =
gpdat)
the model fitted using:
library(Matrix)
library(nlme)
m_sum <- rowSums(gpdat[,11:27])
mod1 <- lme(fixed = Y1 ~ X1 + X2 + X3 + X4 + X5 + m_sum,
random = pdBlocked(list(pdIdent(~1), pdIdent (~ X6
-
1),
pdIdent (~...