作者sodboy (伯伯)
看板R_Language
标题[问题] multiple proportion test
时间Tue Oct 22 15:30:26 2019
大大好
如果我今天要检定多组间发生率是否有无差异,使用prop.test()是正确的嘛?
自变项是性别加国籍,
依变项是发生率(十万分之一为单位)
另,如何判断发生率是否为稀有事件呢,我觉得我发生率好像蛮低,prop.test()较正要
用什麽语法呢?
谢谢!
--
※ 发信站: 批踢踢实业坊(ptt.cc), 来自: 1.200.216.214 (台湾)
※ 文章网址: https://webptt.com/cn.aspx?n=bbs/R_Language/M.1571729428.A.289.html
1F:→ andrew43: 发生率应可回推发生次数与总次数,则考虑存在offset项的10/22 17:28
2F:→ andrew43: poisson regression。10/22 17:29
3F:→ andrew43: 其中自变数为国籍。之後国籍间的多重比较也不难写。10/22 17:29
4F:→ andrew43: 至於prob.test就只是个卡方同质检验,也是要使用次数而10/22 17:32
5F:→ andrew43: 非发生率,这细节可看 ?prob.test10/22 17:34
6F:→ andrew43: 若你只有发生率但没有各国内的总次数,那可能没检验可做10/22 17:35
谢谢
※ 编辑: sodboy (1.200.216.214 台湾), 10/22/2019 17:57:14
7F:→ andrew43: 差不多。可额外考虑纳入年份和性别,以及因子内多重比较10/22 18:16
※ 编辑: sodboy (1.200.216.214 台湾), 10/23/2019 00:40:05
大大好
https://i.imgur.com/r9w8SEN.jpg
https://i.imgur.com/oQVIaUW.jpg
https://i.imgur.com/tEV9zCh.jpg
https://i.imgur.com/oHPZzRg.jpg
变异数大过平均值很多
发生次数画出来也不像poisson分布
还能使用poisson regression来分析嘛
还是须要换别种检定呢
谢谢
※ 编辑: sodboy (110.50.156.241 台湾), 10/28/2019 12:13:21
※ 编辑: sodboy (110.50.156.241 台湾), 10/28/2019 12:14:05
8F:→ andrew43: sandwich修正标准误,或是负二项式回归。 10/30 01:08
9F:→ andrew43: poisson regression之後有goodness-of-fit可以看看 10/30 01:08
10F:→ andrew43: 例如Residual deviance对应其df的右尾卡方检验 10/30 01:10