hongli32 发表于 2012-3-12 08:39 
另外shapiro.test(sample.mean)里的sample.mean应该是一组数据,你那个写法怎么能出来一组大于3个数的数据呢 ...
sample.mean <- numeric (200)
> for(i in 1:200) sample.mean
<- mean(rexp(100))
> shapiro.test (sample.mean)
Shapiro-Wilk normality test
data: sample.mean
W = 0.994, p-value = 0.601
sample.mean <- numeric (2000)
> for(i in 1:2000) sample.mean<-mean(rexp(100))
> shapiro.test(sample.mean)
Shapiro-Wilk normality test
data: sample.mean
W = 0.9926, p-value = 1.699e-08
这又如何解释呢?