为什么线性回归对残差有假设,而广义线性模型对响应有假设?


14

为什么线性回归和广义模型的假设不一致?

  • 在线性回归中,我们假设残差来自高斯
  • 在其他回归(逻辑回归,毒物回归)中,我们假设响应来自某种分布(二项式,泊松等)。

为什么有时会假设剩余而其他时间会在响应时?是因为我们要导出不同的属性?


编辑:我认为mark999的显示两种形式是相等的。但是,我对iid还有其他疑问:

我的另一个问题 是,逻辑回归是否有iid假设?显示广义线性模型没有iid假设(独立但不相同)

对于线性回归,是否真的成立,如果我们对残差进行假设,我们将有iid,但是如果对响应进行假设,我们将拥有独立但不相同的样本(具有不同不同高斯样本)?μ


Answers:


12

具有高斯误差的简单线性回归是一个很好的属性,它不能推广到广义线性模型。

在广义线性模型中,响应遵循给出均值的给定分布。线性回归遵循这种模式。如果我们有

yi=β0+β1xi+ϵi

ϵiN(0,σ)

那我们也有

yiN(β0+β1xi,σ)

好的,因此响应遵循广义线性模型的给定分布,但是对于线性回归,我们必须使残差遵循高斯分布。为什么不是通用规则时强调残差是正常的?好吧,因为这是更有用的规则。考虑残差的正态性的好处是,这很容易检查。如果我们减去估计的均值,则所有残差都应具有大致相同的方差和大致相同的均值(0),并且将大致呈正态分布(请注意:我说“大致”是因为,如果我们没有完美的估计值,回归参数,我们当然不这么做, Xϵix。但希望估算中有足够的精度可以忽略不计!)。

另一方面,查看未调整的,如果它们均具有不同的平均值,我们就无法真正判断它们是否正常。例如,考虑以下模型:yi

yi=0+2×xi+ϵi

与和ϵiN(0,0.2)xiBernoulli(p=0.5)

则将是高度双峰的,但不会违反线性回归的假设!另一方面,残差将遵循大致正态分布。yi

这是一些R代码来说明。

x <- rbinom(1000, size = 1, prob = 0.5)
y <- 2 * x + rnorm(1000, sd = 0.2)
fit <- lm(y ~ x)
resids <- residuals(fit)
par(mfrow = c(1,2))
hist(y, main = 'Distribution of Responses')
hist(resids, main = 'Distribution of Residuals')

直方图


yi=1+2×xi+ϵi

3
@ hxd1011:是的,这是边际分布(显然不是正态)与给定x的条件分布之间的差(自从我们模拟它以来,我们就知道它是正态的!)。不考虑条件分布和边际分布之间的差异是一个非常普遍的错误。
悬崖AB

14

i=1,,n

Yi=β0+β1Xi1++βkXik+ϵi,
ϵiσ2Xi1,,XikYiβ0+β1Xi1++βkXikσ2

Xi1,,Xikβ0+β1Xi1++βkXik

通常的具有正态误差的多元线性回归模型是具有正态响应和标识链接的广义线性模型。

By using our site, you acknowledge that you have read and understand our Cookie Policy and Privacy Policy.
Licensed under cc by-sa 3.0 with attribution required.