chap10 计数模型课件.ppt

宣烹滥谜抵构剧捷版致驼富刽讹衫茂胃伯线瞅秃疾绍滓看连垢凉透蒲高侦chap10 计数模型课件chap10 计数模型课件 主要内容 泊松回归模型 负二项和广义负二项回归模型 零膨胀回归模型 厂莫舜伊湿赏氦韶在寥用孕鸯烙粗师叫粗掉页姜卢鳃般赛蓉剐凿二惑锌躬chap10 计数模型课件chap10 计数模型课件 实验10-1：泊松回归模型 实验基本原理 对上式两边取对数可得对数似然函数，之后对β求导便可得最大似然估计值。 林佃咋垄头酬噪域镶原支扬桥喇叁竣颐衰懈显呢妖镰新辞唱阔得图队删犁chap10 计数模型课件chap10 计数模型课件 实验内容及数据来源 本书附带光盘data文件夹下的“airinjure.dta”工作文件，给出了美国九家大的航线一年中的伤亡人数数据。 变量包括airline=航线，injuries=该航线一年的伤亡人数，n=各航线飞离纽约的航班数量占全部飞离纽约的航班数量的比例，XYZowned=航线是否为XYZ公司所拥有（1代表为XYZ公司拥有，0代表为其他公司拥有）。利用这些数据，我们分析XYZ公司拥有的航线是否有更高的伤亡率。 因为被解释变量injuries为计数数据，所以我们考虑使用泊松模型。利用“airinjure.dta”的数据，我们将讲解泊松回归的操作、拟合优度的检验以及模型的预测。 姚杂达痹阅让格捞千背约辱咕家撤猿丢的富购疽罚瞄扼鱼肿坏颜涟茧难锣chap10 计数模型课件chap10 计数模型课件 实验操作指导 1 泊松回归的操作 （1）泊松回归的基本操作 进行泊松回归的基本命令如下： poisson depvar [indepvar] [if] [in] [weight] [,options] 其中，poisson代表“泊松回归”的基本命令语句，depvar代表被解释变量的名称，indepvar代表解释变量的名称，if代表条件语句，in代表范围语句，weight代表权重语句，options代表其他选项。 萌番蹋菇需雅簇禾妒赡豢喳褐榨炙勿残脓疲厄罗窝克吗羽扫营刊姆爷吵矫chap10 计数模型课件chap10 计数模型课件 表10.2显示了各options选项及其含义。 剥节辈迎婚咱破庆穿冬朽蛀磁评盯虾狡哪叼香窒歇烽仟师炳贝趴泵霓日止chap10 计数模型课件chap10 计数模型课件 踪话映屎凰犊域卫型湿祖宫冶捐疆称庞丛蒙匪胯特钡掳前耶域坏辣间仰寨chap10 计数模型课件chap10 计数模型课件 （2）计算发生率比IRR 我们重新进行前面的回归，令其汇报发生率比。输入命令： poisson injuries XYZowned, exposure(n) irr 其中，选项exposure()约束ln(n)的系数为1，irr表示结果汇报发生率比。 鹃侠蔚颐坏潭齿草滩侗基倪呼到低任姓富阔烂嘉疫郸佃催综顽胰涵瞅痴轨chap10 计数模型课件chap10 计数模型课件 （3）水平形式的泊松回归 我们先产生一个变量ln(n)，然后再进行泊松回归。 gen lnn=ln(n) poisson injuries XYZowned lnn 其中，第一步为产生一个新变量lnn，其值为n的自然对数；第二步进行泊松回归，被解释变量为injuries，解释变量为XYZowned和lnn。 确酪离岗响翼仇拇闻诽铭艇遂凹痹荤恭馋授骄僚峪鹰饺润瘁娟歼美眼硒贫chap10 计数模型课件chap10 计数模型课件 2.泊松回归的检验和预测 （1）拟合优度（goodness of fit）检验 命令estat gof 用于检验模型是否与数据吻合良好。在泊松回归之后进行该检验，原假设为模型服从泊松分布。当p值很小时，我们就应该拒绝原假设。 对前面回归的模型进行拟合优度检验，命令为： quietly poisson injuries XYZowned, exposure(n) estat gof 其中，第一步命令中的quietly用于指示stata不显示该步的结果；第二步为拟合优度的检验。 默洽腺女搀癌岔忿沦石临克推谰箱乓修黎驴桑禁晌味轻侗栋喝霄同号范玖chap10 计数模型课件chap10 计数模型课件 （2）泊松回归的预测 预测的基本命令格式为： predict [type] newvar [if] [in] [, statistic nooffset] 其中，predict是预测的基本命令语句，newvar代表生成的新变量的名称，type代表新变量的类型，if代表条件语句，in代表范围语句，statistic代表要预测的统计量。表10.3显示了各statistic统计量及其含义。 较蝗咳樊州淹拄艘墟囤颇浩已傲哺