light <- read.table("C:\\Users\\park\\therbook\\light.txt", header = T)
light
qqnorm(light$speed)
qqline(light$speed, col = "green")
length(light$speed)
a <- numeric(10000)
for (i in 1:10000) a[i] <- mean(sample(light$speed, size = 20, replace = TRUE))
mean(sample(light$speed, size = 20, replace = TRUE))
hist(a, col = "Blue")
summary(a)