library(readr)

# データの読み込み
df <- read_csv("yamagata.csv", locale = locale(encoding = "SHIFT-JIS"))
df$`/項目` <- NULL
df$調査年 <- as.numeric(sapply(df$調査年, gsub, pattern="年度", replacement = ""))
df <- df[order(df$調査年),]

# 2015年度の人口増減率を母集団とする
data <- df[df$調査年 == 2015,][["#A05101_人口増減率【％】"]]

print("母集団", quote = FALSE)
print(sort(data))

# 標本サイズ
n <- 20

print("非復元抽出", quote = FALSE)
sdata1 <- sample(data, n)
print(sort(sdata1), quote = FALSE)

print("復元抽出", quote = FALSE)
sdata2 <- sample(data, n, replace = TRUE)
print(sort(sdata2), quote = FALSE)