2016-08-17 104 views
0

我不知道我是否制定了正确的标题/问题。也许我的一个问题是我的词汇中缺少术语。抱歉。但让我们尝试:如何在R中的3D图中绘制多个2d图?

我有数据(sleep在这个例子中)我会描述为三维。也许真正的统计学家不会那样做?

我想我想绘制多个二维情节到一个三维的情节。我想并排绘制它们。如果我错了,请纠正我。

enter image description here

我在这里的问题是,只有一条线。 有两组。我想要每组一行。与type='h'同样的数据提供更好的描述,我认为:

enter image description here

你能想象这两条线在这里?我在这个概念中缺少什么? 我们可以使用另一个绘图库进行打印/发布。目前这对我来说无关紧要。也许我总是在错误的地方?

这是代码:

require('mise') 
require('scatterplot3d') 
mise() # clear the workspace 

# example data 
print(sleep) 

scatterplot3d(x=sleep$ID, 
       x.ticklabs=levels(sleep$ID), 
       y=sleep$group, 
       y.ticklabs=levels(sleep$group), 
       lab = c(length(unique(sleep$ID)), 1), 
       z=sleep$extra, 
       type='o') 

而且数据

extra group ID 
1 0.7  1 1 
2 -1.6  1 2 
3 -0.2  1 3 
4 -1.2  1 4 
5 -0.1  1 5 
6 3.4  1 6 
7 3.7  1 7 
8 0.8  1 8 
9 0.0  1 9 
10 2.0  1 10 
11 1.9  2 1 
12 0.8  2 2 
13 1.1  2 3 
14 0.1  2 4 
15 -0.1  2 5 
16 4.4  2 6 
17 5.5  2 7 
18 1.6  2 8 
19 4.6  2 9 
20 3.4  2 10 

回答

1

你可以手工添加的行两步:

# Store the plot in rr 
rr <- scatterplot3d(x=as.numeric(sleep$ID), 
        x.ticklabs=levels(sleep$ID), 
        y=sleep$group, 
        y.ticklabs=levels(sleep$group), 
        z=sleep$extra) 
# find all that belong to group one 
idx = sleep$group == 1 
# add the first line 
rr$points3d(x = sleep$ID[idx], y = rep(1, each = sum(idx)), z = sleep$extra[idx], type = 'l', col = 'red') 
# add the second line 
rr$points3d(x = sleep$ID[!idx], y = rep(2, each = sum(!idx)), z = sleep$extra[!idx], type = 'l', col = 'blue') 

所以,而不是加线的丝带事情变了一点。特别是,这些色带是用polygon函数绘制的。但是,该函数只能处理二维坐标,所以我们需要使用函数rr$xyz.convert将我们的三维坐标转换为二维坐标。

rr <- scatterplot3d(x=sleep$ID, 
        x.ticklabs=levels(sleep$ID), 
        y=sleep$group, 
        y.ticklabs=levels(sleep$group), 
        z=sleep$extra) 
idx = sleep$group == 1 
# draw first group 
mat = matrix(c(rep(sleep$ID[idx], 2), 
       rep(c(1, 1.05), each = sum(idx)), # 1.05 determines width 
       rep(sleep$extra[idx], 2)), ncol = 3) 
ll = rr$xyz.convert(mat) 
polygon(x = ll$x[c(1:10, 20:11)], 
     y = ll$y[c(1:10, 20:11)], col = 'red') 
# draw second group 
mat = matrix(c(rep(sleep$ID[!idx], 2), 
       rep(c(2, 1.95), each = sum(!idx)), # 1.95 determines width 
       rep(sleep$extra[!idx], 2)), ncol = 3) 
ll = rr$xyz.convert(mat) 
polygon(x = ll$x[c(1:10, 20:11)], 
     y = ll$y[c(1:10, 20:11)], col = 'blue') 
+0

这并不完全适合我的主要问题,但我如何将这两条线画成丝带?海事组织这里的解决方案看起来不太好,并没有真正帮助理解数据。 – buhtz

+0

你能举一个你想要情节怎么样的例子图片吗?也许像这样但只有两行? http://i.stack.imgur.com/kBINP.png – Vandenman

+0

是的,这看起来不错。 – buhtz