4

facet_grid()たとえば、ggplot2 を使用して 2 つのヒストグラムを取得しています。ここで、これらのヒストグラムにポアソン曲線 (2 つのヒストグラム プロット/グリッドの平均が異なる) と別の分布の 2 番目の曲線 (値の確率関数を手動で提供したい) を重ね合わせたいと考えています。これはどのように行うことができますか？

例の作成:

library(ggplot2)

value<-c(rpois(500,1.5))

group<-rep(c("A","B"),250)

data<-data.frame(value,group)

g1<-ggplot(data,aes(value))

g1+geom_histogram(aes(y=..count..),binwidth=1,position="identity")+facet_grid(.~group)

次は何？

または、格子パッケージを使用して実行できますか?

4

2 に答える 2

6

簡単な方法は、カウントの代わりに密度をプロットし、stat_function() を使用することです。

library(ggplot2)
value<-c(rpois(500,1.5))
group<-rep(c("A","B"),250)
data<-data.frame(value,group)
ggplot(data,aes(value)) + 
        geom_histogram(aes(y=..density..), binwidth=1,position="identity") + 
        facet_grid(.~group) + 
        stat_function(geom = "line", fun = dpois, arg = list(lambda = 1.5), colour = "red", fill = NA, n = 9)

カウントが必要な場合は、dpois の密度を「カウント」に変換する必要があります

ggplot(data,aes(value)) + 
        geom_histogram(aes(y=..count..), binwidth=1,position="identity") + 
        facet_grid(.~group) + 
        stat_function(geom = "line", fun = function(..., total){dpois(...) * total}, arg = list(lambda = 1.5, total = 250), colour = "red", fill = NA, n = 9)
于 2010-08-16T08:40:27.860 に答える
2

最近、同様の問題 (分布の比較) に直面したとき、どこから始めるべきかについてのアイデアを提供する透明なオーバーラップ ヒストグラムのコードを書きました。

于 2010-08-18T11:27:15.743 に答える