rm(list = ls())
install.packages(“dplyr”) # jalankan sekali jika belum ada
library(dplyr)
data <- read.csv2(file.choose()) data
set.seed(123) B <- 2000 # jumlah simulasi n1 <- 10 n2 <- 30 mean_n10 <- replicate(B, mean(sample(data\(pendapatan, n1, replace = TRUE))) mean_n30 <- replicate(B, mean(sample(data\)pendapatan, n2, replace = TRUE))) sampling_df <- data.frame(mean_n10 = mean_n10, mean_n30 = mean_n30) sampling_df hist(mean_n10, main=“Distribusi Sampling Mean Pendapatan (n=10)”, xlab=“Mean Sampel”) hist(mean_n30, main=“Distribusi Sampling Mean Pendapatan (n=30)”, xlab=“Mean Sampel”)