我和R.
考虑到我生成的下列随机数据,我能够使用这些数据创建一个应急表:
library(memisc)
library(dplyr)
set.seed(123)
v1 <- c("2010-2011","2011-2012", "2012-2013", "2013-2014", "2014-2015")
v2 <- c("A", "B", "C", "D", "E")
v3 <- c("Z", "Y", "X", "W" )
v4 <- c("data_1", "data_2", "data_3", "data_4" )
dates <- as.factor(sample(v1, 1000, replace=TRUE, prob=c(0.5, 0.2, 0.1, 0.1, 0.1)))
types <- as.factor(sample(v2,1000, replace=TRUE, prob=c(0.3, 0.2, 0.1, 0.1, 0.1)))
types2 <- as.factor(sample(v3, 1000, replace=TRUE, prob=c(0.3, 0.5, 0.1, 0.1)))
names <- as.factor(sample(v3, 1000, replace=TRUE, prob=c(0.3, 0.5, 0.1, 0.1)))
var = rnorm(1000,10,10)
problem_data = data.frame(var,dates, types, types2, names)
summary <- xtabs(~dates+names+types+types2, problem_data)
t = ftable(summary, row.vars=1, col.vars=2:4)
show_html(t)

可以直接从数据帧生成类似于表的应急表吗?
例如,假设我想创建上面的列列表,,但不是用“计数”填充这个表,而是用"var"的平均值填充该表。使用"dplyr“库,我可以创建一个数据框架,其中包含这个应急表所需的所有值:
library(dplyr)
contingency_table = data.frame(problem_data %>% group_by(dates,names, types, types2) %>% summarise(mean_value = mean(var)))
head(contingency_table)
dates names types types2 mean_value
1 2010-2011 W A X -10.128687
2 2010-2011 W A Y 9.552724
3 2010-2011 W A Z 9.686354
4 2010-2011 W B W -4.411400
5 2010-2011 W B Y 13.624970
6 2010-2011 W B Z 7.008089能否将上面的数据帧制作成一个应急表,然后转换成html可发布的格式?
使用这里的堆栈溢出帖子(是否有一种(容易)的方法将平面应急表(ftable)转换为可挠性的),我尝试使用所提供的函数将这个数据帧转换为一个应急表--但它并没有给出所需的结果(也就是说,它与上面的表不一样):
ftable_to_flextable <- function( x ){
row.vars = attr( x, "row.vars" )
col.vars = attr( x, "col.vars" )
rows <- rev( expand.grid( rev(row.vars), stringsAsFactors = FALSE ) )
cols <- rev(expand.grid( rev(col.vars), stringsAsFactors = FALSE ))
xmat <- as.matrix(x)
cols$col_keys = dimnames(xmat)[[2]]
xdata <- cbind(
data.frame(rows, stringsAsFactors = FALSE),
data.frame(xmat, stringsAsFactors = FALSE)
)
names(xdata) <- c(names(row.vars), cols$col_keys)
ft <- regulartable(xdata)
ft <- set_header_df(ft, cols)
ft <- theme_booktabs(ft)
ft <- merge_v(ft, j = names(row.vars))
ft
}
library(flextable)
library(magrittr)
ftable(contingency_table, row.vars = 1:2, col.vars = 3:4) %>% ftable_to_flextable()

是否可以制作一个应急表,而不是使用每个组的变量"var“的平均值来代替计数?使用xtabs()函数做这件事更好,而且我使用“数据框架方法”对此过于复杂了吗?有人能告诉我怎么做吗?
谢谢!
发布于 2022-01-12 04:02:25
通用交叉制表方法在R中称为tapply,您可以阅读更详细的这里。例如,你可以
tapply(problem_data$var, problem_data[, -1L], mean)这给出了类似于xtabs(~dates + types + types2 + names, data = problem_data)的东西,只不过单元格中的所有值现在都是组方法。你可以做这样的事
smy <- tapply(problem_data$var, problem_data[, -1L], mean)
x <- ftable(smy, row.vars = c(4L, 1L), col.vars = 2:3) # var 1 is the first var you see in names(problem_data[, -1L])
y <- ftable(smy, row.vars = 1L, col.vars = c(4L, 2:3))而memisc::show_html(x)给出

memisc::show_html(y)给出

https://stackoverflow.com/questions/70675871
复制相似问题