我想创建一个带有线段或箭头的图。
假设我有这个玩具示例
temp <- data.frame(posi=c(1,2,3,3,2,1,5), from=c("A", "B", "C", "D", "D", "B", "A"),
to=c( "C", "D", "D", "C", "A", "A", "B"))
posi from to
1 A C
2 B D
3 C D
3 D C
2 D A
1 B A
5 A B
我想绘制从“from”定义的点及其起始位置到“to”和position+1定义的点的线段或箭头。
我必须用 ggplot 创建绘图
ggplot() + geom_point(data=temp, aes(x=posi, y=from, size=10),show.legend = FALSE) +
geom_point(data=temp, aes(x=posi+1, y=to, size=10),show.legend = FALSE) +
geom_segment(data=temp, aes(x=posi, y=from, xend=posi+1, yend=to),
arrow=arrow(type="closed", angle=10), size=1.5, color="blue") + theme_bw()
但我想用莱迪思来做，因为我的真实数据集要大得多（数百万行）并且莱迪思工作得更快。 我知道我可以减少删除重复项的行数，但那是另一个故事了。
我怎样才能用格子做到这一点？
我一直在研究，我认为我需要使用 panel.segments 或 lsegments，但事情似乎要复杂得多，而且很难找到例子。
xyplot(from ~ posi , type="p", col="black", data=temp, pch=16,
panel = function(x, y, ...){ panel.segments(x,y,) })
我不知道要在函数内部或面板内部写入哪些参数。
您可以使用以下代码：
xyplot(from ~ posi , type="p", col="black", data=temp, pch=16, xlim = c(0,7),
panel = function(...){
panel.dotplot(x = (temp$posi+1), y = temp$to, col="black", cex=1.4)
panel.dotplot(x = temp$posi, y = temp$from, col ="black", cex=1.4)
panel.arrows(x0 = temp$posi, y0 = temp$from, x1 = temp$posi+1, y1 = temp$to, lwd=3, col="blue", )
}
)
产生下图：
请告诉我这是否是您想要的。
我发布了一个关于 @skan 在评论中识别和描述的问题的问题：当“极端”级别（如“D”）不存在于temp$from那么“D”将不会成为图表的一部分，即使稍后需要“D”temp$to。可以找到@Konn的答案的问题here.
据我现在的理解，我们需要一个factor那是ordered，我们需要添加指定的代码drop.unused.levels = FALSE在通话中xyplot。在示例中，我们显示了“from”中包含“extremes”的完整集合，以及作为不存在极端“D”的子集：完整代码为：
l <- c("A", "B", "C", "D")
temp <- data.frame(posi = c(1, 2, 3, 3, 2),
from= factor(c("A", "B", "C", "D", "D"), levels = l, ordered = TRUE),
to = factor(c("C", "D", "D", "C", "A"), levels = l, ordered = TRUE)
)
xyplot(from ~ posi , type="p", col="black", data=temp, pch=16, xlim = c(0,7),
drop.unused.levels = FALSE, ## the added code
panel = function(...){
panel.dotplot(x = temp$posi, y = temp$from, col ="green", cex=1.6)
panel.dotplot(x = (temp$posi+1), y = temp$to, col="black", cex=1.)
panel.arrows(x0 = temp$posi, y0 = temp$from, x1 = temp$posi+1, y1 = temp$to, lwd=2, col="blue" )
})
temp <- temp[1:3, ]
xyplot(from ~ posi , type="p", col="black", data=temp, pch=16, xlim = c(0,7),
drop.unused.levels = FALSE, ## the added code
panel = function(...){
panel.dotplot(x = temp$posi, y = temp$from, col ="green", cex=1.6)
panel.dotplot(x = (temp$posi+1), y = temp$to, col="black", cex=1.)
panel.arrows(x0 = temp$posi, y0 = temp$from, x1 = temp$posi+1, y1 = temp$to, lwd=2, col="blue" )
})
产生以下图片：
我想我们已经解决了这个问题。