#-----------------------数据导入-----------------------------
setwd("C:\\Users\\Administrator\\Desktop\\脑卒中论文")
data<-read.csv("mytable.csv")

#-----------------------转化数据格式-------------------------
str(data)

data[,c(5:23)]<-lapply(data[,c(5:23)],as.factor)
str(data)
summary(data)
#注意：结局变量不要转化为因子型


table(data$zaifa)#分类数据


#-------------------------1.基线信息表---------------------------
#连续变量报告为【正态】平均值±标准差（SD）或【非正态分布】中位数四分位距（IQR）
#分类变量表示为例数(百分比)

##----正态性检验----

#SW检验(样本量<5000 P<0.05是非正态)

lapply(data[,c(24:76)],shapiro.test)

##----绘制基线信息表----
library(tableone)
tab1<-CreateTableOne(
  vars = colnames(data[,c(5:76)]),#描述的所有变量
  factorVars = colnames(data[,c(5:23)]),#分类变量
  strata = "zaifa",#分组
  data = data,
  addOverall = TRUE
)
tab1

tab<-print(tab1,
           nonnormal = colnames(data[,c(24:76)]),##非正态对应的变量
           noSpaces = T,##是否删除为对齐增加的空格
           showAllLevels = TRUE)#显示所有水平

write.csv(tab,"C:\\Users\\Administrator\\Desktop\\脑卒中论文\\脑卒中再发基线信息表.csv")