我正在尝试创建一个圆形图并且卡在某一点上:
dat1 <- data.frame (xvar = 1:10, y = 6, ymin = 4, ymax = 4.5)
使用这些数据,我可以在ggplot2中生成一个圆形带状图
require(ggplot2)
ggplot(dat1, aes(x=xvar, y=y)) + geom_ribbon(aes(ymin=ymin, ymax=ymax),
col = "blue", fill = "blue2") + ylim (c(0,6)) + coord_polar()
然而我想要更多。
我想使用以下数据用不同的颜色和标签填充功能区的片段。
filld <- data.frame (start = c(1, 4, 6, 7.5, 8, 9), end = c(4, 6, 7.5, 8, 9, 10),
label = c("A", "B", "C", "A", "C", "D"))
filld
## start end label
## 1 1.0 4.0 A
## 2 4.0 6.0 B
## 3 6.0 7.5 C
## 4 7.5 8.0 A
## 5 8.0 9.0 C
## 6 9.0 10.0 D
色带将通过标签变量填充不同的颜色。例如,段A将从1开始并以4结束。然后段B将以6开始和结束并填充不同的颜色。具有相同标签的段(例如A和C)将按行连接。
结果图将如下所示:
答案 0 :(得分:19)
以下是一个例子:
filld$p <- rowMeans(subset(filld, select = c(start, end)))
ggplot(filld, aes(xmin = start, xmax = end, ymin = 4, ymax = 5, fill = label)) +
geom_rect() +
geom_segment(data = subset(filld, label %in% label[duplicated(label)]),
aes(x = p, y = 0, xend = p, yend = 4, colour = label),
size = 2, show_guide = FALSE) +
geom_text(aes(x = p, y = 4.5, label = label), colour = "white", size = 10) +
coord_polar() +
scale_y_continuous(limits = c(0, 5))
更新
我不建议这样的事情:
filld <- data.frame (start = c(1, 4, 6, 7.5, 8, 9), end = c(4, 6, 7.5, 8, 9, 10),
label = c("A", "B", "C", "A", "C", "D"))
filld$p <- rowMeans(subset(filld, select = c(start, end)))
filld <- merge(filld, ddply(filld, .(label), summarize, p2 = mean(p)))
lnd <- subset(filld, label %in% label[duplicated(label)])
lnd <- ddply(lnd, .(label), function(x) {
x <- seq(x$p[1], x$p[2], length = 100)
y <- 4.5 + ((x - mean(x))^2 - (x[1]-mean(x))^2) / (x[1]-mean(x))^2 * 3 + sin(x*3*pi) * 0.1
data.frame(x, y)
})
p <- ggplot(filld, aes(xmin = start, xmax = end, ymin = 4, ymax = 5, colour = label, fill = label)) +
geom_line(aes(x, y, xmin = NULL, ymin = NULL, xmax = NULL, ymax = NULL), data = lnd, size = 2) +
geom_rect() +
geom_text(aes(x = p, y = 4.5, label = label), colour = "white", size = 10) +
coord_polar() +
scale_y_continuous(limits = c(0, 5))
p
也许,你想要的东西超出了ggplot2的范围。
答案 1 :(得分:7)
看看ggbio。这个包将ggplot2和图形语法扩展到序列数据(生物信息学),但它们似乎解决了你的问题(例如见here)。看一下包的source code可能会引导您找到更通用的解决方案。
答案 2 :(得分:6)
添加ymin等以填充d
fd <- data.frame(filld, ymin = 4, y =6, ymax = 4.5)
然后使用geom_rect
列label
作为fill
和colour
美学
ggplot(fd, aes(x=start,xmin = start, xmax = end, y=y)) +
geom_rect(aes(ymin=ymin, ymax=ymax, fill = label )) +
ylim (c(0,6)) +
coord_polar()
添加以下行:
## calculate the midpoint for each segment
fd$xmid <- apply(fd[,1:2],1,mean)
## get the replicate id for the labels (what occurence is this)
library(plyr)
library(reshape2)
fd1 <- ddply(fd, .(label), mutate, id = 1:length(xmid))
## reshape to wide, subsetting only with more than one rep
.lines <- na.omit(dcast(fd1, label~id, value.var = 'xmid'))
## add a mid point between the mid points
.lines$mid <- apply(.lines[,2:3],1,mean)
## reshape to long, and add some y values
ld <- data.frame(arrange(melt(.lines,value.name = 'x'), label, x), y = rep(c(4,3,4),2))
ggplot(fd) +
geom_rect(aes(x=start,xmin = start, xmax = end, y=y,ymin=ymin, ymax=ymax, fill = label )) +
ylim (c(0,6)) +
coord_polar() + geom_path(data = ld, aes(x=x,y=y,colour = label))
线条很难看,但是那里!