R语言 在10次抛硬币中有统计意义的正面数

kcwpcxri  于 2023-03-05  发布在  其他
关注(0)|答案(1)|浏览(141)

我有一个问题,我不知道如何回答在一个问题表我已经得到了:
“一个实验者抛硬币,收集最初的一批N=10,如果这些结果在零假设检验中不显著,那么她将继续下去,直到N=20。
显著性水平为alpha = 0.05,硬币是公平的(p=0.5)。对于公平硬币的10次翻转,正面的数量是多少具有统计学显著性?
在对N=20重复此操作后,N=10和N=20的哪些结果组合具有显著性?
我知道你可以使用dbinom()函数来得到概率,然后你可以像这样把概率加在一起:

dbinom(0, 10, 0.5) + dbinom(10, 10, 0.5)  #probability of getting 0 heads and 10 heads for N=10

但我不知道如何得到显著的正面数的统计显著性,我想一旦N=10的情况解释给我听,我就可以解出问题的最后两个部分,在R中如何解出?

unftdfkk

unftdfkk1#

二回答第一部分问题:
dbinom将给予概率密度函数(pdf),表示观察到随机实验特定结果的可能性。

# What's the probability of observing exactly 5 heads when you flip the coin 10 times?
dbinom(5, size=10, p=0.5)

然后,你可以为第一个参数插入所有可能的值,并观察哪个结果符合你的p阈值。
但更优雅的是,这就是累积分布函数(cdf)的作用,在R中,它是通过qbinom函数实现的。

# Assuming the null-hypothesis,
# what is the number of heads X such that when repeating this experiment

# you will with 95% probability not observe a higher number of heads
qbinom(0.95, size=10, prob=0.5) # X = 8

# or: you will with 95% probability not observe a lower number of heads
qbinom(0.05, size=10, prob=0.5) # X = 2

为了加深理解,你可以反过来测试:

# What's the probability of observing exactly 8 heads when you flip the coin 10 times?
binom(x=8, size=10, prob=.5) # 0.04394531
binom(x=2, size=10, prob=.5) # 0.04394531

正如预期的那样,假设显著性水平为5%,得出的概率是显著的。

相关问题