14 lines
502 B
R
14 lines
502 B
R
nooutput=sapply(1:3,function(i){
|
|
n=10
|
|
mu=10
|
|
sigma=4
|
|
sem=sigma/sqrt(n)
|
|
|
|
# Compute 1000 times the mean a random sample and get the standard deviation of the mean
|
|
means_sd=sd(sapply(1:1000,function(dummy){mean(rnorm(n,mean=mu,sd=sigma))}))
|
|
|
|
# See how SEM is a good approximation of the mean standard deviation
|
|
message(paste("----- Experiment",i,"-----"))
|
|
message(paste0("Means SD: ", round(means_sd,digits=2)))
|
|
message(paste0("SEM ", round(sem,digits=2)))
|
|
})
|