1

我正在尝试在 R 中编写一个脚本，该脚本允许通过模拟 Pearson Chi Squared 测试的临界值（p 值）来近似，并采用不同的 alpha 值。

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

例如：

请在http://www.biostat.wisc.edu/~kbroman/teaching/stat371/comp21.R查看代码（我不知道如何正确放置代码）

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

非常感谢你！

4

2 回答 2

2

如果您的意思是显着性水平，那么任何统计测试的 p 值的计算（无论方法：经典、自举）都与 alpha 值无关。在决定接受或拒绝原假设时，您需要 alpha 值（如果 p 值小于选择的 alpha，则拒绝原假设）。

于 2012-06-05T13:05:57.407 回答
1

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

quantile(xsqsim,1-alpha)

如果你有一个小样本，你必须要小心一点，因为临界值应该是检验统计量q的值，使得观察值大于或等于 q的概率等于 alpha ...

于 2012-06-05T13:28:52.670 回答