r - 绘制薄板样条模型时出现的错误

Question

library(mgcv)
#Input dataset
vo2max<-read.csv('C:/Users/Dell/Desktop/Paper(2009)/vo2max.csv')
model1 <- gam(VO2max ~ s(Load, Speed, m=4,  bs=c("tp", data = vo2max) 
x <- range(vo2max$Load)
x <- seq(x[1], x[2], length.out=14)
y <- range(vo2max$Speed)
y <- seq(y[1], y[2], length.out=14)
z <- outer(x,y,
           function(Load,Speed)
             predict(model1, data.frame(Load,Speed)))
p <- persp(x,y,z, theta=30, phi=30,
           col="yellow",expand = 0.5,shade = 0.2, ticktype = "detailed",
           xlab="Load", ylab="Speed", zlab="VO2max")
obs<- trans3d(vo2max$Load, vo2max$Speed,vo2max$VO2max,p)
pred<- trans3d(vo2max$Load, vo2max$Speed,fitted(model1),p)`enter code here`
points(obs, col="red",pch=16)
segments(obs$x, obs$y, pred$x, pred$y)
predict(model1)

数据集包含 14 个值。我需要绘制一个薄板样条曲线并预测模型，但我在平滑构造中遇到了错误。我是否必须使用任何其他库包。帮助我提供适当的代码。

标签： rmgcv