微信公众号搜"智元新知"关注
微信扫一扫可直接关注哦!

按组使用 describeBy将描述性统计数据导出到 LaTeX

如何解决按组使用 describeBy将描述性统计数据导出到 LaTeX

我无法将按变量分组的描述性统计表导出到 LaTeX。

我的真实数据中有很多组变量的情况(并且比下面的示例要描述的变量更多),但这里的示例数据应该足以说明我的问题。

假设我有以下数据:

structure(list(election = c("USA_2008","BRA_2014","TWN_2008","BRA_2006","MEX_2012","BRA_2010","PER_2000","CHL_1999","USA_2008","USA_2012","URY_2009","FRA_2012","BRA_2002","USA_1996","PHL_2016","TWN_2012","ROU_2009","PER_2016","MEX_2006","ARG_2015","ROU_2014","ROU_2004","PHL_2010","MEX_2000","TWN_2004","PHL_2004","TWN_1996","PER_2011","CHL_2005","KEN_2013","PER_2001","USA_2004","CHL_2009","PER_2006","USA_2008"),ideol_self = c(2,8,NA,10,5,6,7,3,9,4,2,1,8),ideolparty_A = c(2,0),ideolparty_B = c(8,2.5,8.5,10),ideolparty_C = c(NA,NA
),ideolparty_D = c(NA,NA),ideolparty_E = c(NA,ideolparty_F = c(NA,ideolparty_G = c(NA,ideolparty_H = c(NA,ideolparty_I = c(NA,ideol_leader_A = c(2L,6L,8L,0L,10L,4L,5L,9L,3L,1L,1L),ideol_leader_B = c(8L,7L,2L,ideol_leader_C = c(NA,ideol_leader_D = c(NA,ideol_leader_E = c(NA,ideol_leader_F = c(NA,ideol_leader_G = c(NA,ideol_leader_H = c(NA,ideol_leader_I = c(NA,ideol_Voted_PR_1 = c(2,ideol_elected_PR_1 = c(2,kNowledge_adj = c(1,0.666666666666667,0.333333333333333,NaN,0.75,0.5,0.25,1),abs_growth2 = c(4.6408656404,0.420588274400036,10.8843736,13.3688332656,12.2270156288,12.3130505024,3.81362268409997,11.6114849984,4.6408656404,8.57400049520001,19.567637804,0.0300030000999874,9.0691721405,10.8450427508,22.2297293888,3.95091969999999,6.13615700000001,12.0974731072,12.6843322465,7.27983820530003,9.14544143360003,19.455863534,13.3857564175,13.38966386,16.091099135,17.3167134808,18.0004356286,18.0597050009,4.86624550219998,8.62645503200001,7.4901184544,19.9626013952,4.6408656404),compulsory = c(5,5),closest = c(2,10)),row.names = c(NA,-175L),class = "data.frame")

按变量“选举”分组的汇总统计

library(psych)
sum1 <- describeBy(df,df$election) 
sum1

由此创建的输出正是我正在寻找的。但我无法将其导出到 LaTeX。

首先尝试使用 here 提供的解决方案将 describeBy 生成的对象转换为数据帧,以便我可以运行 stargazer 导出到 LaTeX:

library(stargazer)
sumdf<-do.call("rbind",sum1)
stargazer(sumdf)

然后我尝试了一个 describe 选项,使 outpub 成为一个矩阵,而不是一个列表:

sum1 <- describeBy(cses1,cses1$election,mat = TRUE) 

以上选项均无效。两者都返回 LaTeX 代码,用于完整数据集的摘要或类似内容,但不会像我预期的那样返回按组变量拆分的描述性统计数据。

我也尝试使用其他软件包而不是 describeBy,但也没有成功。也许可以用 ddply 做一些事情(比如通过变量“选举”分割数据集,然后运行一个函数从所有数据集中获取描述性统计数据)。

版权声明:本文内容由互联网用户自发贡献,该文观点与技术仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 dio@foxmail.com 举报,一经查实,本站将立刻删除。