修改描述性参数
table1(strata, labels, groupspan=c(1, 3, 1),
render.continuous=c(.="Mean (CV%)",
.="Median [Min, Max]",
"Geo. mean (Geo. CV%)"="GMEAN (GCV%)"))
另外还可以根据需要修改为N, NMISS, MEAN, SD, CV, GMEAN, GCV, MEDIAN, MIN, MAX, IQR, Q1, Q2, Q3, T1, T2, FREQ, PCT等参数。
改变表的外观
使用内置的样式
该package包括有限数量的内置样式，包括：
zebra: 交替的阴影行和非阴影行(斑马条纹)
grid: 显示所有网格线
shade: 将标题行涂成灰色
times: 使用serif字体
center: 将所有列居中，包括第一列（行标签列）
table1(~ age + sex + wt | treat, data=dat, topclass="Rtable1-zebra")
table1(~ age + sex + wt | treat, data=dat, topclass="Rtable1-grid")
table1(~ age + sex + wt | treat, data=dat, topclass="Rtable1-grid Rtable1-shade Rtable1-times")
添加p值
在本例中，分类变量使用独立卡方检验，连续变量使用t检验(如果需要，还可以使用其他检验)。
pvalue <- function(x, ...) { # Construct vectors of data y, and groups (strata) g y <- unlist(x) g <- factor(rep(1:length(x), times=sapply(x, length))) if (is.numeric(y)) { # For numeric variables, perform a standard 2-sample t-test p <- t.test(y ~ g)$p.value } else { # For categorical variables, perform a chi-squared test of independence p <- chisq.test(table(y, g))$p.value } # Format the p-value, using an HTML entity for the less-than sign. # The initial empty string places the output on the line below the variable label. c("", sub("<", "<", format.pval(p, digits=3, eps=0.001))) }
table1(~ age + sex + wt | treat, data=dat, overall=F, extra.col=list(`P-value`=pvalue))
Reference
Using the table1 Package to Create HTML Tables of Descriptive Statistics (r-project.org)
