r - 求解期望值泊松回归

Question

我有每对标记为 Sij 的大脑区域之间在大脑中运行的轨道数量。我还有每对标记为 gij 的区域之间的轨道距离。i 和 j 代表每个大脑区域。

因此，例如这个 if 文件：

Sij    gij
331     15.2
428     11.1
797     45
313     54
142     12

我正在尝试使用泊松回归模型调整低于 12 距离的 gij 值的偏差。

我想要做的是解决得到 alpha0 和 alpha1 给定这个具有日志链接功能的泊松模型：

log(μ(Sij|gij))=α0+α1gij

但是这里的问题是我不确定如何找到这个值：μ(Sij|gij)，它等效地称为期望值 E(x)。

我正在考虑这样做是r代码：

summary(m1 <- glm(Sij$file ~ gij$file, family=poisson(link=log), data=p))

但我的理解是，从中我会得到 alpha0 和 alpha1 但我不确定如何获得预期值

标签： rloggingregressionglmpoisson

Answer 1

对于期望值，有一个命令，它是加权平均值：

weighted.mean(Sij, gij)