R - ggplot2 geom_area 1つのカテゴリのみが表示されますか？

Question

私は許してください...

私はggplot2を使ってRでtimeseriesプロットをしようとしています。私は、与えられた日付の値の下に記された2つの折れ線グラフを表示したいと思います。私はgeom_area（position = "identity"）関数でこれをしようとしています。

ただし、グラフには1色しか表示されません（両方とも凡例に表示されます）。私はmelt（）を使ってデータを溶かし始めました.3つのカラム（X =時間、変数=地下水、値=地下水の標高）で作業しています。以下は私のコードの簡略化されたバージョンと、私が得るもののスクリーンショットです。

Bank01MWtest<-data.frame(X=(c(1,2,2,1)),variable=(c("MW-01A","MW-01A","MW-01B","MW-01B")),value=(c(576,571,584,580))) 

ggplot(data=Bank01MWtest, aes(x=X, y=value,group=variable))+geom_area(position="identity", aes(fill=variable))+geom_line(aes(color=variable))+coord_cartesian(ylim=c(570,590)) 

2色を表示したい。 MW.01Aライン以下の1色とMW.01Bライン以下の1色。

任意のヘルプ？

Answer 1

私はがgeom_ribbonのggplot2に置き換えられていると考えていますので、後者は私の解決策で使用します。このソリューションでは、データを長めから広範囲に再構成する必要があり、それぞれの凡例カテゴリに独自の列を付ける必要があります。 reshape2パッケージ内のdcast関数でこれを行います。

異なるymaxのレイヤーを追加し、fillオプションで凡例ラベルを割り当て、次にscale_fill_manual機能を使用して凡例を追加します。

library(ggplot2) 
library(reshape2) 

Bank01MWtest<-data.frame(X=sample(c(1,1,2,2)), 
        variable=sample(c("MW01A","MW01A","MW01B","MW01B")), 
        value=sample(c(576,571,584,580))) 

### Note above I modified your category labels by getting rid of the "-" sign 
### so that they can be used as variable names below. 

dat = dcast(Bank01MWtest, X~variable) 

ggplot(data=dat, aes(x=X)) + 
    geom_ribbon(aes(ymin=0, ymax=MW01A, fill="MW01A")) + 
    geom_ribbon(aes(ymin=0, ymax=MW01B, fill="MW01B")) + 
    scale_fill_manual("", values=c("green", "blue")) + 
    coord_cartesian(ylim=c(570,590)) 

Answer 2

いくつかの合成的に生成Bank01MWtestデータセットで、geom_areaでこれを試してみてください：

head(Bank01MWtest) 
    Time variable value 
1 2016-07-01 MW-01A 582.5482 
2 2016-07-02 MW-01A 580.5652 
3 2016-07-03 MW-01A 582.3305 
4 2016-07-04 MW-01A 583.3122 
5 2016-07-05 MW-01A 576.3432 
6 2016-07-06 MW-01A 584.4086 

tail(Bank01MWtest) 
     Time variable value 
195 2016-10-03 MW-01B 573.8355 
196 2016-10-04 MW-01B 575.3218 
197 2016-10-05 MW-01B 570.8007 
198 2016-10-06 MW-01B 572.3415 
199 2016-10-07 MW-01B 575.3291 
200 2016-10-08 MW-01B 578.0055 

ggplot(data=Bank01MWtest, aes(x=Time, y=value,group=variable))+ 
    geom_area(position='identity', aes(fill=variable), alpha=0.2)+ 
    scale_x_date(date_breaks= "1 month", date_minor_breaks = "15 days", date_labels = "%b", 
       limits = c(min(Bank01MWtest$Time),max(Bank01MWtest$Time))) + 
    geom_line(aes(color=variable))+coord_cartesian(ylim=c(570,590))