从R中的几个lme对象中提取固定效果和标准误差

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

试试这个

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功能。

Answer 3

在tex.stackexchange.com上查看Alan Munn's对我同一个问题的回答。

\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}