Displaying 3 results from an estimated 3 matches for "varcor".
Did you mean:
arcor
2013 Nov 12
1
Getting residual term out of lmer summary table
...1); output # ok, there is the variance component of the random effect group and the residual term
## Now I'd like to get the variance components of the random effect "Group" and of the residual term "Residual" in order
# to do further calculations with these numbers
output$varcor[1] ## reveals the variance of the random effect "Group"
output$varcor[2] ## does not reveal the residual term! what other command do I need to use then?
2008 Aug 29
3
extract variance components
HI,
I would like to extract the variance components estimation in lme function
like
a.fit<-lme(distance~age, data=aaa, random=~day/subject)
There should be three variances \sigma_day, \sigma_{day %in% subject } and
\sigma_e.
I can extract the \sigma_e using something like a.fit$var. However, I cannot
manage to extract the first two variance components. I can only see the
results in
2017 Dec 26
1
identifying convergence or non-convergence of mixed-effects regression model in lme4 from model output
....m)$devcomp$dims[1])
? #add model name
? mod.data$model = name
? return(mod.data)
}
lmer.ranef.data.extract = function(lmer.mod,
name=deparse(substitute(lmer.mod))){
? #extract random effect variance, standard error, correlations between
slope and intercept
? mod.data.ef = as.data.frame(VarCorr(lmer.mod))
? mod.data.ef$n.subj = as.numeric(summary(lmer.mod)$ngrps[1]) #number
of subjects
? mod.data.ef$n.item = as.numeric(summary(lmer.mod)$ngrps[2]) #number
of items
? #add number of datapoints
? mod.data.ef$N = as.numeric(summary(incr.best.m)$devcomp$dims[1])
? #add model name
?...