统计学系列条目
迴歸分析
模型
線性回歸 简单线性回归 普通最小二乘法（OLS） 多项式回归 一般线性模型
廣義線性模式 离散选择（英语：Discrete choice） 对数几率回归 多项罗吉特（英语：Multinomial logit） 混合罗吉特 波比（英语：Probit model） 多项式波比（英语：Multinomial probit） 排序性模型（英语：Ordered logit） 有序波比（英语：Ordered probit） 泊松回归
等级线性模型 固定效应（英语：Fixed effects model） 随机效应（英语：Random effects model） 混合模型（英语：Mixed model）
非线性回归 非参数 半参数 稳健 分位数迴歸 保序回归 主成分 最小角 局部（英语：Local regression） 分段
含误差变量（英语：Errors-in-variables models）
估计
最小二乘法 普通最小二乘法 线性 偏最小二乘回归 总体（英语：Total least squares） 广义 加权 非线性 非负（英语：Non-negative least squares） 重复再加权（英语：Iteratively reweighted least squares） 脊迴歸（嶺迴歸） LASSO
最小绝对值导数法（英语：Least absolute deviations） 贝叶斯（英语：Bayesian linear regression） 贝叶斯多元
背景
回归模型驗證（英语：Regression model validation） 平均响应和预测响应（英语：Mean and predicted response） 误差和残差 拟合优度 学生化残差（英语：Studentized residual） 高斯-马尔可夫定理
概率与统计主题
查 论 编

在统计学上，泊松回归（英語：Poisson regression）是用来为计数资料（英语：Count data）和列联表建模的一种回归分析。泊松回归假设因变量（英语：response variable）Y是泊松分布，并假设它期望值的对数可由一组未知参数进行线性表达。当其用于列联表分析时，泊松回归模型也被称作对数-线性模型。

泊松回归模型是广义线性模型（GLM）的一种，以对数变化作为连接函数（link function），该模型的假设之一是其被解释变量服从泊松分布。

泊松回归模型

${\displaystyle \mathbf {x} \in \mathbb {R} ^{n))$ 代表由一组相互独立的变量组成的向量，其泊松回归的模型形式为:

${\displaystyle \log(\operatorname {E} (Y\mid \mathbf {x} ))=\alpha +\mathbf {\beta } '\mathbf {x} ,}$ ${\displaystyle \alpha \in \mathbb {R} }$， ${\displaystyle \mathbf {\beta } \in \mathbb {R} ^{n))$.

亦可简洁表示为：${\displaystyle \log(\operatorname {E} (Y\mid \mathbf {x} ))={\boldsymbol {\theta ))'\mathbf {x} ,\,}$

此处，${\displaystyle \mathbf {x} }$ 是 n+1维的向量，由n个独立变量（自变量向量）一个常向量（元素取值全为1）构成，用一个θ 代表第一个表达式当中的 α 和 β。

因此，当已知泊松回归模型当中的 θ和解释变量 ${\displaystyle \mathbf {x} }$, 其满足泊松分布的被解释变量的期望值可以由下式来预测：

${\displaystyle \operatorname {E} (Y\mid \mathbf {x} )=e^((\boldsymbol {\theta ))'\mathbf {x} }.\,}$

Y_i 是被解释变量的观测值，相应的解释变量为 x_i ，可由极大似然估计（MLE）的方法来估计参数θ。 极大似然估计不能通过解析表达式获得解析解，是由其对数似然函数为凸函数的特性，可通过Newton–Raphson或其他基于梯度下降的思想方法来进行参数估计。

极大似然估计

如上所述，已知泊松回归模型当中的 θ和解释变量 ${\displaystyle \mathbf {x} }$, 其回归表达式为：

${\displaystyle \operatorname {E} (Y\mid x)=e^{\theta 'x}\,}$,

泊松分布的概率密度函数为：

${\displaystyle p(y\mid x;\theta )={\frac {[\operatorname {E} (Y\mid x)]^{y}\times e^{-\operatorname {E} (Y\mid x))){y!))={\frac {e^{y\theta 'x}e^{-e^{\theta 'x))}{y!))}$

现已知解释变量的观测值为由 m个向量组成 ${\displaystyle x_{i}\in \mathbb {R} ^{n+1},\,i=1,\ldots ,m}$, 对应 m 个被解释变量的观测值，${\displaystyle y_{1},\ldots ,y_{m}\in \mathbb {R} }$. 若同时已知θ, 则该组观测值所对应的联合概率可由下式表达：

${\displaystyle p(y_{1},\ldots ,y_{m}\mid x_{1},\ldots ,x_{m};\theta )=\prod _{i=1}^{m}{\frac {e^{y_{i}\theta 'x_{i))e^{-e^{\theta 'x_{i)))){y_{i}!)).}$

极大似然方法估计 θ的核心思想是，去找到能使得基于当前观测值的联合概率尽可能达到最大的θ。（可理解为：变量的取值当前观测值，与取值为其他任何数值相比，是发生概率最高的事件）。 既然目标是寻找到最优的θ，可以先将上式的等号左边简单表达为关于θ 的表达式：

${\displaystyle L(\theta \mid X,Y)=\prod _{i=1}^{m}{\frac {e^{y_{i}\theta 'x_{i))e^{-e^{\theta 'x_{i)))){y_{i}!))}$.

注意等号右边的表达式并未改写，但通常难于付诸计算，因而采用其对数变化后的表达式（ log-likelihood）即：

${\displaystyle \ell (\theta \mid X,Y)=\log L(\theta \mid X,Y)=\sum _{i=1}^{m}\left(y_{i}\theta 'x_{i}-e^{\theta 'x_{i))-\log(y_{i}!)\right)}$.

由于 θ 仅出现在似然函数的前两项，因而在极大化似然函数的运算过程中，可以只考虑前两项。可以删去第三项y_i!，待优化的似然函数可以简洁表达为：

${\displaystyle \ell (\theta \mid X,Y)=\sum _{i=1}^{m}\left(y_{i}\theta 'x_{i}-e^{\theta 'x_{i))\right)}$.

为了找到极大值，需要求解方程：

${\displaystyle {\frac {\partial \ell (\theta \mid X,Y)}{\partial \theta ))=0}$

可以通过对其似然函数取负值 （negative log-likelihood）, ${\displaystyle -\ell (\theta \mid X,Y)}$是一个凸函数, 标准的凸优化方法可以考虑来求解 θ的最优值。统一的方法是Newton-Raphson 与Iterative Weighted Least Square（IWLS）算法。 给θ一组初始值，IWLS 是通过多次迭代更新直到θ 收敛。

泊松回归的应用

泊松回归常用于被解释变量为计数（Count）形式时，包括事件发生的次数，比如：客服中心接到的电话次数。其满足相互独立的假设。在此例子中，即为：拨打客服电话的人们之间不存在相互关联。不会因为甲拨打了客服，而影响乙拨打的可能性。但在建模时，需要考虑统计该事件发生的时期，比如目标变量统计的是一天接到的电话次数，还是一个星期，或者一个月。这个时期的数据作为回归模型中的抵消值，在下面解释。

"曝光量"（Exposure） 与 偏移量 (trade off)

泊松分布也可以适用于比率数据，即事件发生次数与其测量时间或测量范围的比值。比如生物学家测量某森林中树木种类的数目， 比率变量即为每平方千米的树木种类数。人口学家关注的是每个人口年（person-year）的人口死亡数。通常来说，比率变量表达的是单位时间内该事件发生的次数。这些例子中，平方米”，“人口年”这些变量就是所谓的"曝光量"（Exposure）。泊松回归中将其视为偏移量放在等式右边。

${\displaystyle \log {((\operatorname {E} (Y\mid x)/({\text{exposure))))}=\theta 'x}$

which implies

${\displaystyle \log {(\operatorname {E} (Y\mid x))}-\log {({\text{exposure)))}=\log {\left({\frac {\operatorname {E} (Y\mid x)}{\text{exposure))}\right)}=\theta 'x}$

在R中运行广义线性模型时，可用offset()来指定表示“曝光量”的变量:

glm(y ~ offset(log(exposure)) + x, family=poisson(link=log) )

过度离势和零膨胀

服从泊松分布的变量,具有期望与方差相等的特征。若观测样本的方差远大于期望值的时，则认为存在过度离势，当前的模型不合理。其常见的原因是缺失重要的解释变量。解决该问题的方法，通常采用准似然估计（quasi-likelihood） 或者负二项分布来估计。^[1][2]

泊松回归的另一个常见的问题是零膨胀zero-inflated model。标准的泊松分布其定义域为非负整数，被解释变量y取值为0的概率为：

${\displaystyle p(y=0\mid x;\theta )=e^{-e^{\theta 'x))}$

但如果观测样本中添加大量的0，则取值为0的频率远大于理论概率，此时不适宜直接采用泊松回归。比如观测一组人在一小时内的吸烟情况，目标变量是每人吸了多少根烟。但当观测人群中有大量的非吸烟者，就会有过多的目标变量为0， 这就是零膨胀。可以采用其他的广义线性模型，比如负二项分布负二项分布来建模，或者零膨胀模型zero-inflated model 来解决。

參考文獻

查 论 编 社會調查方法 資料蒐集（英语：Data collection） 蒐集方法（英语：Survey data collection） 人口普查 調查抽樣（英语：Survey sampling） 隨機抽樣 問卷調查 訪談（英语：Interview (research)） 結構式 半結構式（英语：Semi-structured interview） 無結構式（英语：Unstructured interview） 伴侶（英语：Couple interview） 非介入性研究 比較研究法（英语：Comparative research） 內容分析法 调查方法 人口普查 采样框架（英语：Sampling frame） 样本 调查抽样（英语：Survey sampling） 随机抽樣 简单随机抽样 Quota sampling 分层抽样 Nonprobability sampling 样本量确定 研究设计（英语：Research design） 纵向研究 队列研究 橫切面研究（英语：Cross-sectional study） Cross-sequential study 資料分析 類別資料 類別資料（英语：Contingency table） 測量尺度 描述統計學 探索式資料分析（英语：Exploratory data analysis） 多變量分析 心理統計學 推論統計學 機率模型 圖模式 泊松迴歸 結構（英语：Structural equation modeling） 應用 應用領域 收視率調查 人口學 市場調查 舆论调查 社會統計學 研究計畫 非洲晴雨表（英语：Afrobarometer） 美國選舉研究（英语：American National Election Studies） 亚洲晴雨表（英语：Asian Barometer Survey） 歐洲晴雨表 歐洲社會調查（英语：European Social Survey） 普通社會調查（英语：General Social Survey） 國際社會調查計畫 拉丁美洲晴雨表（英语：Latinobarómetro） 美國家戶調查列表（英语：List of household surveys in the United States） 國家健康與營養調查（英语：National Health and Nutrition Examination Survey） 紐西蘭社會態度與價值觀研究（英语：New Zealand Attitudes and Values Study） 世界價值觀調查（英语：World Values Survey） 學會 美國公共輿論研究學會 歐洲輿論與市場研究學會（英语：European Society for Opinion and Marketing Research） 蓋洛普 國際統計學會 皮尤研究中心 世界民意研究學會 分類