今回、

1班はThalassiosira属の珪藻を暗条件で培養したものを

2班はPhaeodmctylum属の珪藻を暗条件で培養したものを

3班はThalassiosira属の珪藻を明条件で培養したものを

4班はPhaeodmctylum属の珪藻を明条件で培養したものを

測定してもらいました。

クロロフィルの測定時に希釈している場合は、それが反映されていない可能性があるので、教えてください。

1 亜硝酸態窒素 時系列

# コードを一つだけ示す。汚コード。
nitrite_data <- read_excel(path = "./teaching_assistant/data_student.xlsx",
                          sheet = "data") %>% 
  filter(nutrient == "nitrite") %>%
  select(-one_of("f/2"))

nitrite_data <- nitrite_data %>% 
  pivot_longer(cols = 3:6,
               names_to = "date",
               values_to = "Abs") %>% 
  mutate(date = as.Date(date, format = "%Y%m%d")) %>%  
  mutate(concentration　 = 23.357 * Abs + 0.3968) %>% 
  group_by(group, date) %>%
  summarise(.groups = "drop",
            mean_concentration_nitrite = mean(concentration), 
            std_error = sd(concentration) / sqrt(n()))

g_nitrite <- ggplot(data = nitrite_data, 
                    aes(x = date,
                    y = mean_concentration_nitrite,
                    group = group,
                    color = group)) +
  geom_point() +
  geom_errorbar(aes(ymin = mean_concentration_nitrite - std_error, ymax = mean_concentration_nitrite + std_error), width = 0.2) + 
  labs(title = "亜硝酸態窒素",
       x = "日付",
       y = "濃度 (μmol/L)") +
  scale_color_hue(labels = c(group1 = "T 暗", group2 = "P 暗", group3 = "T 明", group4 = "P 明")) +
  theme_bw(base_family = "HiraKakuProN-W3")

g_nitrite

2 硝酸態窒素 時系列

3 アンモニア態窒素 時系列

4 リン酸態リン 時系列

5 溶存ケイ酸 時系列

6 クロロフィル 時系列