我如何绘制一个locfit?

时间:2014-06-25 20:32:57

标签: r

我有变量x,y,z,当绘制时形成点云:

plot3d(x,y,z)

我在他们身上做了 locfit

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

据我所知,这会给我一条穿过最密集的太空区域的曲线。

如何在plot3d / RGL中绘制此曲线?

1 个答案:

答案 0 :(得分: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图也很有用。