プログラミングによる確率分布の理解

二項分布

表が出る確率がpのコインをn回投げて,表がk回出る確率の分布.

Rでは,dbinom(k,n,p)という関数を使うことで,\(B(n,p)\) において\(k\)が出る確率を求めることができる(参考ウェブサイト).

dbinom(k,n,p)

Rでは,rbinom(size,n,p)という関数を使うことで,\(B(n,p)\) からsize個の乱数をサンプリングすることが出来る.

rbinom(size,n,p)

発展問題1:パラメータがそれぞれ \(n=50, p=0.3\) の 二項分布\(B(50,0.3)\) において,表が出る回数を \(k\) と表す.このとき, \(k=20\)となる確率を求めよ.

発展問題2:パラメータがそれぞれ \(n=50, p=0.3\) の 二項分布\(B(50,0.3)\) から,20個の乱数をサンプリングせよ.

解答

# EIP 第10回課題
# URL:https://fnakai.web.nitech.ac.jp/?p=589
# 解答

#問1:以下の四則演算を実行し,計算結果が出るかどうか確かめてください.

7+17
17-4
17*28
28/4
0.9^4

##計算結果の代入
eip <-((10^2)+7)/(1+(2*3))
((10^2)+7)/(1+(2*3)) -> eip
eip =((10^2)+7)/(1+(2*3))
assign("eip",((10^2)+7)/(1+(2*3)))

#eipの中身を見る
eip

#問3:4と5の二つの要素からなるベクトルを作成し,オブジェクト「e1」に代入せよ.

e1 <- c(4,5)
e1

#発展問題1 
dbinom(20,50,0.3)

#発展問題2
rbinom(20,50,0.3)