如何在r中的雙向y軸圖中繪制標準偏差？

我的資料如下所示：

dput(NDVI)
structure(list(date = structure(c(19091, 19091, 19092, 19092, 
19093, 19093, 19094, 19094, 19095, 19095, 19096, 19096, 19097, 
19097, 19098, 19098, 19099, 19099, 19100, 19100, 19101, 19101, 
19102, 19102, 19103, 19103, 19104, 19104, 19105, 19105, 19106, 
19106, 19107, 19107, 19109, 19109, 19110, 19110, 19111, 19111, 
19112, 19112, 19113, 19113, 19114, 19114), class = "Date"), target = c("grass", 
"tree", "grass", "tree", "grass", "tree", "grass", "tree", "grass", 
"tree", "grass", "tree", "grass", "tree", "grass", "tree", "grass", 
"tree", "grass", "tree", "grass", "tree", "grass", "tree", "grass", 
"tree", "grass", "tree", "grass", "tree", "grass", "tree", "grass", 
"tree", "grass", "tree", "grass", "tree", "grass", "tree", "grass", 
"tree", "grass", "tree", "grass", "tree"), cNDVI = c(0.500272817021574, 
0.696083797270416, 0.489691804085761, 0.692795403987576, 0.570071858811957, 
0.698997926130435, 0.510441047568798, 0.693961708101642, 0.497010486993716, 
0.677617144059372, 0.489045582825731, 0.672753166123319, 0.488193790481098, 
0.663609030602538, 0.488885606064579, 0.644550512172325, 0.49094954276713, 
0.625607907143415, 0.477036686933829, 0.606892335456545, 0.46818825238034, 
0.598763281698822, 0.464516342335005, 0.594977691269838, 0.463896574547792, 
0.611135472706338, 0.60754090771395, 0.657544930586677, 0.483871290250419, 
0.641341913760879, 0.455301838422566, 0.625854396853959, 0.461233936815514, 
0.620543555588439, 0.467959523995231, 0.615756464377253, 0.441352232545973, 
0.597588129996248, 0.443084476959249, 0.588500193948734, 0.447515720605795, 
0.569604003059877, 0.445609012507936, 0.562613166422831, 0.495220373276007, 
0.58718847627638), TempAmb_Avg = c(13.16, NA, 7.929, NA, 12.29, 
NA, 10.37, NA, 10.91, NA, 10.14, NA, 9.15, NA, 11.25, NA, 9.17, 
NA, 11.94, NA, 11.26, NA, 9.45, NA, 9.09, NA, NA, NA, 6.447, 
NA, 9.14, NA, 8.02, NA, 10.54, NA, 10.12, NA, 11.56, NA, 12.3, 
NA, 10.82, NA, 11.17, NA), sd_NDVI = c(0.00206275442936463, 0.00483529298045758, 
0.00143646028945461, 0.00547488157335041, 0.0107421213523848, 
NA, 0.0057345829257413, 0.00239673774468361, 0.0122234699109748, 
0.0041293261382514, 0.00159520722600412, 0.00400695247999584, 
0.0014958900631327, 0.00207345410768067, 0.00142529228044798, 
0.00188141954676642, 0.00254072590985734, 0.00374571468034617, 
0.00162042982764391, 0.00371980695360771, 0.00149205232583476, 
0.00860693732895818, 0.00195420993834865, 0.0135370828066329, 
0.00267586232456702, 0.00366685362066545, 0.00548912106770043, 
0.0115656879529442, 0.00892427653015911, 0.00876020639322138, 
0.00793602431786174, 0.00308048721534673, 0.0022395875936781, 
0.00464228152203431, 0.010403456112334, 0.01197424728207, 0.0203370856390955, 
0.00842098123326728, 0.00179785152483323, 0.00238532744345886, 
0.00174703441109131, 0.0016114072840488, 0.00150060553756978, 
0.00245323735514388, 0.0276795273807586, 0.0187674033605546)), row.names = c(NA, 
46L), class = "data.frame")

我正在嘗試在雙向 y 軸上繪制cNDVI 上的標準差sd_NDVI 。

plot <- ggplot(NDVI) 
  geom_bar(aes(x=date, y=TempAmb_Avg),stat="identity",colour="blue", fill="blue") 
  geom_point(aes(x=date, y=30*cNDVI,color=target),stat="identity") 
  geom_line(aes(x=date, y=30*cNDVI,color=target),stat="identity") 
  geom_errorbar(aes(ymin=30*cNDVI-30*sd_NDVI, ymax=30*cNDVI 30*sd_NDVI, width=.2, position = position_dodge(0.05)),color=target, stat="identity") 
  labs(x="TIME",y="ToC")   theme_bw()
plot

我收到此錯誤：

Error in layer(data = data, mapping = mapping, stat = stat, geom = GeomErrorbar,  : 
  object 'target' not found

盡管顯示了繪圖，但 cNDVI 中的標準偏差sd_NDVI并未顯示。如何顯示這些條件下的標準偏差？

先感謝您！

uj5u.com熱心網友回復：

從技術上講，這是您的代碼，其中包含更正的語法錯誤。但是，我想引起您注意的是，tree類別幾乎沒有資料，也沒有明顯需要為 繪制點和線30*sd_NDVI。過濾資料集并洗掉未使用的 geom 可能會更好。您還需要調整圖例。

您可以根據需要調整所需的輸出。

NDVI %>% 
  mutate(target = as.factor(target)) %>% 
ggplot() 
  geom_bar(aes(x=date, y=TempAmb_Avg, fill = target, color = target),
           stat="identity") 
  geom_point(aes(x=date, y=30*sd_NDVI, color=target),
             stat="identity") 
  geom_line(aes(x=date, y=30*sd_NDVI, color=target),
            stat="identity") 
  geom_errorbar(aes(x = date, 
                    y = TempAmb_Avg,
                    ymin=TempAmb_Avg-30*sd_NDVI, 
                    ymax=TempAmb_Avg 30*sd_NDVI, 
                    width=.2,,
                    color = target
                    ),# position = position_dodge(0.05), 
                stat="identity")  
  scale_color_manual(values = c(adjustcolor("blue", 0.5), adjustcolor("red", 0.5)), breaks = c("grass", "tree"))  
  scale_fill_manual(values = c(adjustcolor("blue", 0.3), adjustcolor("red", 0.3)), breaks = c("grass", "tree")) 
  labs(x="TIME",y="ToC")   
  theme_bw()

你會得到這樣的東西：

已編輯

在編輯問題和評論后：

NDVI %>% 
  mutate(target = as.factor(target)) %>% 
ggplot() 
  geom_bar(aes(x=date, y=TempAmb_Avg, fill = target, color = target),
           stat="identity") 
  geom_point(aes(x=date, y=30*cNDVI , color=target), 
             #position = position_dodge(),
             stat="identity") 
  geom_line(aes(x=date, y=30*cNDVI, color=target), 
            #position = position_dodge(),
            stat="identity") 
  geom_errorbar(aes(x = date, 
                    y = cNDVI ,
                    ymin=30*cNDVI -30*sd_NDVI, 
                    ymax=30*cNDVI  30*sd_NDVI, 
                    width=.2,,
                    color = target
                    ), 
                #position = position_dodge(), 
                stat="identity")  
  scale_color_manual(values = c(adjustcolor("blue", 0.5), adjustcolor("red", 0.5)), breaks = c("grass", "tree"))  
  scale_fill_manual(values = c(adjustcolor("blue", 0.3), adjustcolor("red", 0.3)), breaks = c("grass", "tree")) 
  labs(x="TIME",y="ToC")   
  theme_bw()

標籤：r图表2标准差风水

上一篇：R：使用手動顏色在geom_sf中設定某些值

下一篇：提取ggplot平滑函式并保存在資料框中