7

これはRの簡単なコマンドだと思いますが、何らかの理由で解決策を見つけるのに苦労しています。

Rで（table（）コマンドを使用して）一連のクロスタブを実行しようとしていますが、各タブには2つの列（治療と治療なし）があります。列間の違いがすべての行で互いに有意に異なるかどうかを知りたいです (行は調査からの回答の選択肢の一部です)。全体的な有意性には興味がありません。治療と治療なしを比較するクロスタブ内のみです。

この種の分析は SPSS では非常に簡単ですが (以下のリンクで私が話していることを説明しています)、R ではうまくいかないようです。

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

EDITED：これは私が何を意味するかについての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

私はこのようなテーブルを持っています ^ (治療変数の列ごとにパーセンテージ)。治療 0 から治療 1 までの各質問の選択肢 (行) に有意差があるかどうかを確認したいと思います。したがって、上記の例では、 4 と 2 (行 1)、3 と 3 (行 2)、1 と 3 (行 3) の間に有意差があるかどうかを知りたいとします。したがって、この例では、質問 1 の選択肢は、選択肢 1 と選択肢 3 で有意差がある可能性があります (差が 2 であるため) が、差がゼロであるため、選択肢 2 の差はそうではありません。最終的に、私はこのタイプの重要性を判断しようとしています。それが役立つことを願っています。

ありがとう！

4

2 に答える 2

9

あなたが探している機能はpairwise.prop.test(). 例については、を参照?pairwise.prop.testしてください。

于 2011-11-28T20:58:07.547 に答える
5

あなたの例を使用すると、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

警告に注意してください。これらのテストは、このような少数の場合には必ずしも適切ではありません。

于 2011-11-28T21:53:48.153 に答える