ggplotを使用した時系列の複数変数の平均のプロット

Question

aからkまでの複数の変数の時系列データを含むファイルがあります。

変数a〜kの平均値を時系列以上にプロットしたグラフを作成したい場合、その平均線は毎日の最大変動と最小変動を表す平滑化領域を追加します。

信頼区間のようなものですが、平滑化されたバージョンです。ここで

は、データセットです： https://dl.dropbox.com/u/22681355/co.csv

、ここでは、私が持っているコードは、これまでのところです：

library(ggplot2) 
library(reshape2) 
meltdf <- melt(df,id="Year") 
ggplot(meltdf,aes(x=Year,y=value,colour=variable,group=variable)) + geom_line() 

Answer 1

これは、ブートストラップの95％信頼区間を示している：

ggplot(meltdf,aes(x=Year,y=value,colour=variable,group=variable)) + 
    stat_summary(fun.data = "mean_cl_boot", geom = "smooth") 

これは、彼はすべての変数+ -1SDのすべての値の平均：

ggplot(meltdf,aes(x=Year,y=value)) + 
    stat_summary(fun.data ="mean_sdl", mult=1, geom = "smooth") 

あなたは年が変数より手段とSDを計算する前に意味が、私はあなたにそれを残して計算する場合があります。

しかし、私は、分布が明らかに対称ではないので、ブーストラップの信頼区間がより賢明であると信じています。それはまた狭くなります。 ;）

もちろん、値をログに変換することもできます。