-
Notifications
You must be signed in to change notification settings - Fork 0
Open
Description
In #6, there was a wrapper added to perform the extraction. However, we extracted a bit too many terms. In particular, we only want to extract elements corresponding to the random-effects, not all effects (p
). c.f.
Line 152 in f649ced
T1 <- VarCorr(x)[[group_term]][1:p,1:p] |
Consider revising given:
Lines 145 to 152 in f649ced
# Figure out the class term: | |
term_labels <- attr(mf, "names") | |
# Retrieve the _LAST_ variable (assume it is a group variable.) | |
# To Do: What is the best way to make this approach robust ? | |
group_term <- term_labels[length(term_labels)] | |
T1 <- VarCorr(x)[[group_term]][1:p,1:p] |
Issue arose with:
m <- lmer(DV ~ P + P:H + (1+P||id), data = data)
No data was given, so an equivalent test case should be:
fm1 <- lmer(Reaction ~ Days + (1+Days || Subject), data = sleepstudy)
Metadata
Metadata
Assignees
Labels
No labels