R中的自定义Fisher精确检验

Question

前面是初学者的问题！
（花了很多时间之后，找不到直接的解决方案..）

在尝试了所有相关帖子后，我似乎找不到答案，也许是因为我的问题非常基本。

我想对我的数据运行 fisher.test（无论数据，对我来说并不重要 - 我的是来自 QR33 的 Rubin 的儿童电视工作坊 - http://www.stat.columbia.edu/~cook/qr33.pdf）它必须模拟完全随机的实验。

我的假设是，在这种情况下，RCT 意味着所有单位都有相同的概率被分配到治疗（1/N）。 （当然，如果我错了，请纠正我。谢谢）。

我被要求创建一个自定义函数，我的函数必须包含以下参数：

1. 治疗观察（矢量）
2. 对照观察（向量）
3. 一个标量，表示尖锐零假设的值，例如零；和
4. 函数应运行的模拟实验数。

在 R 的 fisher.test 中挖掘时，我发现我可以指定 X、Y 和许多其他参数，但我不确定以下内容：

1. Y 是什么意思？ （即，a factor object; ignored if x is a matrix. 根据统计意义没有提供信息）。
2. 如何指定我的零假设？ （即，如果我不想使用 0。）我看到 there is a class "htest" with null.value 但如何在函数中使用它？
3. Reg 模拟次数，我的计划是通过循环运行所有内容 - 听起来很昂贵 - 有什么想法可以更好地编写它吗？

感谢您的帮助 - 我相信这不是一件容易的事，希望对许多人有用。

干杯，

NB - 以下解释令人不满意：

https://www.r-bloggers.com/contingency-tables-%E2%80%93-fisher%E2%80%99s-exact-test/

https://stats.stackexchange.com/questions/252234/compute-a-fisher-exact-test-in-r

https://stats.stackexchange.com/questions/133441/computing-the-power-of-fishers-exact-test-in-r

https://stats.stackexchange.com/questions/147559/fisher-exact-test-on-paired-data

Answer 1

我对我来说并不是完全清楚的，因为你想做的事情一定是正确的东西（这对统计数据来说是一个很好的问题），但我会解决R问题。 p >

1. 如“详细信息”部分的开始，R提供两个 EM>指定您的数据。

   您可以1.提供参数x 987654322 @ A的次次计数表（省略任何用于y），或者您可以向个人提供观察，作为指示行和列类别的两个向量（它没有这是哪个）;每个传染媒介包含x和y。 [我不确定为什么它也不会让你在x作为计数矢量和y作为数据框架，但它很容易转换]

2. 具有fisher测试，下核假设（有条件地对边缘）观察类别成为可交换的是独立性，但您可以选择使其成为一个或两个尾被的（通过alternative Argument） p >

3. 我不确定我清楚地了解模拟方面，但我几乎从不使用循环进行模拟（不是效率，而是为了清晰度和简洁）。 replicate非常适合进行模拟。我大致每天使用它，有时多次。