关于如何在R中执行chisq测试的空白

特德·莫斯比

我知道我是在某堂课上学到的，但现在我不记得了。

我有这样的数据：

 dput(tbl)
structure(c(160L, 7094L, 0L, 0L, 3287L, 373L, 164L, 2406L, 0L, 
0L, 33L, 0L, 0L, 0L, 0L, 122L, 20775L, 0L, 0L, 0L, 0L, 0L, 0L, 
417L, 0L, 1709L, 0L, 0L, 471L, 0L, 499L, 0L, 0L, 0L, 1029L, 4399L, 
3413L, 0L, 890L, 57L, 3185L, 0L, 0L, 1137L, 103L, 105L, 899L, 
0L, 0L, 7L, 0L, 0L, 0L, 0L, 69L, 8852L, 0L, 0L, 0L, 0L, 0L, 0L, 
53L, 0L, 776L, 0L, 0L, 222L, 0L, 193L, 0L, 0L, 0L, 312L, 1889L, 
1417L, 0L, 352L), .Dim = c(39L, 2L), .Dimnames = structure(list(
    c("ARSON", "ASSAULT", "BAD CHECKS", "BRIBERY", "BURGLARY", 
    "DISORDERLY CONDUCT", "DRIVING UNDER THE INFLUENCE", "DRUG/NARCOTIC", 
    "DRUNKENNESS", "EMBEZZLEMENT", "EXTORTION", "FAMILY OFFENSES", 
    "FORGERY/COUNTERFEITING", "FRAUD", "GAMBLING", "KIDNAPPING", 
    "LARCENY/THEFT", "LIQUOR LAWS", "LOITERING", "MISSING PERSON", 
    "NON-CRIMINAL", "OTHER OFFENSES", "PORNOGRAPHY/OBSCENE MAT", 
    "PROSTITUTION", "RECOVERED VEHICLE", "ROBBERY", "RUNAWAY", 
    "SECONDARY CODES", "SEX OFFENSES, FORCIBLE", "SEX OFFENSES, NON FORCIBLE", 
    "STOLEN PROPERTY", "SUICIDE", "SUSPICIOUS OCC", "TREA", "TRESPASS", 
    "VANDALISM", "VEHICLE THEFT", "WARRANTS", "WEAPON LAWS"), 
    c("Weekday", "Weekend")), .Names = c("", "")), class = "table")

我尝试做a，chisq.test(tbl)但结果以NA的形式返回，很可能是因为零。有人有见识吗？我正在计算工作日和周末之间的差异，犯罪类型可以组合为全部犯罪。

太空人

好吧，如果您认为其零值，请尝试不使用零值：

> chisq.test(tbl[tbl[,1]!=0,])

    Pearson's Chi-squared test

data:  tbl[tbl[, 1] != 0, ]
X-squared = 194.13, df = 16, p-value < 2.2e-16

这似乎产生了一些数字。

一对零足以扔掉它吗？

> chisq.test(rbind(tbl[tbl[,1]!=0,],c(0,0)))

    Pearson's Chi-squared test

data:  rbind(tbl[tbl[, 1] != 0, ], c(0, 0))
X-squared = NaN, df = 17, p-value = NA

Warning message:
In chisq.test(rbind(tbl[tbl[, 1] != 0, ], c(0, 0))) :
  Chi-squared approximation may be incorrect

是的。显然，犯罪总数为0，无法使工作日或周末变得更糟。我想您可以提交一个请求，使其删除零并给出警告，以确保这样做，但是我看不到该请求已实现。

本文收集自互联网，转载请注明来源。

如有侵权，请联系 [email protected] 删除。