r - locfit を plot3d するにはどうすればよいですか?

Question

次のようにプロットすると点群を形成する変数 x、y、z があります。

library(plot3D)
plot3d(x,y,z)

私はそれらにロックフィットをしました

myfit = locfit(y~lp(x,z),maxk=200)

これにより、空間の最も密度の高い領域を通過する曲線が得られることを理解しています。

この曲線を plot3d / RGL でプロットするにはどうすればよいですか?

Answer 1

を使用しsurface3dます。x と y は余白のベクトルで、z は行列です。

require(locfit)
fit <- locfit(NOx~lp(E,C,nn=0.5,scale=0), data=ethanol)
plot(locfit)  # there is an ordinary contour plot method for locfit objects.

require(rgl)
open3d()
surface3d( x=seq(0.5, 1.3, by=0.1), y=seq(7.5,18,by=.5) ,
           z= matrix(  predict(fit, newdata= 
                           expand.grid(E=seq(0.5, 1.3, by=0.1), 
                                       C=seq(7.5,18,by=.5) ) ) ),
                      ,nrow= length(seq(0.5, 1.3, by=0.1)) ,
                       ncol= length(seq(7.5,18,by=.5) ) ,
           xlim=c(.5, 1.3) )
 # grab and spin

実際には、等高線図の方が有益だと思いますが、3D プロットも役立ちます。