r - ggplot2 を使用して積み上げ領域プロットの欠損値を修正する

Question

積み上げ棒グラフと面グラフの組み合わせでこの投稿を再現しようとしています。ただし、欠損値にはいくつか問題があります。

ここに私のデータがあります: https://www.dropbox.com/sh/pnkspwnn1qslm6u/JapTKCwqMS

私が実行しているのは;

    wa=read.table('wa_class.txt', sep="", header=F, na.string="0")
    names(wa)=c("Class","Jan","Feb","Mar","Apr","May","Jun","Jul","Aug","Sep","Oct","Nov","Dec")
    wam=melt(wa)
wam$variablen=as.numeric(wam$variable)

どのように見えるか

    > head(wam)
                   Class variable     value variablen
1         Actinobacteria      Jan 38.115163         1
2          Flavobacteria      Jan        NA         1
3        Sphingobacteria      Jan  3.640469         1
4    Alphaproteobacteria      Jan 13.631663         1
5 Betaproteobacteria_b28      Jan  3.718671         1
6     Betaproteobacteria      Jan 14.732354         1

ggplot(na.omit(wam[,c("Class","value","variablen")]), aes(wam,x=variablen, y=value, fill=Class)) + geom_area(color="black") + geom_linerange(aes(ymax=value), position="stack") + scale_x_continuous(breaks=1:max(wam$variablen)) + labs(title="Water", x="Month", y="Relative abundance (%)")

...だから、プロットする変数で na.omit を使用して欠損値を修正しようとしました。ただし、たとえば互いにオーバーライドしているプロット内のレイヤーを取得します（ドロップボックスフォルダーを参照）。

それを修正したこの投稿（ドロップボックスフォルダーを参照）を見つけましたが、1つしかないようでした。また、データへのリンクが切れているため、再現できませんでした。

どんな助けでも本当に感謝します!

ありがとう、

ジョー