Добавьте 3D аблайн к облачному графику в решетчатом пакете R - PullRequest
Новая
       38

Добавьте 3D аблайн к облачному графику в решетчатом пакете R

3 голосов
/ 09 октября 2011

Я хочу добавить трехмерную аблайн к графику рассеяния облака в решетчатом пакете R.Вот подмножество моих данных (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 Ответ

2 голосов
/ 09 октября 2011

Не элегантный и, возможно, хрупкий, но это, кажется, работает ... 

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 - ваш единственный вариант ... очень мощный, но не такой красивый и с гораздо более примитивным интерфейсом.

...