中央値データを使って複数の時系列プロットをプロットするR

Question

私はちょうどRを学ぶようになったので、これは冗長かもしれないので謝ります。私はこの添付をRでエクセルグラフにしたいと思っています。私のデータには、1951年から2016年までの貯水池の毎日の水位が含まれています。私は3つの系列（1951年から2013年、2014年から2015年、2016年）私はこれらのシリーズの中央値を1日にプロットしたいと思います。そして、中央値からfebuary 29を削除したいと思います。ここでhttps://docs.google.com/spreadsheets/d/1u1Whfp6VHXkZgrC0sVn_mT9XiVxszMhqlAszjZXzM1E/edit?usp=sharing

私がこれまで持っているものされています：ここに私のデータへのリンクです

LL <- read.csv("BSLL.csv") 
str(LL) 
LLpre <- filter(LL, year > "1952" & year <"2014") 
headtail(LLpre, n=3) 
medianLLpre = ddply(LLpre, .(month, day), summarise, level = median(level), na.rm = FALSE) 
LLpost <- filter(LL, year > "2013" & year < "2016") 
headtail(LLpost, n=3) 
medianLLpost = ddply(LLpost, .(month, day), summarise, level = median(level), na.rm = FALSE) 
LL2016 <- filter(LL, year == "2016") 
headtail(LL2016, n=3) 
medianLL2016 = ddply(LL2016, .(month, day), summarise, level = median(level), na.rm = FALSE)  

Excel Plot

Answer 1

library(data.table) 
Data <- fread("BSLL - BSLL.csv") 

は一時的な時間間隔のカテゴリを作成します。1 = 1951 2013,2 = 2014 2015、3 = 2016

Data[,Categories:=ifelse(as.numeric(year)<2014,1, 
         ifelse(as.numeric(year)<2016,2,3))] 

febなしの一時データを作成する29 th

Temp.Data <- Data[!month==2][!day==29] 

カテゴリとプロットで中央値を計算します。

Medians <- Temp.Data[,median(as.numeric(level)),by=.(month,Categories)] 
plot(x=Medians[,month], y=Medians[,V1],type="n", xlab="", ylab="", 
    xaxt="n", yaxt="n", bty="l") 
axis(2, tick=TRUE, labels=FALSE) 
axis(2, tick=FALSE, labels=TRUE, line=-0.5) 
axis(1, tick=TRUE, labels=FALSE) 
axis(1, tick=FALSE, labels=TRUE, line=-0.5) 
lines(x=Medians[Categories==1][,month], 
     y=Medians[Categories==1][,V1],type="l",lwd=1,col="red") 
lines(x=Medians[Categories==2][,month], 
     y=Medians[Categories==2][,V1],type="l",lwd=1,col="green") 
lines(x=Medians[Categories==3][,month], 
     y=Medians[Categories==3][,V1],type="l",lwd=1,col="blue") 
legend('topright','groups', 
     legend=c("1951 to 2013","2014 to 2015","2016"), bty="n",lty=c(1,1), 
     col=c("red","green","blue")) 

どういうことですか？プロットの詳細（サイズ、位置、色、ラベルなど）を細かく調整することができます

P.S今後、フルデータセットを一般公開することはお勧めしません。貢献者は、ここで作業してテストするデータの例が必要です