0

次のグラフを生成する次のコードがあります。

colvec <-c("white", "white","gray85", "gray85", "gray58", "gray58", "gray33", "gray33","black", "black") 

ggplot(nut, aes(Date, Nitrate, group=Wetland, shape=Hydrology)) +
  geom_point(aes(fill=Wetland), colour="black", size=4)+
  scale_fill_manual(values=colvec) +
  scale_shape_manual(values=c(21,22))+
  facet_grid(. ~ Hydrology) +
  ylab ("Nitrate (mg/L) ") +
  theme(legend.position="none",
        panel.background = element_rect(fill='white', colour='white'), 
        panel.grid = element_line(color = NA),
        panel.grid.minor = element_line(color = NA),
        panel.border = element_rect(fill = NA, color = "black"),
        axis.text.x  = element_text(size=10, colour="black"),  
        axis.title.x = element_text(vjust=0.1),
        axis.text.y = element_text(size=12, colour="black"),
        axis.title.y = element_text(vjust=0.3))

ここに画像の説明を入力してください

stat_sum各日付の5ポイントを平均するような使い方がわかりません。プログラムは、それぞれ異なる陰影のあるポイント(湿地に対応)を独自の平均として扱います。湿地によるグループ化を維持し、これをさまざまな色合いで表示するだけでなく、各日付のすべてのy値の平均も表示したいと思います。

データ

     日付湿地水文学硝酸塩
117-6月1日パルス0.2647287              
218-7月1パルス0.1807388             
31-8月1日パルス0.9995910      
415-8月1日パルス0.6566667       
57-11月1日パルス0.2150000      
617-6月2日静的0.2134027      
718-7月2日静的0.1971669      
81-8月2日静的0.4774424       
915-8月2日静的0.3110000      
107-11月2日静的0.3333333       
1117-6月3パルス0.3369253       
1218-7月3パルス0.2056284       
131-8月3パルス0.6731924       
1415-8月3パルス0.5516667       
157-11月3パルス0.1853333      
1617-6月4日静的0.3293668      
1718-7月4日静的0.4664748       
181-8月4日静的0.4555003       
1915-8月4日静的0.3993333       
207-11月4日静的0.1133333       
2117-6月5日静的0.3497963     
2218-7月5日静的0.3618659      
231-8月5日静的0.3721719     
2415-8月5日静的0.2916667      
257-11月5日静的0.2526667      
2617-6月6日パルス0.2779667       
2718-7月6パルス0.7609531      
281-8月6パルス0.7177083       
2915-8月6パルス0.6610000       
307-11月6パルス0.2083333       
3117-6月7日パルス0.2232040      
32 18-Jul Seven Pulsed 0.3655621       
331-8月7日パルス0.7006131       
3415-8月7日パルス0.4753333      
357-11月7日パルス0.3206667    
3617-6月8日静的0.3339319       
3718-7月8日静的0.3286641       
381-8月8日静的0.4390918      
3915-8月8日静的0.3276667       
407-11月8日静的0.2446667       
4117-6月9日静的0.3456627      
42 18-Jul Nine Static 0.2519814      
431-8月9日静的0.3807550      
4415-8月9日静的0.3873333      
457-11月9日静的0.1663333       
46 17-Jun Ten Pulsed 0.4135023      
47 18-Jul Ten Pulsed 0.1921382       
481-8月10パルス0.3898374       
4915-8月10パルス0.2700000       
507-11月10日パルス0.1216667       

dput(nut)

structure(list(Date = structure(c(3L、4L、1L、2L、5L、3L、4L、
1L、2L、5L、3L、4L、1L、2L、5L、3L、4L、1L、2L、5L、3L、4L、1L、
2L、5L、3L、4L、1L、2L、5L、3L、4L、1L、2L、5L、3L、4L、1L、2L、
5L、3L、4L、1L、2L、5L、3L、4L、1L、2L、5L)、. Label = c( "1-Aug"、
「8月15日」、「6月17日」、「7月18日」、「11月7日」)、クラス=「ファクター」)、湿地=構造(c(5L、
5L、5L、5L、5L、10L、10L、10L、10L、10L、9L、9L、9L、9L、9L、
3L、3L、3L、3L、3L、2L、2L、2L、2L、2L、7L、7L、7L、7L、7L、6L、
6L、6L、6L、6L、1L、1L、1L、1L、1L、4L、4L、4L、4L、4L、8L、8L、
8L、8L、8L)、. Label = c( "Eight"、 "Five"、 "Four"、 "Nine"、 "One"、
"Seven"、 "Six"、 "Ten"、 "Three"、 "Two")、class = "factor")、Hydrology = structure(c(1L、
1L、1L、1L、1L、2L、2L、2L、2L、2L、1L、1L、1L、1L、1L、2L、2L、
2L、2L、2L、2L、2L、2L、2L、2L、1L、1L、1L、1L、1L、1L、1L、1L、
1L、1L、2L、2L、2L、2L、2L、2L、2L、2L、2L、2L、1L、1L、1L、1L、
1L)、. Label = c( "Pulsed"、 "Static")、class = "factor")、Nitrate = c(0.264728748、
0.180738787、0.995951021、0.656666667、0.215、0.213402705、0.197166881
0.477442378、0.311、0.333333333、0.33692531、0.205628403、0.67319236、
0.551666667、0.185333333、0.329366831、0.466474791、0.455500298、
0.399333333、0.113333333、0.349796312、0.361865927、0.372171941、
0.291666667、0.252666667、0.277966745、0.760953065、0.717708344、
0.661、0.208333333、0.223203974、0.365562124、0.700613059、0.475333333、
0.320666667、0.333931889、0.328664129、0.439091764、0.327666667、
0.244666667、0.345662714、0.251981433、0.380755049、0.387333333、
0.166333333、0.413502261、0.192138209、0.389837374、0.27、0.121666667
))、. Names = c( "Date"、 "Wetland"、 "Hydrology"、 "Nitrate")、class = "data.frame"、row.names = c(NA、
-50L))
4

1 に答える 1

4

あなたのデータなしで答えるのはちょっと難しいです(私は通常最初にテストするからです!)が、aggregate日付と水文学によって硝酸塩データを計算し、平均を計算してからgeom_point、別のdata.frame.

これは機能しますか...

nutm <- aggregate( Nitrate ~ Date + Hydrology , data = nut , FUN = mean )

p <-ggplot(nut, aes(Date, Nitrate, shape = Hydrology)) +
  geom_point( data = nut , aes( fill = Wetland , group=Wetland ), colour="black", size=4)+
  scale_fill_manual(values=colvec) +
  scale_shape_manual(values=c(21,22))+
  facet_grid(. ~ Hydrology) +
  ylab ("Nitrate (mg/L) ") +
  geom_point( data = nutm , aes( x = Date , y = Nitrate) , color = "red" , fill = "red" , size = 4 ) +
  theme(legend.position="none",
        panel.background = element_rect(fill='white', colour='white'), 
        panel.grid = element_line(color = NA),
        panel.grid.minor = element_line(color = NA),
        panel.border = element_rect(fill = NA, color = "black"),
        axis.text.x  = element_text(size=10, colour="black"),  
        axis.title.x = element_text(vjust=0.1),
        axis.text.y = element_text(size=12, colour="black"),
        axis.title.y = element_text(vjust=0.3))

p

これにより、各日付の赤い平均マークが付いたグラフが表示されます... ここに画像の説明を入力

于 2013-03-19T15:08:51.787 に答える