Table1函数制作表格 加上P值
以下参考这个网址https://cran.r-project.org/web/packages/table1/vignettes/table1-examples.html
根据基因高低表达量两组，
rm(list =ls())
library(MatchIt)
setwd("C:\\Users\\17657\\Desktop\\ska3")
library(table1)
library(xlsx)
baseline<- read.xlsx('卡方检验.xlsx',sheetIndex = 1,header = T,encoding = 'UTF-8')
baseline$SKA3<- factor(baseline$SKA3, levels=c("High","Low","P-value"), labels=c("High","Low","P-value"))
rndr <- function(x, name, ...) {
if (length(x) == 0) {
y <- baseline[[name]]
s <- rep("", length(render.default(x=y, name=name, ...)))
if (is.numeric(y)) {
p <- t.test(y ~ baseline$SKA3)$p.value
} else {
p <- chisq.test(table(y, droplevels(baseline$SKA3)))$p.value
}
s[2] <- sub("<", "<", format.pval(p, digits=3, eps=0.001))
s
} else {
render.default(x=x, name=name, ...)
}
}
rndr.strat <- function(label, n, ...) {
ifelse(n==0, label, render.strat.default(label, n, ...))
}
table1(~ age +T +M +N +stage +ER +PR +HER2+subtype+fustat| SKA3,
data=baseline, droplevels=F, render=rndr, render.strat=rndr.strat)
Missing值不知道如何去除，要注意是否用fisher检验
