作者tcssh611503 (york)
看板R_Language
标题[问题] 用R跑卡方检定,并计算phi值
时间Mon May 29 21:49:50 2017
[软体熟悉度]:
入门(写过其他程式,只是对语法不熟悉)
[问题叙述]:
想利用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()
--
※ 发信站: 批踢踢实业坊(ptt.cc), 来自: 140.112.25.100
※ 文章网址: https://webptt.com/cn.aspx?n=bbs/R_Language/M.1496065793.A.F6F.html
1F:→ celestialgod: 当你的table不是2 x 2就会是NA 05/29 22:02
2F:→ celestialgod: The absolute value of the phi coefficient (only 05/29 22:03
3F:→ celestialgod: defined for 2 x 2 tables). 05/29 22:03
4F:→ celestialgod: 请确定你的table是 2 x 2 05/29 22:03
5F:→ celestialgod: 看你的自由度应该是 3 x 3 05/29 22:03
6F:→ celestialgod: 第二个推文是出自assocstats的manual 05/29 22:04
7F:→ celestialgod: 如果SPSS跑得出来要看一下它的phi是怎麽计算的 05/29 22:05
8F:→ celestialgod: 我看wiki Phi只有在2 x 2有定义 05/29 22:05
9F:→ tcssh611503: 阿,原来如此,感谢!!! 真得只能在2X2用 05/29 22:15