ggplotを特定の順序でプロットして、オブジェクトが重なったときに表示されるものを制御します。各データ行は、2つのジオメトリレイヤーのコンポジットにマップされます。プロットの詳細には、これが必要です。私はこれを行うためにループを使用してきましたが、それは非常に遅いです。もっと良い方法があるのだろうか?例えばggplot2で複数のジオメトリを持つビジュアルオブジェクトのプロット順序を制御する
d = data.frame(x=c(1,5,2,2,4,2), y=c(1,1,4,3,3,5), grp=c(1,1,1,2,2,2))
ggplot(d, aes(x, y, group=grp)) +
geom_polygon(aes(fill = factor(grp))) +
geom_line(size=1)
各ポリゴンの線は、そのポリゴンでプロットしなければならない - そう例えば赤いポリゴンの線は青のポリゴンによって隠されるべきです。 geom_polygon
とgeom_line
の両方が同じデータセットを使用しているときに、ループなしでこれを達成する方法はありますか?
編集:ルーピング方法..ここ
は、私が使用したループメソッドです。パフォーマンスを比較するためのより良いデータセットを追加しました。どちらも私のマシンで約5.6秒かかる。比較すると、典型的なアプローチ(ggplot(d, aes(x, y, fill=factor(grp))) + geom_polygon() + geom_line(size=1)
)は0.45秒かかる。
d = data.frame(x = sample(-30:30,99,rep=T) + rep(sample(1:100,33),each=3),
y = sample(-30:30,99,rep=T) + rep(sample(1:100,33),each=3),
grp = rep(1:33,each=3))
# Method 1 - for loop
p = ggplot()
for(g in unique(d$grp)){
dat = subset(d, grp == g)
p = p + geom_polygon(data=dat, aes(x, y, fill = factor(grp))) +
geom_line(data=dat, aes(x, y), size=1)
}
print(p)
# Method 2 - apply
ggplot() + lapply(unique(d$grp), FUN=function(g){
dat = subset(d, grp == g)
list(geom_polygon(data=dat, aes(x, y, fill = factor(grp))),
geom_line(data=dat, aes(x, y), size=1))
})
あなたが直面している問題は、ラインと関係していると思います。 geom_line()を削除すると、目的の結果が得られます。また、geom_polygon()の前にgeom_lineを置くと、あなたが望む出力が得られるとわかります –
あなたが使っていたfor-loopを表示できますか? – Heroka
@geotheory申し訳ありませんが解決策はありません。私が持っている唯一のアイデアは、あなたがラインを必要としていたときに、両側で少し大きめのグループごとに余分なポリゴンを作成し、それらを最初に黒でプロットし、次にカラーで上書きすることです。しかし、私はこれがあなたのイメージを生成するためにループを実行するのに要するよりもプログラム/テストに時間がかかると思う。 – Heroka