R调用daewr包LGB函数报错‘找不到该函数’的解决求助

问题描述 #

运行无重复部分析因设计的显著效应分析代码时触发报错：

Error in LGB(coef(modc)[-1], rpt = FALSE) : could not find function "LGB"

已安装并加载daewr包，使用RStudio 2023.03.0 Build 386及最新版R，涉及代码如下：

Model.1 <- lm (Yield ~ (.)^2, data = Soup)
modc <- lm(y~(.)^2, data=soup)
library(daewr)
LGB(coef(modc)[-1], rpt = FALSE)

报错原因及解决办法 #

1. 核查daewr包版本与函数可用性 #

• 新版daewr包可能移除或更名了LGB函数，先确认当前包版本：

  packageVersion("daewr")
  

• 若版本较新，可回退至确认包含LGB函数的旧版本（如0.1.8）：

  install.packages("https://cran.r-project.org/src/contrib/Archive/daewr/daewr_0.1.8.tar.gz", repos = NULL, type = "source")
  

• 重装后重新加载包并运行代码：

  library(daewr)
  LGB(coef(modc)[-1], rpt = FALSE)
  

2. 排查代码细节问题 #

• 确认modc模型拟合有效，检查系数向量是否正常：

  # 查看完整模型系数
  print(coef(modc))
  # 确认去掉截距后的向量
  print(coef(modc)[-1])
  

• 注意数据集名称大小写：Soup与soup是不同对象，确保soup数据集已正确加载且包含变量y。

替代工具与包 #

若daewr的LGB函数无法使用，可选择以下方案：

• FrF2包：专注2水平因子设计分析，支持无重复部分析因设计的效应筛选，通过半正态图识别显著效应：

  library(FrF2)
  # 先创建部分析因设计对象（填入实际因子数）
  design <- FrF2(nruns = nrow(soup), nfactors = 4)
  model <- lm(y ~ (.)^2, data = soup)
  DanielPlot(model, alpha = 0.05)
  

• effects包：提取并可视化模型效应，辅助判断显著性：

  library(effects)
  plot(allEffects(model))
  

• 手动绘制半正态图：基础R函数实现效应显著性判断：

  effects <- abs(coef(modc)[-1])
  qqnorm(effects, main = "半正态图")
  qqline(effects, col = "red")
  

内容的提问来源于stack exchange，提问作者user3569147