在阅读了重新设定薄板样条和“字段”包(here,here和here)的帖子后,我有一些跟进问题:
考虑此代码使用光圈数据集
创建薄板样条data(iris)
library(fields)
library(plot3D)
x <- sep.l <- iris$Sepal.Length
y <- pet.l <- iris$Petal.Length
z <- sep.w <- iris$Sepal.Width
test.spline <- Tps(data.frame(x,y),z)
x2 <- seq(min(x), max(x), len=250)
y2 <- seq(min(y), max(y), len=250)
data <- expand.grid(x=x2, y=y2)
fitted <- matrix(predict(test.spline, data), 250,250)
persp3D(z=fitted, facets = FALSE, contour = TRUE, theta = 70, phi=30, ticktype="detailed", xlab="", ylab="", zlab="")
上面的代码创建了一个薄板样条,然后将模型拟合到网格上(我希望)。
谢谢大家的帮助。