原文链接：http://tecdat.cn/?p=24127
最近我们被客户要求撰写关于预测UCI鲍鱼年龄数据的研究报告，包括一些图形和统计输出。
鲍鱼是一种贝类，在世界许多地方都被视为美味佳肴
介绍
养殖者通常会切开贝壳并通过显微镜计算环数来估计鲍鱼的年龄。因此，判断鲍鱼的年龄很困难，主要是因为它们的大小不仅取决于它们的年龄，还取决于食物的供应情况。而且，鲍鱼有时会形成所谓的“发育不良”种群，其生长特征与其他鲍鱼种群非常不同。这种复杂的方法增加了成本并限制了其普及。我们在这份报告中的目标是找出最好的指标来预测鲍鱼的环，然后是鲍鱼的年龄。
数据集
背景介绍
这个数据集 （ 查看文末了解数据获取方式 ） 来自一项原始（非机器学习）研究。
从原始数据中删除了有缺失值的例子（大多数预测值缺失），连续值的范围被缩放用于NA（通过除以200）。在本分析中，我们将通过乘以200的方式将这些变量恢复到其原始形式。
数据集中的观测值总数：4176
数据集中的变量总数：8个
变量列表
变量 | 数据类型 | 测量 | 描述 |
|---|---|---|---|
性别 | 分类（因子） | M、F 和 I（婴儿） | |
长度 | 连续 | 毫米 | 最长壳测量 |
直径 | 连续 | 毫米 | 垂直长度 |
高度 | 连续 | 毫米 | 带壳肉 |
整体重量 | 连续 | 克 | 整只鲍鱼 |
去壳重量 | 连续 | 克 | 肉的重量 |
内脏重量 | 连续 | 克 | 肠道重量 |
外壳重量 | 连续 | 克 | 晒干后 |
鲍鱼的环 | 连续 | 1.5 给出以年为单位的年龄 |
下面是分析
“使用回归预测鲍鱼的年龄”
数据汇总与统计
代码语言：javascript复制balne$Sx <- s.acor(aalne$Sex)
kale(abaoe[1:10,],fomt 'madw')
分类变量
数值变量
看一下数据集的摘要，我们可以看到，数据在雄性、雌性和婴儿这三个因素水平之间的分布是相当均匀的。
因变量
因果变量Rings包含在数据集中。它被测量为切割和检查鲍鱼后观察到的环的数量。虽然它不能直接表示一个给定的鲍鱼的年龄，但它可以或多或少完美地确定它。一个鲍鱼的年龄等于环数 1.5。由于这种关系是可靠的，环数将被视为因变量。数据中测量的环数从1到29不等，大多数鲍鱼的环数在5到15之间。分布也有轻微的正偏斜，但没有问题。(见下面的图)
配对图
代码语言：javascript复制pairs(aalone, es(colour =Sex, aph = 0.)
从配对图中观察到的情况。
首先要注意的是数据的高度相关性。例如，直径和长度之间的相关性非常高（约98.7）。
r语言中对LASSO回归，Ridge岭回归和弹性网络Elastic Net模型实现
01
02
03
04
同样，Whole_weight似乎与其他重量预测因子高度相关，是Shucked_weight、Viscera_weight和Shell_weight之和。
其次，预测因子Sex的分布与所有其他预测因子的因子水平值雌性和雄性非常相似。
对于雌性和雄性的因子水平，分布的形状也是非常相似的。
我们可以考虑重新定义这一特征，将性别定义为婴儿与非婴儿（其中非婴儿=雌性和雄性都是）。
大多数的abalones环都在5到15之间。
数据质量
增加变量。我们将更新鲍鱼数据集，创建名为 "婴儿 "的新变量，它的值将基于性别变量的原始值。当性别变量为I时，它的值为I，否则为NI。 我们还观察到，预测高度的最小值是0，实际上这是不可能的，我们将调查这些观察结果，仔细研究。
代码语言：javascript复制##高度为0的数据质量检查
kable(abloe[aban$Height == 0,])我们看到，有两个观测值的高度可能没有被正确记录，因为其他预测因子似乎都有有效的值。另外，如果我们看一下预测因子Whole_weight，我们会发现这些值与其他观察值相比真的很小，而且低于第一个四分法。这告诉我们，这可能不是一个数据错误，因此我们不能将这些数据从我们的数据集中排除。
我们还将添加一个名为weight.diff的新变量。我们可以在摘要中看到有四种不同的重量测量方法，即Whole_weight、Shucked_weight、Viscera_weight和Shell.weight。Whole_weight是其他重量预测因子的线性函数，在剥壳过程中损失的水/血的质量未知。
代码语言：javascript复制str(aane, give.attr= FASE)
我们看到变量Whole_weight应该是Shucked_weight、Viscersa_weight和Shell_weight的线性函数，我们可以写成Whole_weight = Shucked_weight Viscera_weight Shell_weight 剥壳过程中损失的未知水/血质量。
然而，当我们计算Whole_weight和其他重量变量之间的差异时，我们发现有153个观测值违反了这一规定，也就是说，这似乎不符合逻辑，可能是记录数据时的错误。
如果我们绘制新添加的weight.diff变量的直方图，我们可以看到，当weight.diff为负数时，有一些观察结果。
代码语言：javascript复制#确定没有正确记录的观察结果
#不符合逻辑的观察结果的柱状图
ggplt(aalone, as(x=weight.diff)) ,ill=rb(1,.4,0,.7), bins = 30)
我们来看看其中的一些记录。
代码语言：javascript复制nrow(ablon[abaoneweihtdff < 0,])
请注意，总共有153个观测值的综合权重超过了Whole_weight。当我们看了10个这样的观测值时，似乎其他的值都是正确的，没有任何相似之处，所以我们确信这可能是一个数据输入错误。因此，我们将保留这些观察结果，以便进一步分析。
我们将首先在训练和测试中潜入我们的数据集。数据集将以70/30的比例在训练和测试之间进行分割，并随机选择观测值。
训练和测试拆分
代码语言：javascript复制set.ee(4)
#使用70/30方法在训练和测试中分割数据
ndxes <-spl(1:owabaone, size= 0.3 nrw(bone))
aboetrai <- ablon[-indxs,]
abetest <- abloneindxe,]我们已经开始用所有的变量拟合一个加法模型，并将研究参数的重要性。在此基础上，我们将修改我们的模型。现在我们将使用变量Sex的原始值，它的因子水平为F、I和M。
加性多元线性回归模型
代码语言：javascript复制summary(abneadd)
在第一个加性模型中，注意因子水平雌性是性别变量的参考水平。
在用所有预测因子拟合加性模型后，我们可以看到，除了长度之外，测试统计显示所有变量都是显著的。正如我们之前从配对图中看到的那样，长度和直径的预测因子是高度相关的。我们还看到，不同重量的预测因子也是显著的，尽管它们应该是彼此的线性函数。
RMSE 分数
代码语言：javascript复制kable(rmse(aaloe_ad,"Aditve odel"))
我们将计算方差膨胀因子，以发现数据集存在的多重共线性问题。
多重共线性
代码语言：javascript复制vif
我们看了所有变量的变量膨胀系数，似乎所有的预测因子都有多重共线性问题，除了我们之前在配对图中看到的性别和身高。预测因子Whole_weight的VIF值最高，因为它是其他体重的线性函数。
Whole_weight & Rings 之间的偏相关系数：我们将首先计算Whole_weight 变量和因变量（Rings）的偏相关系数 。
#检查高共线性关系变量的异方差性
wole_wigt_it <- lm(holweight ~Sx LnhDametr eit Sucked_ght Visrwght Shl_wegh data=alotrin)变量添加图
同样地，变量添加图将这些残差相互之间的关系可视化。将因变量的残差与预测的残差进行回归，并将回归线添加到图中，也是有帮助的。
代码语言：javascript复制cor(resid(whole_weight),resid(addwtouwolwigh))
cre_plot(baead_itht_whe_eght,wleeghtfit)
没有Whole_weight的加法模型的方差膨胀因子
但直径和长度的VIF还是很高。
代码语言：javascript复制vif(abaln_ddithu_whoeeiht)
直径和环之间的偏相关系数
我们现在将 在模型中Diameter 没有Whole_weight变量的情况下计算变量和因变量（环） 的偏相关系数 。
mete_i <- lm(Diameter ~ Sex Length Height Shucked_weight Viscera_weight Shell_weight)
abaoned_sal <- lm(Rings ~ Sex Length Height Shucked_weight Viscera_weight Shell_weight)这两个残差的相关性接近于零，这意味着未被性别、长度、高度、去壳重量、内脏重量和贝壳重量解释的环的变化与未被性别、长度、高度、去壳重量、内脏重量和贝壳重量解释的直径的变化的相关性很小。因此，在模型中加入直径可能没有什么好处。
代码语言：javascript复制cor(resid(damer_it),resid(abonead_mll))
creaevarlt(ablone_d_smaldiaete_fi)
没有 Whole_weight & Diameter 的加法模型的方差膨胀因子
代码语言：javascript复制vif(ablonadd_mll)
现在的VIF要低得多。我们将同时使用abalone_add和abalone_add_small进行分析。
abalone_add_small 的 RMSE 分数
代码语言：javascript复制kable(rmse(abalone_add_small
方差分析 F 检验
代码语言：javascript复制anova(small,ablone_add)
根据Anova检验，我们可以拒绝无效假设。
现在让我们尝试对加性模型进行AIC和BIC的参数选择。
代码语言：javascript复制#在加性模型上运行AIC和BIC
step(abae_add, dtonbackward" trace0)
step(abalon_ad, direction="backward"
我们使用之前Anova F测试中的最佳模型（abalone_add）运行了AIC和BIC方法，AIC和BIC都选择了没有Length预测器的相同模型。
我们选择了BIC的模型（因为两个模型都是一样的），并将绘制拟合与残差和QQ正态图。
加性模型假设
代码语言：javascript复制assumptions(aln_odad_bic,"baln_meddbic")这里我们看到，拟合图与残差图表明，模型违反了恒定方差，似乎有一种模式，即残差随着拟合值的增加而增加。
正态QQ图也有肥尾，表明误差可能不是正态分布。
表明我们可能需要寻找改进这个模型。尽管多重共线性对预测没有影响，但这看起来是一个极端的案例，有一个明显的模式违反了模型的假设。这可以通过几种技术来实现，例如变量选择和转换。让我们来看看变量选择方法，看看哪些变量可以用来创建另一个不违反恒定方差和正态性假设的模型。
穷举搜索
代码语言：javascript复制#穷举搜索
allabaone_add<- sumr(ruetsings Sex Legth Diamter Hight Whole_eght Shllweigh Shucke_weght Viscea_weigh , data=ablontra))
for(i in c(1:8)){
vr_nm[i]=sum(all_abwh[i,])-1
}
plot(var_num,all_a)
(besr <- which.max(adjr2))
alabaoe_ad$hch[bsj2,]
#画出模型参数与AIC的关系图
n * log(a_aln_dd$rs / n) 2 * (2:p)
plot(aloe_mo_ac ~ I(2:), ylab = "AIC"
")
现在我们看到了一些有趣的结果。之前我们看到t检验显示一些预测因子是不显著的，但是当我们进行穷举搜索时，它表明我们确实需要所有的预测因子来创建AIC值最低的模型。从图中可以看出，AIC值随着8个参数的模型大小而下降，并且是最小的。我们将再次使用数据集中的所有预测因子来创建模型，并寻找变量转换技术。
接下来，为了稳定恒定的变化，我们将进行一些因变量和预测变量的转换。
因变量转换
Box-Cox 变换
稳定方差的方法之一是使用对数转换因变量。为了得到正确的顺序，我们使用了boxcox方法，该方法建议使用$0的值。因为在0的值上，对数可能性最大，而且区间非常接近。因此，我们将使用log(Rings)形式的转换，用于我们的加性模型。
代码语言：javascript复制boxcox(abloe_ad lambda = seq(-0.1, 0.1, by = 0.1))
Additive 模型与对数因变量转换
代码语言：javascript复制summary(abaone_dd_log)
将因变量进行对数转换后，我们看到t检验是显著的，它也增加了先前加法模型的调整r平方值。我们还看到，在这个模型中，几乎所有的预测因子都是显著的。让我们检查一下假设。
模型假设
下面的拟合与残差图和Q-Q图显示，对因变量进行对数转换后，结果有了很大的改善。
代码语言：javascript复制assumptionsba
均方根分数
代码语言：javascript复制kable(log_rmse(abalo)
然而，我们没有看到RMSE分数有任何改善。恒定方差问题似乎得到了改善，QQ图也看起来不错。
下一步，我们将对预测器进行一些转换，并评估模型，看看这是否有助于进一步提高预测的准确性。
Predictor 转换
回归分析
为了使我们能够进行任何预测器的转换，首先让我们看看每个预测变量和因变量的关系。转换将取决于数据的形状以及预测因子和因变量之间的关系。
代码语言：javascript复制scatter(abale_tra$Lngt,abaone_train$Rngs,"Lenth""Rngs"),
我们可以看到环和预测指标长度、直径、高度的关系几乎是线性的。我们还可以看到，重量预测指标之间的关系并不是真正的线性关系，而是可以从多项式转换中受益。因此，让我们使用高阶多项式创建一个模型，即所有重量预测指标Whole_weight、Viscera_weight、Shucked_weight和Shell_weight。
多项式
在模型中使用二阶项后，模型假设相同。
代码语言：javascript复制asumptons(abloe_dd_oly2,"Poly2 Log Model")
均方根分数
代码语言：javascript复制kable(log_rmse(abaoly2,"Poly2 Log Moel)
均方根分数
在这里，我们进行了一些变量转换。首先，我们按照Boxcox方法的建议对因变量进行了对数转换，并按照对数图的建议对权重预测因子进行了多项式转换。在拟合模型后，我们看到rmse比以前的模型要低，与以前拟合的加性模型相比，它也有更好的恒定方差和Q-Q图。由于我们已经进行了程度为2的多项式转换，让我们尝试拟合程度为3的另一个模型并检查其意义。
方差分析 F 检验
代码语言：javascript复制anova(abaloe_addpoy2,aalon_add_oy3)
均方根分数
代码语言：javascript复制kable(log_rmse(abaloe_dd_pol4
方差分析 F 检验
代码语言：javascript复制anova
均方根分数
代码语言：javascript复制kable(log_rmse(abloneaddpoly5
方差分析 F 检验
代码语言：javascript复制anova
- 我们再次看到测试对于较低的 rmse 是显着的。让我们尝试拟合度数为 6 的模型。
均方根分数
代码语言：javascript复制kable(log_rmseaban_dd_poly6
方差分析 F 检验
代码语言：javascript复制anova
现在在用多项式次数为 6 进行拟合后，我们看到即使 F 检验表明它很重要，但检验的 RMSE 上升了。这表明我们现在可能已经开始过度拟合数据，即我们的模型非常接近地拟合数据，这是我们不希望发生的。
在此之前，我们看到多项式次数为 5 和 4 的测试和训练 RMSE 之间存在非常细微的差异。测试 RMSE 几乎相同。因此，我们愿意牺牲相对于更简单模型的 RMSE 非常微小的改进（第三个小数点）。因此我们选择多项式次数为 4 的模型，即模型 abalone_add_poly4。
代码语言：javascript复制for(d in um_poly){
abalone_add_polyestmodel(d)
rmse=g_log_mse(balone_ad_poly)
train_rmse[d]rmse$tran
test_re[d]=rse$st
}
plot(train_rmse
我们看到多项式次数为 5 和 4 的测试和训练 RMSE 之间存在非常细微的差异。测试 RMSE 几乎相同。因此，我们愿意牺牲相对于更简单模型的 RMSE 非常微小的改进。因此我们选择多项式次数为 4 的模型，即模型 abalone_add_poly4。
既然我们已经选择了模型，让运行 AIC 和 BIC 方法进一步选择合适的模型，看看我们是否可以做进一步的改进。
现在让我们计算和比较高阶项的 RMSE，并绘制训练和测试数据的均方根误差。
多加法模型上的 AIC 和 BIC：
- 既然我们已经选择了模型，让我们运行
AIC和BIC方法来进一步选择合适的模型。
step(abane_ad_poy4, directin="backwrd", trac=FALSE)Compare AIC 与 BIC 模型参数
代码语言：javascript复制_aic$call[2]
add_bic$call[2]
Anove F 检验
代码语言：javascript复制anova(abalone_mode
- 选择的模型
BIC中没有预测器Length。Anova F 检验的 p 值很大，因此我们无法拒绝原假设。abalone_model_add_bic模型很重要，因此我们将继续推进并检查模型假设。
模型假设（AIC 和 BIC）：
代码语言：javascript复制model_assumption- 在这种情况下，恒定方差和正态性看起来都不错。
RMSE 分数 - AIC
代码语言：javascript复制kable(log_rmse(abaone_mde_down')
RMSE 分数 - BIC
代码语言：javascript复制kable(log_rmse(abalone_model_add_bic,paste("Additive Model - Degree 4 - BIC")), digits = 4,format = 'markdown')
在使用AIC 和 进行变量选择后 BIC，我们从中选择了模型 BIC 并检查了 t 统计量和假设。有趣的是， BIC 模型丢弃了很少的预测变量，但也具有与我们开始使用的原始模型（多项式次数为 4 的模型）相似的测试 RMSE。这表明我们可以删除一些变量并仍然保持较低的 RMSE。这将我们带到下一个修改和引入BIC 上述模型选择的变量之间的交互项 。
来自BIC 模型的模型假设 也看起来更好。
接下来，我们将介绍交互项，并将尝试使用BIC 方法建议的预测变量来拟合模型 。
交互模型
代码语言：javascript复制log(Rings) ~ Height Diameter poly(Whole_weight, 4)
poly(Viscera_weight, 4) poly(Shucked_weight,4) poly(Shell_weight, 4) Sex Diameter:poly(Shucked_weight, 4) poly(Shucked_weight, 4):SexRMSE 分数
方差分析 F 检验
代码语言：javascript复制anova
在拟合交互模型并使用最佳可加模型执行 F 检验后，我们看到该检验表明交互模型是一个显着模型，具有改进的调整 r 平方值。RMSE 也变低了 因为它更好地解释了可变性，我们现在将选择交互模型并尝试在交互模型上运行 AIC 和 BIC。
同时，我们会比较交互模型的多个度数，以计算和比较高阶项的 RMSE，并绘制训练和测试数据的均方根误差。
代码语言：javascript复制for(d in num_poly){
ablone_int_poly=test_itmodel(d)
rmse=et_lg_rmseaaloneint_poly)
trainrse_int[d]=rmse$train
test_mse_it[d]=rme$tst
}
plot(tran_rse_n
- 我们可以看到，随着多项式次数的增加，RMSE 越来越低。尽管对于此分析而言，RMSE 的这种改进非常微小，我们可以忽略这种对模型简单性的改进。考虑到这一点，我们可以看到多项式次数为 4 的模型性能更好，因此我们将继续使用该模型。
交互模型上的 AIC 和 BIC
代码语言：javascript复制step(aalone_int, diretin="backar", trac=FALSE)
step(aalone_nt, diection="bacward", =loce=ALE)
RMSE 分数 - BIC/AIC
代码语言：javascript复制kable(log_rmse(abalone_model_int_bic,paste("Interaction Model - Degree 4 - BIC")), digits = 4,format = 'markdown')
在我们的交互模型上运行 AIC 和 BIC 后，我们看到该模型选择了相同的模型。由于这是我们迄今为止看到的最好的模型之一，具有合理的复杂性，我们将把它视为我们比较的候选模型之一，作为本分析的最佳拟合模型。
谈到候选模型，在详尽的搜索过程中，我们已经看到，当我们使用所有预测变量时，模型附带了最低的 AIC。我们可以尝试构建一个模型，其中包含所有具有交互作用和多项式次数的预测变量，并与我们选择的第一个候选模型进行比较，看看它的表现如何。因此，让我们拟合一个包含所有预测变量的模型。
在最初的数据分析中，我们发现Sex 因子水平为 female 和 的分类变量的分布 male 极其相似。因此，我们决定将这两个因子水平合并为一个，并且总因子水平为 2 infant 和 non-infant。我们创建了新变量 Infant。这里 non-infant 代表 female 和 male 两者。我们也通过这种方法进行了分析（可以在本报告的附录部分找到）。
让我们看看 Infant 模型分析，看看这个模型如何与我们上面选择的模型相抗衡。
婴儿模型分析
我们讨论过针对此分析采用不同的方法。我们引入了一个新的分类预测变量名称 Infant。我们使用现有的Sex 具有 3 个因子水平的分类预测变量，并创建了一个具有 2 个因子水平的新分类预测变量 。我们这样做是因为我们从原始分类预测变量female 和 中 确定了 2 个因子水平上的相似分布 male。新的因素水平现在是 I （婴儿 = 雌性和雄性组合）和 NI （非婴儿）。
这个新分类的分析与上面的分析完全一样，所以我们将用最少的解释和细节快速进行这个分析。
代码语言：javascript复制summary(abae_d_nf)
均方根分数
代码语言：javascript复制rmse <- funcin(atual predicted) {
sqrt(mean((actual - predicted ^ 2))
}
加性模型假设
代码语言：javascript复制model_assumptionBox-Cox 变换
代码语言：javascript复制boxcox(abon_adinf,lmda seq(-0.1, 0.1, by = 0.1))
具有对数因变量转换的附加婴儿模型
代码语言：javascript复制summary(abaln_ad_log)
均方根分数
代码语言：javascript复制kable(log_mseabalone_d_lo_f,"")
加性模型假设
代码语言：javascript复制model_assumptions多项式模型
代码语言：javascript复制lm(log(Rings) ~ Diameter Length Height poly(Whol_weiht,3) poly(Viscra_weight,3) polyShuckwegt,3) poly(Selweight,3) Infant, data=abaone_rain均方根分数
代码语言：javascript复制kable(log_rmse(abalnead_poly
多项式模型假设
代码语言：javascript复制model_assumption方差分析 F 检验
代码语言：javascript复制``````
anova交互模型
代码语言：javascript复制lm(log(Rings) ~ Height Diaeter poly(Whole_weight, 4)
poly(Viscera_weght, 4) poly(Shucke_weght,4) poly(Shel_weight, 4) Infant Diameter:poly(Shckedweight, 4) poly(Shckd_weight, 4):Infant, data=aboe_tain)均方根分数
代码语言：javascript复制kable(log_rms(abaoe_int_infan,"
交互模型假设
代码语言：javascript复制model_assumption方差分析 F 检验
代码语言：javascript复制anova
在这里，我们选择性能最佳的交互模型作为我们的候选模型之一。该模型具有较低的复杂性以及最低的 rmse。
我们在这里看到的是，与前两个完整模型相比，该模型在未见数据（测试数据）上的 rmse 最低。除此之外，这个模型也没有前两个复杂，这些特性使得这个模型成为我们选择的候选模型之一。
现在，由于我们有一些我们喜欢的候选模型，让我们将它们相互比较并寻找最佳模型。
候选模型
从上面的分析中，我们选择了多项式次数为 4 的交互模型，因为我们已经看到它具有最低的测试 rmse 之一，并且具有良好的模型假设。我们将把它作为第一选择的候选模型之一。该模型包含一些预测变量，而不是具有高度多重共线性的预测变量。这个模型是 abalone_model_int_bic。
另一个选定的模型是所有预测变量都具有相似较低 rmse 的交互模型。虽然这个模型将有更多的参数，因此是一个比早期模型更复杂的模型，但由于我们的目标是预测，我们不关心模型的复杂性。这个型号是 abalone_int_full
另一种候选模型将是具有 2 个因子水平的分类变量的交互模型。在开始或我们的分析中，我们看到该预测变量的分布在级别之间相似 male ， female 因此我们决定引入一个新变量，该变量基本上有 2 个级别，婴儿和非婴儿，其中非婴儿 = 雄性和雌性。该模型也是多项式次数为4的交互模型 。这个模型是 abalone_int_infant。注：（附录部分增加了使用婴儿预测器的加法模型和相关分析）。使用分类预测变量和具有3 值和 2 值的因子水平进行分析的整个过程 是相同的。
候选模型的 RMSE 比较
- 候选模型1（abalone_model_int_bic）
``````
kable(log_rmseabaone_modlint_bic
- 候选模型2（abalone_int__candidate_full）
kable(log_rmse(aalone_int_l,pte
kable(log_rmse(abalone_nt_nan- 候选模型 3 (abalone_int_infant)
- 从上面我们可以看到， 候选模型 3 (abalone_int_infant) 在所有三个模型中表现最好。
候选模型的模型假设
- 候选模型1（abalone_model_int_bic）
候选模型的模型假设
- 候选模型1（abalone_model_int_bic）
model_assumptions(abalone_model_int_bic,"Candidate model 1")
- 候选模型 2 (abalone_int_full)
``````
model_assumption
- 候选模型 3 (abalone_int_infant)
model_assumption- 基于这些图，我们可以看到我们所有的候选模型都很好地满足了恒定方差和正态性假设。
- 接下来，我们将比较所有这些候选模型的预测和置信区间，并观察数据中的任何影响点。让我们看看它。
异常观察样本：高杠杆、离群点和影响点
我们将查看每个候选模型的训练数据集中存在的异常观察结果。首先我们将分别寻找高杠杆点、离群点和影响点，然后通过公共数据点得到具有高杠杆、离群点和影响力的数据点。
候选模型1（abalone_model_int_bic）
代码语言：javascript复制## [1] 16- 我们看到有 16 个观察值是高杠杆、异常值和影响点。
候选模型2（abalone_int_full）
代码语言：javascript复制length(comon_mdel_2)## [1] 15- 在这个模型中，我们看到有 15 个观察值是高杠杆、异常值和有影响的点。
候选模型3（abalone_int_infant）
代码语言：javascript复制length(comon_odl_3)## [1] 19- 这个模型有 19 个这样的观察结果，它们是高杠杆、离群点和有影响力的点。
所有候选模型之间的异常观察
- 我们现在已经编制了候选模型中常见的观察结果列表。让我们来看看这些，看看是否有任何相似之处。
``````
intersect(commodl,intersect(commodel1,omomodl_2))
length(cmmondaa)## [1] 13我们立即注意到的第一件事是，我们最初在数据集中确定为不合逻辑或数据输入错误的观察结果很少。这些观察结果具有负权重差异，这是不可能的。这是有道理的，因为在记录此数据时可能会出现一些错误。
另一种模式是其中一些观察结果的整体权重与平均整体权重相比超过 2。此外，这些观察结果具有大量环，即这些是较老的鲍鱼。
由于存在其他符合异常观察条件的观察结果，因此总体而言，它们之间没有特定的模式。因此，我们不能从我们正在考虑进行预测的模型和分析中排除这些。
我们将使用平行坐标图来可视化结果部分中的任何明显模式。
让我们看看一旦我们删除这些数据集，我们的候选模型的表现如何。
删除异常观察后的RMSE 分数
接下来，我们将通过从训练数据集中删除异常观察并计算测试数据集的 RMSE 分数来执行快速实验。
代码语言：javascript复制obs_rmse <- get_log_rmse_obs(lmlog(_updated)
result <- ata.frae( "Mdel" = modelnames,
"RMSE Train" = c(
),
"RMSE Test" = c(
"RMSE Train all Data" = c(
"RMSE Test all data" = c(
get_lo_ms(aaloe_modl_intbic)$test, - 有趣的是，所有在没有异常观察的情况下训练的模型在测试数据上产生了更好的 rmse。在上表中我们可以看到不同RMSE值的比较。
- 由于这些不寻常的观察没有揭示任何特定的模式或任何不寻常的行为，我们无法从我们的数据集中删除这些，并将使用候选模型继续我们的分析。
预测和置信区间
代码语言：javascript复制#model 1
vif(abalone_mdl_ntbic)
#model 2
#model 3让我们从置信区间开始。我们所有候选模型的置信区间都在相似的范围内。即平均而言，在任何一点，所有模型的区间范围似乎都相同。请注意，3 个候选模型中有 2 个使用了所有预测变量，因此具有非常高的共线性。由于我们的模型很好地拟合了数据，因此平均减少了任何点的方差范围。这很好。
此外，与置信区间类似，模型之间的预测区间也在相同的范围内。模型拟合再次对这个较窄的区间范围产生影响。请注意，预测区间比我们预期的置信区间宽得多。由于所有模型的预测区间即使在高度共线性之后也处于相似的范围内，我们会说这些是好的。我们使用默认的 95% 公式来计算所有候选模型的置信区间和预测区间。
基于我们对所有这些候选模型的分析，我们认为候选模型 3 abalone_int_infant 将是预测鲍鱼年龄的最佳选择。与其他两个候选模型相比，该模型更小，复杂度更低，但 rmse 最低。此外，模型之间的预测和置信区间非常相似。由于这些原因，我们将选择模型 abalone_int_infant 作为我们分析的最佳模型。
预测鲍鱼的年龄
- 现在让我们使用我们选择的模型预测鲍鱼的年龄。我们将从测试数据中随机选择 5 个观察值并尝试预测年龄。请注意，该模型实际上预测了鲍鱼中的环数。我们需要在环数中加上 1.5 才能确定鲍鱼的真实年龄。所以本质上，我们将比较模型预测的环数与测试数据集中的实际环数。
``````
#从测试数据中随机抽取5个观察值
sample <- sample_n(test_data, 5)
exp(predict( ta=sample )我们可以从上面的例子中看到，我们从测试和模型中计算了环数和年龄，用于少数观察。我们看到我们的模型正确预测了鲍鱼的年龄，或者非常接近鲍鱼的实际年龄。这些是模型以前从未见过的测试数据中的 5 个样本观察值。我们可以多取一些并执行相同的过程，看看我们的模型对鲍鱼年龄的预测效果如何。
在这份分析报告中，我们应用了各种多元回归技术和方法，但这并不意味着这个问题不能使用其他更先进的技术来解决。我们尝试应用一些超出本报告范围的先进方法，只是为了评估这些技术是否可以进一步改进我们的预测。
使用先进技术
- 在本节中，我们超越了传统的常见回归技术，并尝试应用一些先进的技术，看看我们是否可以创建一个具有较低测试 RMSE 的模型。
Elastic Net 弹性网络是一种正则化回归方法，它线性组合了 lasso 和 ridge 方法的 L1 和 L2 惩罚。
- 我们将尝试使用弹性网络来进一步优化我们的预测。
- 与 lasso 和 ridge 一样，截距不会受到惩罚，并且 glment 在内部负责标准化。
- 我们将使用与候选模型 3 (abalone_int_infant) 中相同的预测变量。
train(
Rings ~ke_ight, 4) pol(Shell_weigh, 4) Infant Diameter:poly(Shucked_weight, 4) poly(Shuckedwight, 4:Infant,daa=abalone_trai
method = "glmnet",
trControl = cv_10,
kable(result, digit
- 我们看到训练和测试 rmse 结果都低于我们的候选模型。
随机森林回归
随机森林或随机决策森林是一种用于分类、回归和任务的集成学习方法，它通过在训练时构建大量决策树并输出类别（在分类的情况下）或平均预测来进行操作（在回归的情况下）单个树。随机决策森林纠正了 Doe 决策树过度拟合训练数据集的习惯。
我们在这里实现了随机森林回归模型进行预测，看看我们是否可以进一步改进。随机森林回归模型的 RMSE 结果比候选模型的 RMSE Score 更好。
代码语言：javascript复制 bm(Rig ~, daa aloeran stbton = "gausian", ntre =, inrtindeth 3, shrinkage = 0.01)
以下是随机森林方法的一些优缺点
它不对数据分布做出任何假设。
自动拟合高度非线性的交互。
通过接近度很好地处理缺失值。
即使对于大型数据集也能快速拟合。
已经观察到随机森林对某些具有嘈杂分类/回归任务的数据集过度拟合。
如果树不小就不容易解释
结果
本报告系统地处理了预测鲍鱼年龄的任务。探索了第一个数据集，发现了不寻常的观察结果，并对变量进行了适当的标准化。
我们从一个加法模型开始，该模型包含所有名为abalone_add的预测变量 。我们继续检查各种假设和测试，以进一步提高模型的准确性并降低 rmse。
然后，我们对因变量和预测变量进行了多次转换。当我们执行各种变量选择程序时，我们看到了 rmse 的改进。
经过我们的分析，我们找到了几个不错的候选模型，我们选择了模型 abalone_int_infant 作为我们的最佳模型。
以下是一些重要的关键点，说明了我们对这一结论的选择。
应用上述技术后，我们得出了 3 个可能的候选模型。
候选模型 1 是 abalone_model_int_bic ，这是对交互模型的各种选择过程的结果。我们看到这个模型 AIC 得分最低，rmse 排名第二。模型包含一些高共线预测变量，但当检查 rmse、预测和置信区间时，这是一个合适的选择。
候选模型 2 是 abalone_int_full。回想一下，详尽的搜索过程表明，当我们所有的预测变量，我们得到最好的模型。所以我们决定添加模型中交互项和多项式次数为 4 的所有预测变量。该模型具有最高的多重共线性，但从 rmse 和预测以及置信区间来看是非常合理的。
候选模型 3 是 abalone_int_infant。正如我们最初的数据分析表明分类变量 Sex 有 3 个因子水平。female， male 和 infant。水平 female 和 male 我们在分布方面极其相似，因此我们将这两个水平合二为一，使预测变量的总水平为 Sex 两个。由于这种级别的降低，与前两个模型相比，该模型可能是最简单的模型。我们对此进行了彻底的分析，最终模型具有相似的 rmse、预测和置信区间。这就是我们选择该模型作为候选模型之一的原因。
在以上 3 个候选模型中，候选模型 3 abalone_int_infant 稍微好一点，同时也不太复杂，因此它是我们进行此分析的首选模型。这有很多原因。让我们讨论一下。
以下是为此分析选择的模型的公式
代码语言：javascript复制 lm(log(Rings) ~ Height Diameter poly(Whole_weight, 4)
poly(Viscera_weight, 4) poly(Shucked_weight,4) poly(Shell_weight, 4) Infant Diameter:poly(Shucked_weight, 4) poly(Shucked_weight, 4):Infant- 这是所有 3 个候选模型的测试/训练 RMSE 的快速比较。
- 下图显示了整个分析报告中所选模型的 rmse（以红色标记）、候选模型和模型的其余部分以及它们的测试 rmse 比较。我们可以清楚地看到，我们选择的模型 abalone_int_infant 也在
int infant下面的图中重命名 ，其测试 rmse 最低。
- 这里一个明确的模式是，当模型更小、更简单时，测试 rmse 很高，但随着我们改进模型并开始变得更加灵活（因此模型长度更大），测试 rmse 下降。
- 以下是所选模型的正态性假设。我们可以清楚地看到，当我们用对数变换和交互项拟合模型时，常数变化和正态性假设变得越来越好。
这里要注意的一个关键点是，多重共线性的存在并没有真正影响预测（这是我们分析的目标），我们能够将其影响降低到一定程度。我们通过拟合与残差图验证了这一点，可以看出我们选择的模型具有以残差为中心的残差， 0 并且大致是恒定的变化。
如果我们使用平行坐标图绘制异常观测值，我们可以看到 Height 这些观测值与鲍鱼的平均高度相比更短。Whole_weight 和 Viscera_weight which as lower也是如此 。此外 weight.diff ，其中一些观察结果是否定的，这是我们在初始数据分析期间发现的。
下面的平行坐标图描绘了这些异常观察的特征，其中婴儿 = 1（橙色线）和非婴儿 [男/女] = 2（蓝线）。
现在，如果我们只是从我们的训练数据中抽取一个样本并使用平行坐标绘图，我们可以看到一般情况下 Height， Whole_weight 并且 Viscera_weight 均匀分布。
- 我们还看到从数据集中删除异常观察进一步降低了测试 rmse，但是我们没有从我们的数据集中删除那些，因为它们似乎是有效的假设而不是极端情况。
- 我们可以看到的均方根误差得分
Elastic Net和Random Forest比我们的候选模型较低，但差别不是很显著。
讨论
让我们从我们正在分析的数据集开始讨论。从最初的汇总统计数据中，我们看到数据中存在非常高的共线性。几乎所有的变量都有很高的 vif 值。
我们所做的另一个重要观察是数据质量。我们发现在某些情况下没有正确记录观察结果。例如，对于 2 个观察，高度为 0，这实际上是不可能的。但是当我们查看整个观察结果和其他变量时，我们没有发现任何异常模式。我们意识到这可能是数据输入错误，因此我们决定将它们保留在我们的数据集中进行分析。以下是这 2 个观察结果。
我们还看到，有些观察结果不符合逻辑。预测器 Whole_weight 被认为是预测的线性函数 Shucked_weight， Viscera_weight， Shell_weight 并在去壳过程中损失的血液/水的加成未知质量。换句话说，所有其他权重预测变量的总和应该等于或小于 Whole_weight。但是我们发现 153 这样的观察违反了这个假设。仔细观察后，我们发现观察结果之间没有相似之处或任何模式，因此我们得出结论，这可能只是数据记录错误，因此这些观察结果并未从数据集中删除。以下是这些 153 观察结果。
我们在探索性分析中发现的下一个有趣趋势是预测Sex 因子水平的分布 。因子水平female 和 male 所有其他预测因子 都非常相似。此外，整体数据也大致均匀地分布在female、 male 和 的 水平值之间 infant。以下是水平之间的分布及其相对于其他预测变量和响应变量的总体分布。
- 在上面的箱线图中，我们可以看到
female和 的分布形状male对于所有预测变量都非常相似。 - 现在让我们来看看它的频率。
同样，它非常相似。这些结果为我们指明了另一个方向，通过将female 和 组合male 成一个级别名称作为 ，我们在数据集中创建了一个新变量 non-infant。所以现在新的分类变量有两个级别 infant 和 non-infant。在我们的一个候选模型中，我们使用这个新的分类而不是原始的来查看其效果。最后，这是我们选择的模型。
初始数据分析期间的下一个重要发现是预测变量之间存在非常高的多重共线性。这是有道理的，因为我们之前发现 Whole_weight 预测在某种程度上依赖于其他重量预测。Length 鲍鱼的 Diameter. 以下是 来自加法模型的模型系数的 变异膨胀因子值。
代码语言：javascript复制## SexI SexM Length Diameter Height
## 1.982824 1.400405 42.874862 44.907429 3.114493
## Whole_weight Shucked_weight Viscera_weight Shell_weight
## 118.819537 29.579686 18.591058 22.182789初始模型违反了恒定方差和正态性假设。使用各种方差稳定技术和转换，我们能够提出一个很好的模型。这些选定的候选模型使用了交互项、预测变量的多项式变换和响应变量的对数变换。
有趣的是，即使存在高度多重共线性，在模型选择过程中，我们最终还是选择了具有所有 8 个参数的模型。即所有的选择程序都表明当使用模型中的所有预测变量时，最低的 AIC 出现。尽管在我们使用各种响应和预测变量转换进行模型选择之后，我们能够在保持相似的测试 rmse 的同时稍微减小模型的大小。
我们看到我们几乎没有对降低测试 rmse 有轻微影响的异常观察，但由于这些观察不代表显着变化或移动模型系数，我们没有将它们从数据集中删除。
在数据的背景下谈论这个分析，我们可以从模型中看到， Length 变量可能不会成为预测鲍鱼年龄的重要因素。我们可以看到我们没有在我们的模型中选择这个。这种情况是有意义的，因为 Length 或多或少地相关 Diameter （即长度越大直径越大），并且由于我们已经将直径带入我们的模型，所以我们不需要鲍鱼数据集中的长度。
预测时另一个非常有趣的事实是，不能单独使用整个体重来预测年龄。我们知道，鲍鱼的去壳过程会使鲍鱼的血液和水分变干，然后再称重，不同阶段（去皮、去壳、干燥）的重量组合也是决定鲍鱼年龄的重要因素。
我们从测试数据中抽取了 10 个样本观察值，并使用我们选择的模型 (abalone_int_infant) 预测了年龄，平均而言，它正确预测了 10 次中的第 3 次。从所选模型 (2.2021) 的测试 rmse 值考虑出响应变量环的范围（1 到 29），我们认为我们的模型非常适合预测鲍鱼的年龄。
