r - scatterplot3d を使用して薄板スプラインをプロットする

Question

スプラインは私にとってまだかなり新しいものです。

統計学習入門 ( http://www-bcf.usc.edu/~gareth /ISL/ISLR%20Sixth%20Printing.pdf )。私は scatterplot3d で作業しています。データを簡単に再現できるようにするために、実際のデータの代わりに「ツリー」データセットを使用します。

初期プロットの設定は簡単です:

data(trees)
attach(trees)
s3d <- scatterplot3d(Girth, Height, Volume,
                 type = "n", grid = FALSE, angle = 70,
                 zlab = 'volume',
                 xlab = 'girth', 
                 ylab = 'height',
                 main = "TREES") # blank 3d plot

フィールド ライブラリの Tps 関数を使用して、スプラインを作成します。

my.spline <- Tps(cbind(Girth, Height), Volume)

そして、スプラインを視覚的に表現し始めることができます:

for(i in nrow(my.spline$x):1) # for every girth . . .
s3d$points3d(my.spline$x[,1], rep(my.spline$x[i,2], times=nrow(my.spline$x)), # repeat every height . . . 
              my.spline$y, type='l') # and match these values to a predicted volume

しかし、高さのアクセスに沿ってハッチング線を交差させてスプラインを完成させようとすると、結果に問題が生じます。

for(i in nrow(my.spline$x):1) # for every height . . .
s3d$points3d(rep(my.spline$x[i,1], times=nrow(my.spline$x)), my.spline$x[,2],  # repeat every girth . . . 
           my.spline$y, type='l') # and match these values to a predicted volume 

そして、結果のプロットを見れば見るほど、my.spline の正しいデータを使用していることさえ確信が持てなくなります。

このプロジェクトは他のビジュアライゼーションに scatterplot3d を使用することに注意してください。そのため、既存のチームの選択の結果として、このパッケージに夢中になっています。どんな助けでも大歓迎です。