数据分享|R语言广义线性模型GLM：线性最小二乘、对数变换、泊松、二项式逻辑回归分析冰淇淋销售时间序列数据和模拟-1

https://developer.aliyun.com/article/1489380

泊松回归

计数数据的经典方法是泊松分布。

泊松分布只有一个参数，这里是 μi，这也是它的期望值。μi 的链接函数是对数，这意味着我必须将指数函数应用于线性模型才能恢复到原始比例。这是我的模型：

再说一遍，虽然我观察到的期望值是实数，但泊松分布只会产生整数，与实际销售额相符。

pos.md <- glm
display(poi.od)

Plot()
lines

这看起来很不错。系数的解释现在应该很清楚了。

从系数中我可以看出，0ºC 时，我预计会卖出 exp⁡(4.45)=94 冰淇淋，而温度每升高 1 度，预计销量会增加 exp⁡(0.076)−1=7.9% .

到现在为止还挺好。我的模型符合我的观察。此外，它不会预测负销售额，如果我使用上述模型给出的平均值从泊松分布进行模拟，我将始终只得到整数。

但是，我的模型还会预测，如果温度达到 32 摄氏度，我应该会卖出 1000 多个冰淇淋：

predict(pmod, newdata)S

二项式回归

好的，让我这样思考这个问题：我有 800 个潜在销售量，我想了解在给定温度下销售的比例。

这表明成功销售数量为 800 次的二项式分布。二项式分布的关键参数是成功概率，即有人购买我的冰淇淋的概率作为温度的函数。

因此，我需要一条将销售统计数据映射到 0 到 100% 之间的概率的 S 形曲线。

一个典型的选择是逻辑函数：

S

有了这个，我的模型可以描述为：

mize <- 1000
icectunity <- marksize - icenits
display(b.glm)

binred <- predict(biglm, type="response")*marsize
basicPlot

随着温度越来越高，该模型将预测销售将达到市场饱和，而迄今为止所有其他模型都将预测越来越高的销售。

我可以使用逻辑函数的倒数来预测 0ºC 和 35ºC 时的销售额：

# 0摄氏度下的销售
plogis(coef(biglm)\[1\])*market.size

# 在35摄氏度下销售
plogis(coef(bnm)\[1\] +  coef(bglm)\[2\]\*35)\*maksie

概括

让我们将所有模型放在一张图中，温度范围为 0 到 35ºC。

p.lm <- predict
po.lm <- exp + 
                 0.5 * sumary(loim)$dispersion)
p.pis <- preict(poiso daaframe(tp=tm, type="response")
p.bn <- predict(biglm, datafrme(emptem), type="espns")*arke.ze 
baPlot

该图表显示了我的四个模型在 0ºC 到 35ºC 的温度范围内的预测。尽管线性模型在 10ºC 到 30ºC 之间看起来还可以，但它清楚地表明了它的局限性。对数变换的线性模型和泊松模型似乎给出了类似的预测，但将预测随着温度的升高，销售额将不断加速增长。我不相信这是有道理的，因为即使是最喜欢冰淇淋的人也只能在非常炎热的一天吃这么多冰淇淋。这就是为什么我会使用二项式模型来预测冰淇淋销量。

模拟

使用以分布为中心的视图来描述我的模型自然会导致模拟。如果模型很好，那么我应该无法从模拟中识别出真实数据。

在我所有的模型中，线性结构都是

或以矩阵表示法

其中 Ai,⋅=[1,xi] 和 v=[α,β]，其中 A 是模型矩阵，v 是系数向量。

话虽如此，让我们模拟原始数据中测量的温度的每个分布的数据，并与实际销售单位进行比较。

n <- nrow(icre)
A <- modl.(uits ~ temp, data=cam)
set.seed(1234)
(rad.nal <- rnorm(n,
                     mean  A %*% cof(li.od),
                     sd = sqrt(sumary(liod)$esion)))

(ranlans <- rlnorm(n,
                         mnog = A %*% coef(.od),
                         sdlog =  sqrt(summary(loiod)$isin)))

(nd.ps <- rpois(n,
                   labd = exp(A %*% coef(piod))))

(ra<- rbinom(n,
                   size = meze,
                   prob = plogis(A %*% coef(b.m))))

bacPlot
cols <- adscor(c
points(iceram$tmp,  pch=19, col=cols\[1\])

该图表仅显示每个模型的一个模拟，但显示了一些有趣的方面。我不仅看到泊松和二项式模型生成整数，而高斯和对数变换的高斯预测实数，我注意到红点处对数正态分布的偏度为 19.4ºC。

此外，线性模型预测高于和低于平均值的可能性相同，在 16.4ºC 时，预测似乎有点低 - 可能是结果。

此外，对数转换和泊松模型在 25.1ºC 时的高销售额预测也不意外。

同样，二项式模型的模拟似乎是最接近现实的。

结论

我希望这篇文章能说明广义线性模型背后的直觉。

将模型拟合到数据需要的不仅仅是应用算法。特别值得思考的是：

• 期望值的范围：它们是有界的还是范围从 -∞ 到 ∞？
  
• 观察类型：我期望实数、整数还是比例？
  
• 如何将分布参数与观测值联系起来