因此,我试图在R中进行卡方检验,比较鸟类种群趋势(稳定或下降)与饮食。两者都是分类变量,我使用长度函数来计算。
con.table <- tapply(britishbirds$trend, list(britishbirds$trend,
britishbirds$diet), length)
con.table
这是输出。
fish generalist invertebrates mammals/birds plants/grains
decline 5 NA 32 4 14
increase/stable 18 10 46 10 24
我知道我不能在卡方测试中使用NA,我希望这个函数输出0但是放了NA。有没有办法将其更改为0,或者我是否需要从数据集中省略此类别?