在R中的显著性检验，确定单变量内一个列的比例是否与另一个列显著不同。

Question

我确定这是R中的一个简单命令，但出于某种原因，我无法找到解决方案。

我试图在R中运行一堆交叉表（使用table()命令），每个表都有两列（治疗和非治疗）。 我想知道所有行的列之间是否存在显着差异（行是来自调查的几个答案选择）。 我不关心总体显着性，只关心比较治疗与非治疗的交叉表内部的显着性。

这种分析在SPSS中非常容易（下面的链接说明了我所说的），但我似乎无法在R中使其工作。 你知道我怎么做吗？

http://help.vovici.net/robohelp/robohelp/server/general/projects_fhpro/survey_workbench_MX/Significance_testing.htm

编辑后: 这是一个关于我所说的R语言示例:

 treatmentVar <-c(0,1,0,1,0,1,0,1,0,1,0,1,0,1,0,1) # treatment is 1 or 0
 question1 <-c(1,2,2,3,1,1,2,2,3,1,1,2,2,3,1,3) #choices available are 1, 2, or 3
 Questiontab <- table(question1, treatmentVar)
 Questiontab

我有这样的表格^（按treatmentVar列百分比），我想知道从treatment 0到treatment 1每个问题选项（行）之间是否存在显著差异。因此，在上面的示例中，我想知道第1行的4和2（行1），第2行的3和3（行2），以及第3行的1和3之间是否存在显著差异。因此，在这个例子中，问题1的选择可能在选择1和3之间存在显著差异（因为差异为2），但是选择2的差异不是因为差异为零。最终，我试图确定这种类型的显著性。希望这可以帮助你。

谢谢！

- Captain Murphy

---

如果您提供一个可重现的示例，包括一个小的数据集，那么回答您的问题将会更容易。 - joran

抱歉，我刚刚用R语言的例子修改了我的问题。 - Captain Murphy

Answer 1

我认为你需要的函数是 pairwise.prop.test()。请参考 ?pairwise.prop.test 中的示例。

- John Colby

---

谢谢John，但那不是我想要的。这是我的错，我应该在问题中更清楚。我现在已经用一个R示例更新了我的问题。希望这能澄清问题。 - Captain Murphy

Answer 2

使用您的示例，要么使用chisq.test或prop.test（在这种情况下等效）:

> chisq.test(Questiontab)

        Pearson's Chi-squared test

data:  Questiontab 
X-squared = 1.6667, df = 2, p-value = 0.4346

Warning message:
In chisq.test(Questiontab) : Chi-squared approximation may be incorrect
> prop.test(Questiontab)

        3-sample test for equality of proportions without continuity
        correction

data:  Questiontab 
X-squared = 1.6667, df = 2, p-value = 0.4346
alternative hypothesis: two.sided 
sample estimates:
   prop 1    prop 2    prop 3 
0.6666667 0.5000000 0.2500000 

Warning message:
In prop.test(Questiontab) : Chi-squared approximation may be incorrect

请注意警告：这些测试不一定适合如此少量的数据。

- Brian Diggs