【问题标题】:in R, plot sucessive/sequence events在 R 中,绘制连续/顺序事件
【发布时间】:2015-05-18 21:51:16
【问题描述】:

已编辑:遵循建议

我正在尝试根据这篇文章创建一个简单的连续事件频率图:ggplot sequence patterns

我想为其他主题(不同的序列)堆积相同的情节 主题"aa", "bb", "cc"

示例数据集如下:

subject <- c("aa","aa","aa","aa","aa","aa","aa", "bb","bb","bb","bb","bb","bb","bb","cc","cc","cc","cc","cc","cc","cc")
event <- c("P","C","D","C","P","E","D","C","P","E","D","P","C","E","P","E","D","P","C","x","x")
freq <- c(3,2,1,4,2,1,0,3,4,1,3,3,1,2,1,3,2,1,4,0,0))

dfx <- data.frame(subject, event, freq)

我得到的结果是:

使用此代码,基于原帖:

library(ggplot2)

dfx$type <- factor(dfx$type)
dfx$ymin <- c(0,1,2)
dfx$ymax <- c(1,2,3)
dfx$xmax <- cumsum(dfx$count)
dfx$xmin <- c(0, head(dfx$xmax, n=-1))

plot_x <- ggplot(dfx, 
            aes(xmin=xmin, xmax=xmax, ymin=ymin, ymax=ymax,fill=type))   +geom_rect(colour="grey40", size=0.5)

png("plot_x.png", height=200, width=800)
print(plot_x)
dev.off()

我已经用我想要绘制的内容创建了这个图像。 (它是用excel手工制作的)。在这种情况下,我们有 3 个主题、4 个事件(C、P、D、E)+1 个虚拟事件(X),这是创建数据框所必需的。如您所见,每个主题的事件总数不一定相等。

【问题讨论】:

  • 您还应该分享“结果”部分的代码。
  • 如果与您提供的示例数据相关,您的目标输出示例会更容易理解。要么给出你所拥有的示例背后的数据,要么根据你所拥有的数据制作一个与数据相匹配的示例
  • 感谢您的推荐!我已经相应地编辑了我的帖子。

标签: r ggplot2 rcharts sequence-diagram


【解决方案1】:

尝试使用facet_gridfacet_wrap

library(ggplot2)
p1  <-  
    ggplot(dfx, aes(xmin = xmin, xmax = xmax, ymin = ymin, ymax = ymax, fill = event))+
    theme_bw() + 
    geom_rect(colour = "black", size = 0.5) + 
    facet_grid(subject ~ .) + 
    theme(axis.ticks.y=element_blank(), axis.text.y=element_blank())

p1

数据:

dfx <- data.frame(
       subject = c("aa","aa","aa","aa","aa","aa","aa", "bb","bb","bb","bb","bb","bb","bb","cc","cc","cc","cc","cc","cc","cc"),
       event = c("P","C","D","C","P","E","D","C","P","E","D","P","C","E","P","E","D","P","C","x","x"),
       freq = c(3,2,1,4,2,1,0,3,4,1,3,3,1,2,1,3,2,1,4,0,0),
       ymin = 0,
       ymax = 1)

TEMP        <-  tapply(dfx$freq,dfx$subject,cumsum)
dfx$xmax    <-  unlist(TEMP)
dfx$xmin    <-  c(0, head(TEMP$aa, -1), 0, head(TEMP$bb, -1), 0, head(TEMP$cc, -1))

【讨论】:

  • 谢谢!我已经对其进行了测试,这是一个非常简单且很好的解决方案。非常感谢!!!
猜你喜欢
  • 2020-04-24
  • 1970-01-01
  • 1970-01-01
  • 1970-01-01
  • 1970-01-01
  • 1970-01-01
  • 2016-09-26
  • 1970-01-01
  • 1970-01-01
相关资源
最近更新 更多