#导入数据文件
data <- read.csv("插补1.0.csv")
#导入包
library(mgcv)
# 定义因变量和自变量
dependent_variables <- c("胎儿预测体重", "新生儿一分钟评分", "新生儿性别")
independent_variables <- c("TEP","TNBP","TCIPP","TPHP","TMPP","TEHP","EHDPP")
# 加入其他变量
other_variables <- c("年龄", "孕周", "胎次","吸烟与否")
# 合并自变量和其他变量
all_variables <- data.frame(all_variables[, independent_variables], all_variables[, other_variables])
# 构建三次样条的GAM模型
gam_model <- gam(cbind(胎儿预测体重, 新生儿一分钟评分, 新生儿性别) ~ s(TEP) + s(TNBP) + s(TCIPP) + s(TPHP) + s(TMPP) + s(TEHP) + s(EHDPP) + 年龄 + 孕周 + 胎次 + 吸烟与否, data = all_variables)
# 查看模型摘要
summary(gam_model)
一运行就找不到对象,这是为什么?没有命名的原因吗?要怎么解决啊? # 构建三次样条的GAM模型 > gam_model <- gam(cbind(胎儿预测体重, 新生儿一分钟评分, 新生儿性别) ~ s(TEP) + s(TNBP) + s(TCIPP) + s(TPHP) + s(TMPP) + s(TEHP) + s(EHDPP) + 年龄 + 孕周 + 胎次 + 吸烟与否, data = all_variables) Error in eval(predvars, data, env) : 找不到对象'胎儿预测体重'