截断的正太分布从某种意义已经不服从正太分布了。我写了一个算法得到的数据是正太分布。
a<-rnorm(1000,100,100)
while(length(which(a<0))>0){
a[which(a<0)]<-rnorm(length(which(a<0)),100,100)
}
shapiro.test(a)
思路:
首先生存1000个服从u=100,sd=100的正态分布数列
然后找到那些小于0的,重新生成数目一样的正态分布数列。
用while重复这个过程,直到没有小于0的
最后用shapiro.test()检验通过。