tips

## クロス表をcaseのデータに

library(vcdExtra)

Titanic

Titanic_freq <- data.frame(Titanic)

Titanic_case <- expand.dft(Titanic_freq, freq="Freq")

Titanic_case

##ブートストラップの基礎##

set.seed(1234) #乱数を固定

x <- rnorm(1:100,50,10) #平均50，標準偏差10の正規分布に従う乱数を発生

B <- numeric(1000) #1000個結果の入る入れ物を用意

for (i in 1:length(B)) {

b <- sample(x, length(x), replace=TRUE) #xから重複を許して，xと同じ数の標本サイズ（100）でリサンプリング

B[i] <- mean(b) #リサンプリングした標本の平均を保存

}

#パラメトリック

mean(x)

sd(x)

se <- sd(x)/sqrt(length(x))

c(mean(x)-qt(0.975,length(x)-1)*se, mean(x)+qt(0.975,length(x)-1)*se)

#ブートストラップ

mean(B) #平均値

quantile(B,p=c(0.025,0.975)) #95%信頼区間