[R] plot - central limit theorem *solved*

[Jörg Groß]

Thanks for your great help!

I have now what I wanted.
For sure it's not well written and you can realize it with much less  
lines.
But it works and it's like how I wanted it to look like;


z1 <- rnorm(4500, mean=20, sd=5)
z2 <- rnorm(3600,  mean=28, sd=5)
z3 <- rnorm(1300,  mean=40, sd=7)
z4 <- rnorm(1200,  mean=10, sd=7)
z5 <- rnorm(1500,  mean=50, sd=5)
z6 <- rnorm(1500,  mean=10, sd=5)
z7 <- rnorm(1400,  mean=37, sd=3)
pop <- c(z1, z2, z3, z4, z5, z6, z7)

  meansN2 <- numeric()
  meansN10 <- numeric()
  meansN40 <- numeric()

for (i in 1:1000) {
			x <- runif(2, min=1, max=15000)
			x <- round(x)
			sample <- c(pop[x])
			meansN2 <- append(meansN2, mean(sample))
        				 }
for (i in 1:1000) {
			x <- runif(10, min=1, max=15000)
			x <- round(x)

			sample <- c(pop[x])
			meansN10 <- append(meansN10, mean(sample))
				 }
for (i in 1:1000) {
			x <- runif(40, min=1, max=15000)
			x <- round(x)
			sample <- c(pop[x])
			meansN40 <- append(meansN40, mean(sample))
			 	}

pop <- table(cut(pop,breaks=seq(from=-10, to=70, length.out=80)))

meansN2 <- table(cut(meansN2,breaks=seq(from=-10, to=70,  
length.out=80)))
meansN10 <- table(cut(meansN10,breaks=seq(from=-10, to=70,  
length.out=80)))
meansN40 <- table(cut(meansN40,breaks=seq(from=-10, to=70,  
length.out=80)))

pop <- smooth(pop)
meansN2 <- smooth(meansN2)
meansN10 <- smooth(meansN10)
meansN40 <- smooth(meansN40)

plot(pop, type="l", lwd=1.5)
lines(meansN2, lty="dotted", lwd=1.5)
lines(meansN10, lty="dashed", lwd=1.5)
lines(meansN40, lty="twodash", lwd=1.5)