不完全な楕円に似た 2D 座標雲を当てはめるつもりでした。R の「splines」パッケージのbs()関数を使用してフィッティング曲線 (多項式スプラインの B スプライン基底行列を生成するために使用される) を取得しようとしたとき、次のようになりました。
しかし、私が意図したのはこのようなものでした。
パラメータを bs() にフィードしてプロットするには、これに従いました
mod <- lm(y ~ bs(x, knots = seq(0.1, 0.9, by = 0.1)))
pdat <- data.frame(x = seq(min(x), max(x), length = 100))
## predict for new `x`
pdat <- transform(pdat, yhat = predict(mod, newdata = pdat))
# plotting
plot(y ~ x)
lines(yhat ~ x, data = pdat, lwd = 2, col = "red")
私が意図したことを達成する方法はありますか？
ありがとう！