次のグラフを生成する次のコードがあります。
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))