用ggplot绘制纵向数据时的错误分组

时间:2017-07-25 14:12:40

标签: r ggplot2 longitudinal

我试图用ggplot绘制我的纵向数据集。   以下是用于说明我的问题的示例数据。

library(MASS)
library(ggplot2)
library(dplyr)

# create data    
dat.tx.a <- mvrnorm(n=250, mu=c(30, 20, 28), 
                        Sigma=matrix(c(25.0, 17.5, 12.3, 
                                       17.5, 25.0, 17.5, 
                                       12.3, 17.5, 25.0), nrow=3, byrow=TRUE))
    dat.tx.b <- mvrnorm(n=250, mu=c(30, 20, 22), 
                        Sigma=matrix(c(25.0, 17.5, 12.3, 
                                       17.5, 25.0, 17.5, 
                                       12.3, 17.5, 25.0), nrow=3, byrow=TRUE))
    dat <- data.frame(rbind(dat.tx.a, dat.tx.b))
    names(dat) = c("measure.1", "measure.2", "measure.3")
    dat <- data.frame(subject.id=factor(1:500), tx=rep(c("A", "B"), each=250), dat)
    rm(dat.tx.a, dat.tx.b)
    dat <- reshape(dat, varying=c("measure.1", "measure.2", "measure.3"), 
                   idvar="subject.id", direction="long")
#ggplot
dat %>%
        group_by(tx,time) %>%
        summarize(mean=mean(measure), sd=sd(measure)) %>% 
        ggplot(aes(x=time, y=mean, color=tx))+
        geom_point()+
        geom_line(aes(group=1))+
        labs(x='Time', y='Cytokine Level')

这是我得到的情节。线条完全拧紧。情节怎么了?谢谢。

enter image description here

1 个答案:

答案 0 :(得分:0)

您需要将group=1替换为group=tx

dat %>%
  group_by(tx,time) %>%
  summarize(mean=mean(measure), sd=sd(measure)) %>% 
  ggplot(aes(x=time, y=mean, color=tx))+
  geom_point()+
  geom_line(aes(group=tx))+
  labs(x='Time', y='Cytokine Level')

enter image description here