Tobit生存模型中T的预测疑问:是否需考虑对数正态分布期望?
exp(xβ̂)不是最优选择? 嘿,你提的这个点戳中了对数线性模型预测的一个常见误区——直接用线性预测值的指数exp(xβ̂)作为T的估计值确实有问题,而且必须考虑exp(e)的对数正态期望特性,下面给你掰明白:
核心逻辑:对数正态变量的期望不等于1
根据你的模型设定:
$$T = \exp(x\hat{\beta}) \times \exp(e)$$
其中$e \sim iid\ N(0,\hat{\sigma}2)$,那$\exp(e)$服从对数正态分布$LN(0,\hat{\sigma}2)$。对数正态分布的期望有个固定公式:
$$E[\exp(e)] = \exp\left(\frac{\hat{\sigma}^2}{2}\right)$$
只要模型存在随机误差(也就是$\hat{\sigma}^2>0$,这在实证研究里几乎是必然的),这个值就大于1,而且误差越大,这个值离1越远。
常规做法的偏差来源
直接用exp(x\hat{\beta})预测T,本质上是默认了$E[\exp(e)]=1$,但实际上T的真实期望是:
$$E[T] = \exp(x\hat{\beta}) \times E[\exp(e)] = \exp\left(x\hat{\beta} + \frac{\hat{\sigma}^2}{2}\right)$$
所以只取exp(xβ̂)会系统性低估T的期望,偏差程度和$\hat{\sigma}^2$正相关——模型拟合得越差(误差越大),低估得越离谱。
该怎么修正?
如果你的目标是预测T的期望,那正确的估计式应该是:
$$\hat{T} = \exp\left(x\hat{\beta} + \frac{\hat{\sigma}^2}{2}\right)$$
不过这里要区分清楚:如果你要的是中位数预测,那exp(xβ̂)是完全没问题的——因为对数正态分布的中位数是1(对应e的中位数是0),此时exp(xβ̂)正好是T的中位数。所以关键看你预测的是期望还是中位数。
补充:关于Tobit的截尾特性
哦对了,你提到的是Tobit模型,不过你的设定是$\ln(T)=xβ+e$,如果这是无截尾的对数线性模型(有时候也被宽泛地称为对数Tobit),上面的结论完全适用。如果是有截尾的情况(比如T有下界,比如T≥0),那还要额外考虑截尾对期望的修正,但你的问题里没提截尾,所以先聚焦对数正态的期望偏差问题。
内容的提问来源于stack exchange,提问作者Phil Viton

