檢測數據是否符合正態分布（R語言）

norm.test<- function(input.data,alpha=0.05,pic=TRUE){
    if(pic==TRUE){#畫圖形
           dev.new()
           par(mfrow=c(2,1))
           qqnorm(input.data,main="qq圖")
           qqline(input.data)
           hist(input.data,frep=F,main="直方圖和密度估計曲線")
           lines(density(input.data),col="blue") #密度估計曲線
           x<- c(round(min(input.data)):round(max(input.data)))
           lines(x,dnorm(x,mean(input.data),sd(input.data)),col="red") #正態分布曲線
    }
    sol<- shapiro.test(input.data)
    if(sol$p.value>alpha){
        print(paste("success:服從正態分布,p.value=",sol$p.value,">",alpha))    
    }else{
         print(paste("error:不服從正態分布,p.value=",sol$p.value,"<=",alpha))
    }
    sol
}