1 |
fit |
|
form |
|
varname |
|
varpattern |
|
data |
|
... |
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 | ##---- Should be DIRECTLY executable !! ----
##-- ==> Define data, use random,
##-- or do help(data=index) for the standard data sets.
## The function is currently defined as
function (fit, form, varname, varpattern = vname, data = fit@frame,
...)
{
model.mat <- model.matrix(form, data, ...)
ret <- data[rownames(model.mat), varname, drop = F]
fe <- fixef(fit)
effect.names <- grep(varpattern, names(fe), value = T)
ret$comp <- c(model.mat[, effect.names] %*% fe[effect.names])
attr(ret, "predictors") <- effect.names
ret
}
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.