#Data Populasi
x <- c(1,2,3,4,5,6,7,8)
x
## [1] 1 2 3 4 5 6 7 8
#Distribusi Probabilitas dari Data Populasi
table(x)/length(x)
## x
## 1 2 3 4 5 6 7 8
## 0.125 0.125 0.125 0.125 0.125 0.125 0.125 0.125
#Rata-rata populasi
mean(x)
## [1] 4.5
#Standar deviasi data populasi
sd(x)*sqrt((length(x)-1)/length(x))
## [1] 2.291288
#################
#Kombinasi sampel
sampel <- combn(x, 2)
sampel
## [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [,11] [,12] [,13] [,14]
## [1,] 1 1 1 1 1 1 1 2 2 2 2 2 2 3
## [2,] 2 3 4 5 6 7 8 3 4 5 6 7 8 4
## [,15] [,16] [,17] [,18] [,19] [,20] [,21] [,22] [,23] [,24] [,25] [,26]
## [1,] 3 3 3 3 4 4 4 4 5 5 5 6
## [2,] 5 6 7 8 5 6 7 8 6 7 8 7
## [,27] [,28]
## [1,] 6 7
## [2,] 8 8
#Rata-rata sampel
ratarata <- colMeans(sampel)
ratarata
## [1] 1.5 2.0 2.5 3.0 3.5 4.0 4.5 2.5 3.0 3.5 4.0 4.5 5.0 3.5 4.0 4.5 5.0 5.5 4.5
## [20] 5.0 5.5 6.0 5.5 6.0 6.5 6.5 7.0 7.5
#Ubah format desimal
ratarata <- format(ratarata, digits = 4)
ratarata
## [1] "1.5" "2.0" "2.5" "3.0" "3.5" "4.0" "4.5" "2.5" "3.0" "3.5" "4.0" "4.5"
## [13] "5.0" "3.5" "4.0" "4.5" "5.0" "5.5" "4.5" "5.0" "5.5" "6.0" "5.5" "6.0"
## [25] "6.5" "6.5" "7.0" "7.5"
table(ratarata)
## ratarata
## 1.5 2.0 2.5 3.0 3.5 4.0 4.5 5.0 5.5 6.0 6.5 7.0 7.5
## 1 1 2 2 3 3 4 3 3 2 2 1 1
table(ratarata)/length(ratarata)
## ratarata
## 1.5 2.0 2.5 3.0 3.5 4.0 4.5
## 0.03571429 0.03571429 0.07142857 0.07142857 0.10714286 0.10714286 0.14285714
## 5.0 5.5 6.0 6.5 7.0 7.5
## 0.10714286 0.10714286 0.07142857 0.07142857 0.03571429 0.03571429
prob <- prop.table(table(ratarata))
barplot(prob)
plot(prob)
#Rata-rata Distribusi Sampling Rata-rata Sampel
ratarata
## [1] "1.5" "2.0" "2.5" "3.0" "3.5" "4.0" "4.5" "2.5" "3.0" "3.5" "4.0" "4.5"
## [13] "5.0" "3.5" "4.0" "4.5" "5.0" "5.5" "4.5" "5.0" "5.5" "6.0" "5.5" "6.0"
## [25] "6.5" "6.5" "7.0" "7.5"
mode(ratarata)
## [1] "character"
ratarata <- as.numeric(ratarata)
mode(ratarata)
## [1] "numeric"
mean(ratarata)
## [1] 4.5
sd(ratarata)*sqrt((length(ratarata)-1)/(length(ratarata)))
## [1] 1.5