glht输出到R中乳胶的结果

glht输出到R中乳胶的结果,r,latex,R,Latex,我想将R中的glht对象的结果导出到LaTeX表中 例如,使用“stargazer”库,可以生成一个lme对象的格式非常好的LaTeX表 我希望从glht对象的摘要输出中自动创建一个LaTeX表,例如使用创建的摘要 >summary(glht(dataModel)) Linear Hypotheses: Estimate Std. Erro

我想将R中的
glht
对象的结果导出到LaTeX表中

例如,使用“stargazer”库,可以生成一个
lme
对象的格式非常好的LaTeX表

我希望从
glht
对象的摘要输出中自动创建一个LaTeX表,例如使用创建的摘要

>summary(glht(dataModel))
Linear Hypotheses:
                                                                        Estimate Std.     Error z value Pr(>|z|)    
Group1 - Group2 == 0   -0.14007    0.01589  -8.813   <0.001 "***"
Group1 - Group3 == 0    -0.09396    0.01575  -5.965   <0.001 ***
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 
(Adjusted p values reported -- single-step method)
>摘要(glht(数据模型))
线性假设:
估计标准误差z值Pr(>z)

Group1-Group2==0-0.14007 0.01589-8.813R-help的Rich给出了帮助线索:

延迟glht对象的诀窍是认识到它们非常复杂。 有必要首先隔离所需的零件, 然后Hmisc中的latex()函数工作得很好

此示例基于glht中的一个示例

library(Hmisc)
library(multcomp)

### set up a one-way ANOVA
amod <- aov(breaks ~ tension, data = warpbreaks)

### set up all-pair comparisons for factor `tension'
### using a symbolic description (`type' argument
### to `contrMat()')
amod.glht <- glht(amod, linfct = mcp(tension = "Tukey"))

latex(confint(amod.glht)$confint, dec=3)
库(Hmisc)
图书馆(multcomp)
###建立一个单向方差分析

amod您可以在下面的代码示例中找到答案:

multcomp:::print.summary.glht 

x<-glht(...)    
pq<-summary(x)$test

mtests <- cbind(pq$coefficients, pq$sigma, pq$tstat, pq$pvalues)
error <- attr(pq$pvalues, "error")
pname <- switch(x$alternativ, 
                 less = paste("Pr(<", ifelse(x$df ==0, "z", "t"), ")", sep = ""), 
                 greater = paste("Pr(>", ifelse(x$df == 0, "z", "t"), ")", sep = ""), 
                 two.sided = paste("Pr(>|", ifelse(x$df == 0, "z", "t"), "|)", sep = ""))                                                                   
colnames(mtests) <- c("Estimate", "Std. Error", ifelse(x$df ==0, "z value", "t value"), pname)

xtable(mtests)
multcomp:::print.summary.glht

X输出看起来不够规则,无法全部放入表中。你不应该将常规部分与标题和注释分开吗?我不是问如何设置表格的格式,我想问的是,是否有一个库为我提供创建特定对象的代码,并请求帮助查看该对象的打印方法,以查看如何修改该对象以执行您所请求的任何操作,这在我看来是不精确的。仅供参考,交叉发布到R-help是一种糟糕的形式——选择SO或R-help。否则答案和澄清请求会变得混乱。如果您能详细说明一下您的答案,那就太好了。
confint(amod.glht)$confint
是一个普通的矩阵,所以您不需要
Hmisc
xtable
就足够了)。