我在数据框中有以下数据:
**x** in (0,1)
**y** in [0,1]
**z** in [0,1]
例如:
X,Y,Z
0.1, 0.2, 0.56
0.1, 0.3, 0.57
...
我想把它们画在这类图表上:
我试了R,但我只能得到一个not-so-fancy 3d scatterplot。
我也读过关于lattice 3d wireframe的文章,但我无法理解。
我该怎么做才能在R中得到一个类似Matlab的线框呢?
涉及哪些数据转换?
这是文档中的示例代码:
x <- seq(-pi, pi, len = 20)
y <- seq(-pi, pi, len = 20)
g <- expand.grid(x = x, y = y)
g$z <- sin(sqrt(g$x^2 + g$y^2))
wireframe(z ~ x * y, g, drape = TRUE,
aspect = c(3,1), colorkey = TRUE)
我觉得不是特别清楚。
编辑:
persp3d
功能运行良好,我能够生成一个单一颜色的三维绘图如何设置相对于z值的色阶?谢谢你的提示,
多酮
最佳答案
使用outer
创建z值,然后使用persp
绘制:
z <- outer(x,y, function(x,y) sin(sqrt(x^2+y^2)))
persp(x,y,z)
有颜色和设置视角的选项,请参见
?persp
关于Matlab风格的着色,请参见第四个示例。对于交互式绘图,请考虑在
persp3d
包中使用rgl
:require(rgl)
persp3d(x,y,z,col="blue")
编辑
要添加颜色,与
persp
中的方法略有不同,因为颜色与顶点而不是镶嵌面的中心有关,但这样做更容易。jet.colors <- colorRampPalette( c("blue", "green") )
pal <- jet.colors(100)
col.ind <- cut(z,100) # colour indices of each point
persp3d(x,y,z,col=pal[col.ind])
帮助文件建议添加参数
smooth=FALSE
,但这取决于个人偏好。