数据分享

您所在的位置:网站首页 残差图怎么画 数据分享

数据分享

2023-03-22 06:59| 来源: 网络整理| 查看: 265

原文链接:http://tecdat.cn/?p=27464 

在拟合 GLM(并检查残差)之后,可以使用 z 检验一一检验估计参数的显着性,即将估计值与其标准误差进行比较。

相关视频

GLM 模型拟合和分析示例 示例 1. 小鼠数据的 GLM 建模(剂量和反应)

a) 我们输入数据(查看文末了解数据获取方式)并拟合逻辑回归模型。

9306553154246da4c7004f3dd160c5f9.png

63da70894df29cc523d72817a89efc11.png

61e91dce67175f60b9477f6f602a6d67.png

> summary(it1.lt)

5394e0327ff1861ffe34a82e24120364.png

1-pchisq(17.6,24)

0a36b4986ede051bf50f16fd76f5a4d3.png

模型:0db92035bae9b575ad41795cc26e512a.png

可以与完整模型进行比较。与偏差值 17.639 相关的 P 值 0.82(> 0.10)意味着没有显着证据拒绝拟合模型。

anova(fi.lgi)

2a0c902b0c0cb03b5d50b3647377a16b.png

5fc100dcb984c631913c1573a0b9a07a.png

1-pchisq(35.8-17.69, 25-24)

314a6f54fb0de5a04f21e18a199b9b31.png

空模型不包含预测变量,在 25 个自由度 (df) 上的偏差为 35.89。当协变量 x 添加到空模型时,偏差的变化是 35.890-17.639=18.25。与自由度为 25-24=1 的卡方分布相比,其 P 值为 1.93 × 10 -5 非常显着。

因此模型不能通过删除 x 来简化。x 的系数的 t 检验也很重要(P 值 0.0065 plotx, itte(fi1log,typ"> pot(,y

ee91e72d51e342575f42923a92619431.png

图 1:逻辑回归的小鼠数据和拟合值。

点击标题查阅往期内容

73a3e6fe92ff9f45ea22c46956e9dc52.jpeg

R语言贝叶斯广义线性混合(多层次/水平/嵌套)模型GLMM、逻辑回归分析教育留级影响因素数据

outside_default.png

左右滑动查看更多

outside_default.png

01

a0d11ec6bc350d437745de9fcf3c55de.jpeg

02

496fe4ff07ad370235aae3423f1d9666.png

03

7cbcfcb142030f967d801f632a45fd18.jpeg

04

28d876e9ed45808ed4bc3516375b705f.jpeg

b)我们拟合一个带有概率链接的模型。a430caa575ef279f85c19b9aec3c4601.png

7403772c9970826660a1699f804e60a3.png

> summary()

6aea3aa47d055f66863ebbe232f8e3df.png

配套模型:

e1bce76afe091c43c6b776084ff2035b.png

同样,这两个参数都很重要(P 值 anova

4ce5bd61de4aaf474f285aed1a219169.png

> 1 - pchisq(35.89-17.49 25-24)

bef4884df32154c5a7f90356ecfb2b87.png

> lines(x, fitte

添加 x 时偏差的变化是显着的(P 值 =c4e366146701d0bac9cd25f47dd6743a.png)。模型不能通过删除 x 来简化。

73431f4cfcc96ba5e0fd997e085a8a14.png

图 2:小鼠数据和拟合值(虚线:概率链接)。使用 probit 链接的模型略好于使用 logit 链接的模型,因为偏差更小。在两个模型中,x 的系数都很显着(P 值 summary(it2.it)

7b0ae69d88b8dcec8fd1891dd1966ba9.png

0f05fc31b226c075452f02e7098329d1.png

1-pchisq(13.63,6)

c3cb919a8c5a80713c1c84b99558f5c7.png

与偏差值 13.633 相关的 P 值为 0.034 plot(x, reid(it2.it))

91e1f291c4c45bfcfc75e07e23ac4aa8.png

图 3:仅带有 x 的拟合模型的残差图。

plt(fitdft2lit reid(fi2lot))

所以我们将 x2 添加到模型中。

> summary(ft2qlt)

4583a0b819ffe4f413051e8ee5f4bf0c.png

30d2d35350a7214727806838d6818a6d.png

> 1-pchisq(5.1, 5)

e5a37bd3261474ed4efab63a2c47076a.png

偏差从 13.633 减少到 5.107,不显着(P 值=0.403>0.05)。因此,我们不能通过偏差的证据来拒绝这个模型。

plot(fitted(fit2logit), resid(.logit))

f95424e6a6159a24d835ca1f612a7f9a.png

图 4:带有 x2 的拟合模型的残差图。

残差现在看起来是随机的。拟合模型为a87bd35b661f5da2d8c90e045c38d279.png

并且所有参数估计值都很显着(5%)。对数几率06d9284c610b448f4487638dce220e98.png以二次方式依赖于 x。

 示例 3. 艾滋病数据,泊松

a) 我们输入数据并使用默认对数链接拟合泊松回归模型。

dea1e853a1ab67b27f9a123407b1a7d9.png

726382b9c824c600300c4d1bc4719928.png

> smary(fit.lg)

d965ba5273a203226ca00ba018d936ca.png

> 1-pchisq

065e707c7f0be0ed5ff6a39f1ad630f8.png

偏差 29.654 的 P 值为 0.005 summary(fi3.lg)

8a402b05a65b6a8544ceb4935e0a4f18.png

> 1-pchisq \[1\] 0.1279

16.371 的偏差(P 值为 0.1279>0.10)并不显着。拟合模型e80006f9e2d2f6766fccc7724d6d7055.png

不能以偏差为由拒绝。但残差图只显示了比以前稍微随机的模式。 9dc124286b4b251eda0783efe40c1e60.png

图 6:拟合模型 3 与 x2 的残差图

该模型可以通过使用非规范链接进行改进。

> summary(ft3st)

b5a9b2cda4dee98748012ab3971853ce.png

> 1 - pchisq(16.9, 12)\[1\] 0.153

拟合模型的残差 y = (-0.27571 +0.49277x)2 + e 显示出更加随机的模式。12 df 上的偏差 16.905 略高于之前模型的 16.371(df=11),但仍然不显着(P 值=0.1532>0.10)。AIC 较小,为 73.833



【本文地址】


今日新闻


推荐新闻


CopyRight 2018-2019 办公设备维修网 版权所有 豫ICP备15022753号-3