R语言——第八课

"一些统计分布"
# 二项分布
pbinom(5,8,0.9)
pbinom(54,100,0.5)-pbinom(45,100,0.5) # P(45<X<55)
qbinom(0.05,16,0.5) # 5%置信区间
rbinom(100,20,0.2) # 生成100个服从b(20,0.2)的随机数

# 正态分布
pbnorm(1.5)
qnorm(0.05)
qnorm(0.05,1,2) 
rnorm(50,4,2) # 生成服从N(4,2)的随机数



"抽样"
set.seed(1) # 生成种子数确保每次生成的相同
sample(1:6,replace = TRUE) # 从1到6放回抽样，默认不放回抽样
cards = paste(rep(c('A',2:10,'J','Q','K'),4),c('H','D','S','C'))
sample(cards,5)

k = 100
numcases = 100000
min = 1
max = 6
sumd = rep(0,numcases)
for(i in 1:k){
  dice = sample (1:6, numcases, replace = TRUE)
  sumd = sumd + dice
}
x.bar = sumd/k
hist(x.bar, prob = T, xlim = range(min-0.5,max+0.5),main=paste(numcases,
"rolls of a pair of dice"))