Rでの有意性テスト、1つの列の割合が単一変数内の他の列と大幅に異なるかどうかを判断します
-
25-10-2019 - |
質問
これはRの簡単なコマンドだと確信していますが、何らかの理由で、解決策を見つけるのに苦労しています。
RのCrosStabs(Table()コマンドを使用)の束を実行しようとしています。各タブには2つの列(治療と治療なし)があります。列間の違いがすべての行で互いに大幅に異なるかどうかを知りたいと思います(行は、調査からの回答の選択肢のほんの一握りです)。私は全体的な重要性に興味がありません。治療と治療なしを比較するCrosstab内でのみです。
このタイプの分析はSPSSで非常に簡単です(以下のリンクが私が話していることを説明します)が、Rで機能することはできません。
編集:私が意味することについての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
私はこのような表を持っています ^(TreatentVarの列でパーセント)。そして、各質問の選択肢(行)が治療0から治療までの大きな違いがあるかどうかを見たいです。したがって、上記の例では、私は4と2(行1)、3と3(行2)、1と3(行3)の間に大きな違いがあるかどうかを知りたい。したがって、この例では、Question1の選択肢は選択1と3の有意差である可能性があります(差は2であるため)が、選択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
警告に注意してください。これらのテストは、そのような少数に必ずしも適切ではありません。
他のヒント
あなたが探している機能はそうだと思います pairwise.prop.test()
. 。見る ?pairwise.prop.test
例として。