小编gen*_*iva的帖子

我正在尝试编写一个创建动画图形的函数(不使用动画包),用户可以在其中控制演示中心极限定理的输入(样本大小和分布等).理论上这就是我想要的,但是如上所述,在编写用户可以实际控制输入的功能方面遇到了麻烦.

msample <- NA # set up empty vector
ns <-3 # sample size
for(i in 1:500){
sam <- runif(ns) * 10 # draw sample
 msample[i] <- mean(sam) # save mean of sample
h <- hist(msample, breaks=seq(0,10, len=50), # histogram of all means
xlim=c(0,10), col=grey(.9),
xlab="", main="Central Limit Theorem", border="blue", las=1)
points(sam, rep(max(h$count), length(sam)),
pch=16, col=grey(.2)) # add sampled values
points(msample[i], max(h$count), # add sample mean value
col="red", pch=15)
text(10, max(h$count), paste("sample no", i))
hist(msample[i], breaks=seq(0,10, len=50), # ovelay sample …

这就是我所拥有的,但它会返回所有50个州

grep("[a]|[i]|[o]|[u]",state.name, value=TRUE)