有数百列的数据想分析每一列数据的残差正态性，可是手动输入基本不可能实现，因为数据太多。想用for循环语句进行每一列的数据分析，也在网上找了一些代码，但根本做不出来，请教下各位老师。
**#p_otu<-read.csv("门相对丰度.csv", row.names=1, nrows = 32)
N1=c("Phylum")
P=c("残差正态性")
data<-p_otu[,1:135]
for(r in 1:ncol(data)){
N3=colnames(data)[r]
P1=shapiro.test(lm(p_otu[,r]~OTU,data=p_otu)$residuals)
N1=c(N1,N3)
P=c(P,P1)
}
Q=data.frame(N1,P)