统计--正态性检验、样本sd与峰度的关系

结论：无关

rm(list=ls())
#样本sd与峰度无关
i=1 
p   <- seq(1,100,1) 
b   <- seq(1,100,1)
while (i<100) {
  a   <- rnorm(100,10,i)
  out1<-basicStats(a)
  
  b[i]   <- out1[15,]
  
  i=i+1
}
resu <- as.data.frame(cbind(p,b))
plot(resu$p,resu$b,ylim=c(-1,1))
write.csv(resu,'resu.csv')

#正态P与峰度无关
i=1 
p   <- seq(1,100,1) 
b   <- seq(1,100,1)
while (i<100) {
  a   <- rnorm(100,10,800)
  out <- shapiro.test(a)
  
  p[i]   <- out$p.value
  out1<-basicStats(a)
  
  b[i]   <- out1[15,]
  
  i=i+1
}
res <- cbind(p,b)
res <- as.data.frame(res)
plot(res$p,res$b,ylim=c(-1,1),xlim=c(0,1))