ベイズ推論に必要な要素

RによるベイズデータAnalysisの基礎

Rasmus Bååth

Data Scientist

生成モデルとは？

ゾンビ薬の生成モデル

# Parameters
prop_success <- ???
n_zombies <- ???

ゾンビ薬の生成モデル

# Parameters
prop_success <- 0.15
n_zombies <- 13

ゾンビ薬の生成モデル

# Parameters
prop_success <- 0.15
n_zombies <- 13
# Simulating data
data <- c()
for(zombie in 1:n_zombies) {
  data[zombie] <- ???
}

ゾンビ薬の生成モデル

# Parameters
prop_success <- 0.15
n_zombies <- 13
# Simulating data
data <- c()
for(zombie in 1:n_zombies) {
  data[zombie] <- runif(1, min = 0, max = 1) < prop_success
}

ゾンビ薬の生成モデル

# Parameters
prop_success <- 0.15
n_zombies <- 13
# Simulating data
data <- c()
for(zombie in 1:n_zombies) {
  data[zombie] <- runif(1, min = 0, max = 1) < prop_success}

data

FALSE FALSE FALSE TRUE FALSE FALSE FALSE FALSE TRUE FALSE TRUE FALSE FALSE

ゾンビ薬の生成モデル

# Parameters
prop_success <- 0.15
n_zombies <- 13
# Simulating data
data <- c()
for(zombie in 1:n_zombies) {
  data[zombie] <- runif(1, min = 0, max = 1) < prop_success
}
data <- as.numeric(data)

ゾンビ薬の生成モデル

# Parameters
prop_success <- 0.15
n_zombies <- 13
# Simulating data
data <- c()
for(zombie in 1:n_zombies) {
  data[zombie] <- runif(1, min = 0, max = 1) < prop_success
}
data <- as.numeric(data)
data

0 0 0 1 0 0 0 0 1 0 1 0 0

ゾンビ薬の生成モデル

# Parameters
prop_success <- 0.15
n_zombies <- 13
# Simulating data
data <- c()
for(zombie in 1:n_zombies) {
  data[zombie] <- runif(1, min = 0, max = 1) < prop_success
}
data <- as.numeric(data)
data

0 0 1 0 0 0 0 0 0 0 0 0 0

ゾンビ薬の生成モデル

# Parameters
prop_success <- 0.15
n_zombies <- 13
# Simulating data
data <- c()
for(zombie in 1:n_zombies) {
  data[zombie] <- runif(1, min = 0, max = 1) < prop_success
}
data <- as.numeric(data)
data

0 1 0 1 1 0 0 1 0 1 0 0 0

ゾンビ薬の生成モデル

# Parameters
prop_success <- 0.15
n_zombies <- 13
# Simulating data
data <- c()
for(zombie in 1:n_zombies) {
  data[zombie] <- runif(1, min = 0, max = 1) < prop_success
}
data <- as.numeric(data)
data

0 0 0 0 0 0 0 1 0 0 0 0 0

ゾンビ薬の生成モデル

# Parameters
prop_success <- 0.15
n_zombies <- 13
# Simulating data
data <- c()
for(zombie in 1:n_zombies) {
  data[zombie] <- runif(1, min = 0, max = 1) < prop_success
}
data <- as.numeric(data)
data

0 0 0 0 1 0 0 0 0 1 0 1 0

このモデルを試してみましょう！

RによるベイズデータAnalysisの基礎