Displaying 1 result from an estimated 1 matches for "subjectname1".
Did you mean:
subjectname
2018 Jan 03
0
start values for random effects in nlme
...others in future.
I have been struggling with the same issue. When analysing the source code (https://svn.r-project.org/R-packages/trunk/nlme/R/nlme.R)
I found out that the E matrix needs rownames referring to the random effect of each group in the dataset, something like:
rownames(E) <- c("subjectname1","subjectname2","subjectname3",.....).
Then your code should run perfectly OK.
Example:
----------
library(nlme)
data <- data.frame(a=c(1,2,3,4,5,1,2,3,4,5,1,2,3,4,5),
b=c(2,3,4,5,7,3,4,5,7,8,3,5,6,8,9),
ID=c(rep('X',5),rep...