Displaying 1 result from an estimated 1 matches for "sumlmer".
Did you mean:
sumlme
2010 Mar 22
0
using lmer weights argument to represent heteroskedasticity
...- function(o){
tab <- summary(o)$tTable
r <- c(tab[1,1:2], tab[2,1:2], as.numeric(VarCorr(o))[3:4])
names(r) <- c("b0","b0sd","b1","b1sd","sdgroup","sdresid")
return(r)
}
## function to collect pieces from lmer() output
sumLMER <- function(o){
tab <- summary(o)@coefs
r <- c(tab[1,1:2], tab[2,1:2],
as.numeric(attributes(VarCorr(o)$groupid)$stddev),
as.numeric(attr(VarCorr(o), "sc")))
names(r) <- c("b0","b0sd","b1","b1sd","sdgroup","sdresi...