作者wahaha2 (傀儡:影武者)
看板Statistics
标题[问题] 关於poisson regression medel的建立以及检定
时间Sat Mar 26 18:04:10 2011
想要研究某种疾病在不同地区以及不同年龄层的发生率是否有差异
假设 分成南区和北区 年龄分三组(<35, 35~74, >75)
每个地区每个年龄层的疾病发生率也收集好了 (总共有6笔资料：住在南区且年龄<35的发
生率, 住在南区且年龄35~74的发生率, ....依此类推)
现在想要建立poisson regression model
我的想法是先单独看某年龄层，然後建立该年龄层的model如下
λ(发生率) = exp (β0 + β1X) (北区X=1, 南区X=0)
而因为有三个年龄层所以会有三个model
现在问题来了
1. 要先检定模型适合度，跟据我现在只拥有6笔资料，要怎麽检定?
2. 假设检定结果是适合的，要如何检定年龄、地区对於发生率是否有影响，可不可以说
一下基本原理 谢谢
http://www2.cmu.edu.tw/~biostat/online/teaching_corner_050-2.pdf
上面是使我产生困惑的文章 里面只有说用SAS这种软体分析 我想知道详细方法
--
※ 发信站: 批踢踢实业坊(ptt.cc)
◆ From: 120.126.71.150
1F：→ andrew43:为什麽要分成三个model? 03/26 18:23
2F：→ wahaha2:所以只要加上B2X 然後X=0, 1, 2 分别代表年龄 就好了吗? 03/26 18:31
3F：→ andrew43:至於基本原理，你可以查google，看看英文的网页。非常多 03/26 18:31
4F：→ andrew43:那篇文章里不是交待得很清楚吗?其实我不懂你的问题是什麽 03/26 18:33
5F：→ wahaha2:我不懂位什麽这个时候要用poisson而不用logistic 03/26 18:35
6F：→ wahaha2:还有就是beta0 1 2求出来後如何帮助判断是否有差异 03/26 18:37
7F：→ wahaha2:是只要非0就有差异吗 还是有什麽检定方法 03/26 18:38
8F：→ andrew43:嗯。用logistic的话，还要知道多少人没病。这点差很多。 03/26 18:54
9F：→ yhliu:如果只有 "发生率", 基本上你能做的不是 Poisson model. 03/26 21:04
10F：→ yhliu:Poisson regression model 用的资料是 "发生数". 03/26 21:04
11F：→ yhliu:另外, 事实上可以说没办法检查 Poisson model 是否合适. 03/26 21:05
12F：→ yhliu:事实上用 rate 建模是对的, 但只有 rate 是不够的. 03/26 21:07
13F：→ yhliu:而 model 很简单, ln(rate) = GM + Region + Age 03/26 21:07
14F：→ yhliu:或 ln(rate) = GM + Region + Age + Region*Age 03/26 21:08
15F：→ yhliu:Region 已是 0-1 变数, Age 考虑用两个 dummy 或给予代表值. 03/26 21:09
16F：→ yhliu:若 Age 用两个 dummy variable 表现, 则 03/26 21:10
17F：→ yhliu:有 Region*Age 的 model 是饱和模型, 没有残差. 03/26 21:10
18F：→ wahaha2:谢谢楼上 请问模型建立後 要检定年龄或地区是否造成显着 03/26 21:59
19F：→ wahaha2:是不是要看该变项之exp(系数) 是否等於1 03/26 22:00
20F：→ wahaha2:如果是的话 那麽exp(系数)之信赖区间要怎麽算 03/26 22:01
21F：→ yhliu:如果你只有 rate, 基於各 rate 的误差不等, 我无可建议. 03/26 23:52
22F：→ yhliu:如果你有算出 rate 的人口基数, 用 Poisson regression 的 03/26 23:53
23F：→ yhliu:rate model, 请参阅 "类别资料分析" 或 "广义线模" 或 03/26 23:54
24F：→ yhliu:"Poisson regression model" 的专书. 03/26 23:55
25F：→ kingbee:先列出地区和年龄的列联表 03/28 09:29
26F：→ kingbee:把六笔资料用两个变项去coding 再去跑SAS的genmod 03/28 09:30
27F：→ kingbee:人口用的发生率 是给定人口有一个单位 在计算上作为offset 03/28 09:32