次のようなデータがあります。
type1 <- c(rep(1,49), rep(2, 30), rep(3, 20), rep(4,1))
type2 <- c(rep(1,10), rep(2, 20), rep(3, 40), rep(4,20), rep(5,10))
type3 <- c(rep(5,49), rep(4, 30), rep(3, 20), rep(2,1))
dat2 <- data.frame(dens = c(type1, type2, type3), lines = rep(c("a", "b", "c"),
each = 100))
ヒストグラムは次のとおりです。
require(ggplot2)
ggplot(dat2, aes(x = dens, fill = lines)) + geom_bar(position="dodge")+ theme_bw()
これは離散的に測定されるため、密度プロットは正確には機能しません。
ggplot(dat2, aes(x = dens, fill = lines)) + geom_density(alpha = 0.5)
また、スムーズな理論的分布を生成したいと考えています。次のようなベータ版配布を試しましたが、上記のデータ (dat2) にどのように適合させることができるかわかりません。
dat3 <- data.frame(dens = c(rbeta(4000000, 1, 3, ncp = 0),
rbeta(4000000, 2, 2, ncp = 0), rbeta(4000000, 3, 1, ncp = 0))
, lines = rep(c("a", "b", "c"), each = 4000000))
#Plot.
ggplot(dat3, aes(x = dens, fill = lines)) + geom_density(alpha = 0.5) + theme_bw()
