[问题] 用R跑卡方检定,并计算phi值

楼主: tcssh611503 (york)   2017-05-29 21:49:50
[软件熟悉度]:
入门(写过其他程式,只是对语法不熟悉)
[问题叙述]:
想利用R跑卡方检定,并算出phi值,使用vcd library中的assocstats()计算
但跑任何卡方,跑出来Phi-Coefficient值都是NA,
但在SPSS可以跑出来是,例如下面这个可以算出是0.206
[程式范例]:
#all$newage与all$part2D_6为两个类别变项
assocstats(table(all$newage,all$part2D_6))
执行结果:
X^2 df P(> X^2)
Likelihood Ratio 20.709 4 0.00036158
Pearson 20.667 4 0.00036861
Phi-Coefficient : NA
Contingency Coeff.: 0.202
Cramer's V : 0.206
[环境叙述]:
R version 3.3.2 (2016-10-31)
Platform: x86_64-w64-mingw32/x64 (64-bit)
Running under: Windows >= 8 x64 (build 9200)
[1] vcd_1.4-3
[关键字]:
卡方检定
vcd
assocstats()
作者: celestialgod (天)   2017-05-29 22:02:00
当你的table不是2 x 2就会是NAThe absolute value of the phi coefficient (onlydefined for 2 x 2 tables).请确定你的table是 2 x 2看你的自由度应该是 3 x 3第二个推文是出自assocstats的manual如果SPSS跑得出来要看一下它的phi是怎么计算的我看wiki Phi只有在2 x 2有定义
楼主: tcssh611503 (york)   2017-05-29 22:15:00
阿,原来如此,感谢!!! 真得只能在2X2用

Links booklink

Contact Us: admin [ a t ] ucptt.com