R 用蒙特卡罗方法模拟卡方检验的p值

R 用蒙特卡罗方法模拟卡方检验的p值,r,testing,statistics,simulation,R,Testing,Statistics,Simulation,我正试图用R编写一个脚本,通过模拟Pearson卡方检验的临界值(p值),使用不同的alpha值来近似 我知道“chisq.test”中存在一个选项,但我想知道如何手工进行此模拟 例如: 请在检查代码(我不知道如何正确放置代码) 如果检查最后一部分(“模拟的p值”),您将看到脚本中获得p值的方式。我想这样做,但采用不同的alpha值 多谢各位 任何统计检验的p值计算(无论采用何种方法:经典,自举)与α值无关,如果你是指显著性水平。在决定接受或拒绝空假设时,您需要alpha值(如果p值小于所选的,

我正试图用R编写一个脚本,通过模拟Pearson卡方检验的临界值(p值),使用不同的alpha值来近似

我知道“chisq.test”中存在一个选项,但我想知道如何手工进行此模拟

例如:

请在检查代码(我不知道如何正确放置代码)

如果检查最后一部分(“模拟的p值”),您将看到脚本中获得p值的方式。我想这样做,但采用不同的alpha值


多谢各位

任何统计检验的p值计算(无论采用何种方法:经典,自举)与α值无关,如果你是指显著性水平。在决定接受或拒绝空假设时,您需要alpha值(如果p值小于所选的alpha,则拒绝空假设)。

如果您已按照脚本中所示进行模拟,并已导出模拟值向量
xsqsim
,然后,α水平
alpha
的临界值约为


如果你有一个小样本,你必须稍微小心一点,因为临界值应该是检验统计量q的值,这样观察值大于或等于q的概率等于α…

你的问题太宽泛了,你可以通过阅读一些书得到最佳答案。你有具体的问题吗?那么你的问题是什么?你试过什么。。。您是否深入研究了
chisq.test
的代码?是的,我会说得更具体一些。我试过这个:但是我不知道如何为p值的模拟取不同的alpha值。非常感谢。您可以修改您的问题以反映这一点吗?但该脚本似乎没有指定任何alpha值:相反,它按如下方式计算精确的p值:
mean(xsqsim>=xsq)
。基于此p值,您可以选择任何想要拒绝/失败的阿尔法级别来拒绝无效假设…谢谢!这就是我一直在寻找的答案。
quantile(xsqsim,1-alpha)