извлечение фиксированных эффектов и стандартных ошибок из нескольких объектов Ime в R

Question

Мне интересно, есть ли какой-нибудь умный способ автоматически извлекать несколько различных фиксированных эффектов (а также их стандартные ошибки) из серии объектов lme.Я хотел бы создать что-то вроде латексной таблицы из этих выходных данных, чтобы результаты можно было сравнить друг с другом.Что-то вроде функциональности пакета apsrtable.Спасибо!Антонио.

Answer 1

Это то, что вы имеете в виду?

library(nlme)
## from ?lme:
fm1 <- lme(distance ~ age, data = Orthodont) # random is ~ age
fm2 <- lme(distance ~ age + Sex, data = Orthodont, random = ~ 1)
L <- list(A=fm1,B=fm2)
lapply(L, function(x) summary(x)$tTable)

Answer 2

См. ответ Алана Манна на мой вопрос того же рода на tex.stackexchange.com.

\documentclass{article}
\usepackage[utf8]{inputenc}
\usepackage{booktabs}
\usepackage{caption}
\title{Side-by-side xtables}
\author{}
\date{}
\begin{document}
\maketitle
First some R code to create some data.
<<>>=
myData <- matrix(c(19,89,23,23,74,44,16,39,67),ncol=3,byrow=TRUE)
colnames(myData) <- c("A","B","C")
rownames(myData) <- c("1","2","3")
myData2 <- myData * 2
@

Now we place the data in two side-by-side tables:

\begin{table}[htb]
\begin{minipage}{.45\textwidth}
\centering
<<echo=FALSE,results=tex>>=
library("xtable")
print(xtable(myData),
  floating=FALSE,
  hline.after=NULL,
  add.to.row=list(pos=list(-1,0, nrow(myData)),
  command=c('\\toprule\n','\\midrule\n','\\bottomrule\n')))
@
\captionof{table}{The first table}
\end{minipage}
\begin{minipage}{.45\textwidth}
\centering
<<echo=FALSE,results=tex>>=
print(xtable(myData2),
  floating=FALSE,
  hline.after=NULL,
  add.to.row=list(pos=list(-1,0, nrow(myData2)),
  command=c('\\toprule\n','\\midrule\n','\\bottomrule\n')))
@
\captionof{table}{The second table}
\end{minipage}
\end{table}
\end{document}

Answer 3

Попробуйте это

fm.lmer <- lmer(formula, data=Data)
fm.Means <- data.frame("Mean"=fixef(fm.lmer)[-1])
fm.SEs <- data.frame("SEM"=sqrt(diag(vcov(fm.lmer)))[-1])
fm.MeansSEs <- data.frame("Entry"=gsub('Entry', '', rownames(fm.Means), 
                          fixed = TRUE),"Mean"=fm.Means, "SEM"=fm.SEs)

Другим выбором может быть mtable-lme4 function.