以正态检验为例
设置工作环境后
读取数据:
A <- read.csv("data.csv")
加载包
library(nortest)
查看变量
names(A)
设置一个表格
shapt1 <- c()
for(i in c(3,9:12,13:16,21:22)){#循环
fit <- ad.test(B[,i])#检验
P <- fit$p.value#读值
Met1 <- fit$method#检验方法
result1<-c()#设置内部表格
result1<-data.frame(rbind(result1,cbind(P,Met1)))
result1$variate.name <- names(B)[i]
shapt1<-rbind(shapt1,result1)
}
shapt1#输出