Rで確率分布:ロジスティック分布

# 必要なライブラリを読み込みます
library(ggplot2)
library(dplyr)
library(tidyr)

# 1. 描画範囲となるx軸の値を生成
x_vals <- seq(-10, 15, length.out = 1000)

# 2. 異なるパラメータを持つ分布の確率密度を計算
# ロジスティック分布は dlogis(x, location, scale) を使用
# 標準ロジスティック(μ=0, s=1)の分散は (pi^2 * 1^2)/3
logis_var <- (pi^2 * 1^2) / 3
norm_sd <- sqrt(logis_var) # 対応する正規分布の標準偏差

df <- tibble(
  x = x_vals
) %>%
  mutate(
    `ロジスティック (μ=0, s=1)` = dlogis(x, location = 0, scale = 1),
    `ロジスティック (μ=5, s=1)` = dlogis(x, location = 5, scale = 1),
    `ロジスティック (μ=0, s=2)` = dlogis(x, location = 0, scale = 2),
    `正規分布 (μ=0, σ≈1.8)` = dnorm(x, mean = 0, sd = norm_sd)
  )

# 3. ggplotで描画しやすいように、データを「ロングフォーマット」に変換
df_long <- df %>%
  pivot_longer(
    cols = -x,
    names_to = "distribution",
    values_to = "density"
  ) %>%
  # 凡例の順序を調整
  mutate(distribution = factor(distribution, levels = c(
    "ロジスティック (μ=0, s=1)",
    "ロジスティック (μ=5, s=1)",
    "ロジスティック (μ=0, s=2)",
    "正規分布 (μ=0, σ≈1.8)"
  )))

# 4. 各分布に割り当てる色と線種を定義
manual_colors <- c(
  `ロジスティック (μ=0, s=1)` = "blue",
  `ロジスティック (μ=5, s=1)` = "red",
  `ロジスティック (μ=0, s=2)` = "darkgreen",
  `正規分布 (μ=0, σ≈1.8)` = "black"
)
manual_linetypes <- c(
  `ロジスティック (μ=0, s=1)` = "solid",
  `ロジスティック (μ=5, s=1)` = "solid",
  `ロジスティック (μ=0, s=2)` = "solid",
  `正規分布 (μ=0, σ≈1.8)` = "dashed"
)

# 5. ggplotを使用してチャートを描画
p <- ggplot(df_long, aes(x = x, y = density, color = distribution, linetype = distribution)) +
  geom_line(linewidth = 1.1) +
  scale_color_manual(values = manual_colors) +
  scale_linetype_manual(values = manual_linetypes) +
  labs(
    title = "ロジスティック分布と正規分布の比較",
    subtitle = "形状は似ているが、ロジスティック分布の方がわずかに裾が重い",
    x = "xの値",
    y = "確率密度",
    color = "分布の種類",
    linetype = "分布の種類"
  ) +
  theme_minimal(base_size = 14) +
  theme(
    legend.position = "top",
    plot.title = element_text(face = "bold"),
    plot.subtitle = element_text(color = "gray40")
  ) +
  guides(color = guide_legend(nrow = 2))

# チャートの表示
print(p)