A1、B1、C1 都是實驗 1 的重復。A2、B2、C2 是實驗 2 的重復，依此類推。我一直試圖在相應的時間點對三個點進行 SD 以創建誤差線。我處理過的代碼。我試圖變異、總結和收集，但我一定做錯了什么。
到目前為止的代碼：
Data2 <- select(Data1, Time = c(1), x= c(2:4), y = c(5:7), a = c(8:10), b = c(11:13))
SD =
Data3 <- Data2 %>%
mutate(SD = sd(unlist(select(cur_data(), x1:x3))))
gather(key = "test", value = "value", c(-Time)) %>%
group_by(Data2, Time, c(2:4), c(5:7), c(8:10), c(11:13)) %>%
summarise(sd = sd(value, na.rm=F))
ggplot(aes(x = Time, y = value, colour= test))
geom_point(aes(col = test), alpha = 1)
geom_smooth(aes(col = test, method = "lm"))
geom_errorbar(aes(ymin=value-sd, ymax=value sd), width=5, size = 0.5)
facet_wrap(.~test, ncol = 2)
theme_classic()
xlab("Time (hours)")
ylab("Weight (grams)")
legend(=null)
就像我說的那樣，什么都沒奏效，過去 3 天我一直在試驗。不知道我哪里出錯了，或者在這一點上是什么。
group_by(Time) %>%
summarise(sd= sd(value, na.rm=T))
這段代碼最接近作業，但 SD 我所有的實驗，而不是三個重復。
提前致謝。
資料：
structure(list(Time = c(0, 24, 48, 72, 96, 120, 142), A1 = c(7.687,
NA, NA, 7.687, 7.738, 7.796, 7.747), B1 = c(7.661, 7.68, 7.694,
7.728, 7.707, NA, NA), C1 = c(7.684, 7.991, 8.027, 7.978, 7.966,
NA, NA), A2 = c(7.671, NA, NA, 7.684, 7.679, 7.731, 7.69), B2 = c(7.672,
7.85, 7.838, 7.835, 7.769, NA, NA), C2 = c(7.67, 7.906, 8.02,
8.032, 8.04, NA, NA), A3 = c(7.678, NA, NA, 7.794, 7.731, 7.817,
7.862), B3 = c(7.689, 7.882, 7.937, 7.937, 7.936, NA, NA), C3 = c(7.68,
7.975, 8.003, 7.971, 7.971, NA, NA), A4 = c(7.703, NA, NA, 7.714,
7.752, 7.737, 7.793), B4 = c(7.703, 7.916, 7.992, 7.991, 7.964,
NA, NA), C4 = c(7.706, 7.839, 8.384, 8.391, 8.397, NA, NA)), row.names = c(NA,
-7L), class = c("tbl_df", "tbl", "data.frame"))
uj5u.com熱心網友回復：
不確定您期望的輸出（可能在您的問題中包含一個小草圖？）
這是第一次嘗試/方法
library(tidyverse)
Data1 %>%
pivot_longer(
cols = -1,
names_to = c("try", "exp"),
names_pattern = "(.)(.)") %>%
ggplot(aes(x = Time, y = value))
geom_point()
stat_summary(func = "sd",
geom = "errorbar")
facet_grid(~exp)
轉載請註明出處，本文鏈接：https://www.uj5u.com/shujuku/453826.html
上一篇：ggsurvplot(Survminer)自定義標題（添加斜體）
下一篇：格式不適用于箱線圖中的奇異點