library(ggplot2)
#2つの普通のサイコロを振って出た目を足す
roll01 <- function() {
die <- 1:6
dice <- sample(die, size=2, replace=TRUE)
sum(dice)
}
#2つの、出る目の確率がちょっと違うサイコロを振って出た目を足す
roll04 <- function() {
die <- 1:6
dice <- sample(die, size=2, replace=TRUE, prob=c(1/8, 1/8, 1/8, 1/8, 1/8, 3/8))
sum(dice)
}
#各々を10,000回やってみる
rolls01 <- replicate(10000, roll01())
rolls04 <- replicate(10000, roll04())
#その結果をヒストグラムで描く
# ylim(0, 1700)を指定することで縦の単位が揃う
a <- qplot(rolls01, binwidth = 1) + ylim(0, 1700)
b <- qplot(rolls04, binwidth = 1) + ylim(0, 1700)
a
b
2つのヒストグラムを横に並べる。
#まあ、縦でいいかもしれないけど、横に並べてみたいときもある。
#先ほどの ylim(0, 1700)によって縦軸が揃っている。
#これがないと縦軸が揃わないのでちょっと見にくい。
library(grid)
pushViewport(viewport(layout=grid.layout(1,2)))
print(a, vp=viewport(layout.pos.row = 1, layout.pos.col = 1))
print(b, vp=viewport(layout.pos.row = 1, layout.pos.col = 2))