binomial检验_二项分布检验的R和SAS实现

R代码

### Two-sided binomial test

pvalue_binom2test

pvalue=dbinom(x, n, p)

for (i in 0:n) {

if (dbinom(i, n, p)

else prob=0

pvalue

}

print(pvalue)

}

# left side

pbinom(52, 60, 0.75)

binom.test(52, 60, p=0.75, alternative="less")

prop.test(52, 60, p=0.75, alternative ="less", correct=F)

# right side

pbinom(51, 60, 0.75, lower.tail = F)

binom.test(52, 60, p=0.75, alternative="greater")

prop.test(52, 60, p=0.75, alternative ="greater", correct=F)

#two sides

prop.test(52, 60, p=0.75, alternative ="two.sided", correct=F)

binom.test(52, 60, p=0.75, alternative="two.sided")

pvalue_binom2test(52, 60, 0.75)

### Produce different p-