尝试使用 R 中的 model 预测列中的值

Question

我会预先说明 state 我在 R 方面的经验很少，并且非常感谢人们可以提供的任何帮助。

我有一个值的数据框（我认为）：

在此处输入图像描述

我正在尝试使用我从其他数据创建的 model 以及“ThatwgFAC”和“ThalwgSlop”的值来预测“PredBFwdth”列的每个值。

model的名称为“型号”； 它是一个广义的添加剂 model，所以我会使用 predict.gam。 我认为代码看起来像：

PredBFwdth = predict.gam(Model, ThalwgFAC, ThalwgSlop)

这显然行不通。 我需要使用来自 ThalwgFAC 和 ThalwgSlop 列以及 model 的相应值对数据中的每一行进行预测。 有人可以帮忙吗？ 我四处寻找答案，但其他人问题的代码对我来说是难以理解的。 :(

Answer 1

我强烈建议提供足够的数据和代码来创建可重现的示例。 这意味着您的数据可能是什么样的示例以及您正在尝试使用哪些代码以促进帮助。

首先，我将从 data.frame 中的一些假数据开始

dat <- data.frame(

ThalwgFAC = rnorm(100,3,5),
ThalwgSlop = rnorm(100,3,2),
ThalwgEv = rnorm(100,2,2)
)

接下来，您安装您的 model

library(mgcv)

fit <- gam(ThalwgEv ~ s(ThalwgFAC) + s(ThalwgSlop), data = dat)

现在关于预测新数据的问题，您可以将新数据（作为数据框）提供给predict function 并指定您要使用的 model （上面的fit ZA8CFDE6331BD59EB2AC96F8911C46）。 这假定数据帧包含用于 model 的相同变量名称（例如 ThalwgFAC 和 ThalwgSlop）。

dat2 <- data.frame(
  
  ThalwgFAC = rnorm(100,3,5),
  ThalwgSlop = rnorm(100,3,2),
  ThalwgEv = rnorm(100,2,2)
)

pred <- predict(fit, newdata = dat2)

predict function 将选取fit object 的 class，因此无需指定使用的方法（即 predict.gam）。