在R的网格包中添加3D abline到云图

时间:2011-10-08 20:32:58

标签: r plot lattice scatter3d

我想在R的网格包中为一个云散点图添加一个3D基线。这是我数据的一个子集(3个变量都在0,1之间):

dat <- structure(c(0.413, 0.879, 0.016, 0.631, 0.669, 0.048, 1, 0.004, 0.523, 0.001, 
0.271, 0.306, 0.014, 0.008, 0.001, 0.023, 0.670, 0.027, 0.291, 0.709, 
0.002, 0.003, 0.611, 0.024, 0.580, 0.755, 1, 0.003, 0.038, 0.143, 0.214, 
0.161, 0.008, 0.027, 0.109, 0.026, 0.229, 0.006, 0.377, 0.191, 0.724, 
0.119, 0.203, 0.002, 0.309, 0.011, 0.141, 0.009, 0.340, 0.152, 0.545, 
0.001, 0.217, 0.132, 0.839, 0.052, 0.745, 0.001, 1, 0.273), .Dim = c(20L, 3L))

这是云图:

# cloud plot
trellis.par.set("axis.line", list(col="transparent")) 
cloud(dat[, 1] ~ dat[, 2] + dat[, 3], pch=16, col="darkorange", groups=NULL, cex=0.8, 
    screen=list(z = 30, x = -70, y = 0),
    scales=list(arrows=FALSE, cex=0.6, col="black", font=3, tck=0.6, distance=1) ) 

我想在0,0,0和1,1,1之间添加一条虚线灰色线(即对角线穿过图)。我知道我可以使用“type =”l“,panel.3d.cloud = panel.3dscatter”将点更改为行,但我看不到添加额外点/行的方法情节用这个。

这是我想用scatterplot3d实现的一个例子:

# scatterplot3d
s3d <- scatterplot3d(dat, type="p", color="darkorange", angle=55, scale.y=0.7,
    pch=16, col.axis="blue", col.grid="lightblue") 

# add line 
s3d$points3d(c(0,1), c(0,1), c(0,1), col="grey", type="l", lty=2)

我想用云图来控制我查看绘图的角度(scatterplot3d不允许我让绘图的0,0,0角面对)。感谢您的任何建议。

1 个答案:

答案 0 :(得分:2)

不雅,可能很脆弱,但这似乎有用......

cloud(dat[, 1] ~ dat[, 2] + dat[, 3], pch=16, col="darkorange", 
        groups=NULL, cex=0.8, 
    screen=list(z = 30, x = -70, y = 0),
    scales=list(arrows=FALSE, cex=0.6, col="black", font=3, 
                tck=0.6, distance=1) ,
      panel=function(...) {
        L <- list(...)
        L$x <- L$y <- L$z <- c(0,1)
        L$type <- "l"
        L$col <- "gray"
        L$lty <- 2
        do.call(panel.cloud,L)
        p <- panel.cloud(...)
      })

要记住的一点是,这不会删除隐藏的点/线,因此该线将位于所有点的前面或后面的所有点;在这个(编辑过的)版本中,do.call(panel.cloud,L)是第一个,所以这些点会遮挡这条线而不是相反。如果你想要删除隐藏行,那么我相信rgl是你唯一的选择......非常强大,但不是很漂亮,而且界面更原始。