diff --git a/R/make_boxplot.R b/R/make_boxplot.R
new file mode 100644
index 0000000000000000000000000000000000000000..8613c5673ade605595a9950ddea5b7d9b9a02947
--- /dev/null
+++ b/R/make_boxplot.R
@@ -0,0 +1,111 @@
+library(tidyverse)
+
+#繝代せ謌仙粥邇�
+# CSV繝輔ぃ繧､繝ｫ縺ｮ隱ｭ縺ｿ霎ｼ縺ｿ
+df_epv <- read_csv("~/Desktop/logs/epv2base30half/pass_prob.csv") |>
+  select(pass_success_rate) |>
+  rename(epv = pass_success_rate)
+
+df_base <- read_csv("~/Desktop/logs/base2base30half/pass_prob.csv") |>
+  select(pass_success_rate) |>
+  rename(base = pass_success_rate)
+
+# 繝��繧ｿ繝輔Ξ繝ｼ繝�縺ｮ邨仙粋
+df <- bind_cols(df_epv, df_base)
+
+long_data <- df %>%
+  gather(key = "Category", value = "Value", epv, base)
+
+p <- long_data %>%
+    ggplot(aes(x = Category, y = Value)) +
+    geom_boxplot() +
+    scale_x_discrete(labels = c("HELIOS Base", "EPV繧ｨ繝ｼ繧ｸ繧ｧ繝ｳ繝�")) + 
+    scale_y_continuous(limits = c(0.65, 1)) + 
+    stat_summary(fun = mean, geom = "point", shape = 20, size = 3, color = "red") +
+    xlab("繧ｨ繝ｼ繧ｸ繧ｧ繝ｳ繝�") +
+    ylab("繝代せ謌仙粥邇�")
+
+ggsave("prob.pdf", p, device = cairo_pdf, width = 3, height = 3)
+
+
+# 繝懊�繝ｫ菫晄戟邇�
+# CSV繝輔ぃ繧､繝ｫ縺ｮ隱ｭ縺ｿ霎ｼ縺ｿ
+df_epv <- read_csv("~/Desktop/logs/epv2base30half/EPV.csv") |>
+  select(our_possession) |>
+  rename(epv = our_possession)
+
+df_base <- read_csv("~/Desktop/logs/base2base30half/HELIOS_base.csv") |>
+  select(our_possession) |>
+  rename(base = our_possession)
+
+# 繝��繧ｿ繝輔Ξ繝ｼ繝�縺ｮ邨仙粋
+df <- bind_cols(df_epv, df_base) 
+
+long_data <- df %>%
+  gather(key = "Category", value = "Value", epv, base)
+
+p <- long_data %>%
+    ggplot(aes(x = Category, y = Value)) +
+    geom_boxplot() +
+    scale_x_discrete(labels = c("HELIOS Base", "EPV繧ｨ繝ｼ繧ｸ繧ｧ繝ｳ繝�")) + 
+    scale_y_continuous(limits = c(0.3, 0.9)) + 
+    stat_summary(fun = mean, geom = "point", shape = 20, size = 3, color = "red") +
+    xlab("繧ｨ繝ｼ繧ｸ繧ｧ繝ｳ繝�") +
+    ylab("繝懊�繝ｫ菫晄戟邇�")
+
+ggsave("possession.pdf", p, device = cairo_pdf, width = 3, height = 3)
+
+# 蠕礼せ
+# CSV繝輔ぃ繧､繝ｫ縺ｮ隱ｭ縺ｿ霎ｼ縺ｿ
+df_epv <- read_csv("~/Desktop/logs/epv2base30half/EPV.csv") |>
+  select(our_point) |>
+  rename(epv = our_point)
+
+df_base <- read_csv("~/Desktop/logs/base2base30half/HELIOS_base.csv") |>
+  select(our_point) |>
+  rename(base = our_point)
+
+# 繝��繧ｿ繝輔Ξ繝ｼ繝�縺ｮ邨仙粋
+df <- bind_cols(df_epv, df_base)
+
+long_data <- df %>%
+  gather(key = "Category", value = "Value", epv, base)
+
+p <- long_data %>%
+    ggplot(aes(x = Category, y = Value)) +
+    geom_boxplot() +
+    scale_x_discrete(labels = c("HELIOS Base", "EPV繧ｨ繝ｼ繧ｸ繧ｧ繝ｳ繝�")) + 
+    scale_y_continuous(limits = c(0, 4.5)) + 
+    stat_summary(fun = mean, geom = "point", shape = 20, size = 3, color = "red") +
+    xlab("繧ｨ繝ｼ繧ｸ繧ｧ繝ｳ繝�") +
+    ylab("蠕礼せ")
+
+ggsave("point.pdf", p, device = cairo_pdf, width = 3, height = 3)
+
+
+# 螟ｱ轤ｹ
+# CSV繝輔ぃ繧､繝ｫ縺ｮ隱ｭ縺ｿ霎ｼ縺ｿ
+df_epv <- read_csv("~/Desktop/logs/epv2base30half/EPV.csv") |>
+  select(opp_point) |>
+  rename(epv = opp_point)
+
+df_base <- read_csv("~/Desktop/logs/base2base30half/HELIOS_base.csv") |>
+  select(opp_point) |>
+  rename(base = opp_point)
+
+# 繝��繧ｿ繝輔Ξ繝ｼ繝�縺ｮ邨仙粋
+df <- bind_cols(df_epv, df_base)
+
+long_data <- df %>%
+  gather(key = "Category", value = "Value", epv, base)
+
+p <- long_data %>%
+    ggplot(aes(x = Category, y = Value)) +
+    geom_boxplot() +
+    scale_x_discrete(labels = c("HELIOS Base", "EPV繧ｨ繝ｼ繧ｸ繧ｧ繝ｳ繝�")) + 
+    scale_y_continuous(limits = c(0, 4.5)) + 
+    stat_summary(fun = mean, geom = "point", shape = 20, size = 3, color = "red") +
+    xlab("繧ｨ繝ｼ繧ｸ繧ｧ繝ｳ繝�") +
+    ylab("螟ｱ轤ｹ")
+
+ggsave("opp_point.pdf", p, device = cairo_pdf, width = 3, height = 3)