Description Usage Arguments Details Value Examples
This function loops over the list with unit parameters and returns a dataframe with the random coefficients, centered around zero.
1 | ranef(x)
|
x |
The sema output, a list with unit parameters. |
The ranef
function works similar to lme4's
ranef
. The function uses ldply
function from the
plyr
package to extract all coefficients.
A data frame with random effect coefficients per unit, sorted by entry.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 | ## First we create a dataset, consisting of 2500 observations from 20
## units. The fixed effects have the coefficients 1, 2, 3, 4, and 5. The
## variance of the random effects equals 1, 4, and 9. Lastly the
## residual variance equals 4:
test_data <- build_dataset(n = 2500,
j = 20,
fixed_coef = 1:5,
random_coef_sd = 1:3,
resid_sd = 2)
## Next, we fit a simple model to these data
m1 <- sema_fit_df(formula = y ~ 1 + V3 + V4 + V5 + V6 + (1 + V4 + V5 | id),
data_frame = test_data,
intercept = TRUE)
ranef(m1)
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.