2012-04-05 78 views
7

我試圖繪製不同類型的圖表(線圖和條形圖)彼此之下,它們都具有相同的軸:使用相同的x軸在xyplot下面繪製條形圖?

c1 <- ggplot(data, aes(date, TotalMutObs)) + stat_smooth(se = FALSE) + 
      geom_point() + 
      opts(axis.title.x = theme_blank()) + 
      ylab("Cumulative number of new mutations") 
c2 <- ggplot(data, aes(date, distance)) + stat_smooth(se = FALSE) + 
      geom_point() + 
      opts(axis.title.x = theme_blank()) + 
      ylab("Cumulative mean pairwise distance") 
c3 <- ggplot(data, aes(x = date, y = NbOfHorses)) + 
      geom_bar(stat = "identity") + 
      opts(axis.title.x = theme_blank()) + 
      ylab("Number of horses sampled") 

grid.arrange(c1, c2,c3) 

然而,在x軸上的日期不會排隊等候不同的情節。 enter image description here

下面是一些數據來嘗試一下:

date<-c("2003-03-13","2003-03-25","2003-03-26","2003-03-27","2003-03-28","2003-03-31","2003-04-01","2003-04-02","2003-04-04","2003-04-08","2003-04-09","2003-04-10","2003-04-11","2003-04-14","2003-04-15","2003-04-17","2003-04-19","2003-04-21","2003-04-22","2003-04-28","2003-05-08"); 

NbOfHorses<-c("1","2","1","3","4","5","4","3","3","3","3","4","2","4","1","2","4","1","2","1","2"); 

TotalMutObs<-c("20","30","58","72","140","165","204","230","250","286","302","327","346","388","393","414","443","444","462","467","485"); 

distance<-c("0.000693202","0.00073544","0.000855432","0.000506876","0.000720193","0.000708047","0.000835468","0.000812401","0.000803149","0.000839117","0.000842048","0.000856393","0.000879973","0.000962382","0.000990666","0.001104861","0.001137515","0.001143838","0.00121874","0.001213737","0.001201379"); 

data<-as.data.frame(cbind(date,NbOfHorses,TotalMutObs,distance)); 

乾杯, 約瑟夫

+0

是不是有一個'XLIM()'在GGPLOT2功能?像'c3 < - c3 + xlim(範圍(數據$日期))' – 2012-04-05 14:26:25

+1

'我不認爲這是問題,因爲所有三個x軸跨度完全相同的範圍。 – blJOg 2012-04-05 14:38:02

+0

提供一些數據,我們會對其進行測試。 – 2012-04-05 14:44:33

回答

7

解決這個問題的方法是內GGPLOT2工作,並取得創造性地堆疊的副本你數據,然後將子集發送到您需要的每個geom。

#A version of your data cleaned up 
dat <- data.frame(date = as.Date(date),NbOfHorses = as.numeric(NbOfHorses), 
       TotalMutObs = as.numeric(TotalMutObs),distance = as.numeric(distance)) 

#Create three copies, one for each panel 
# Use informative titles for grp to be panel titles  
fullDat <- rbind(dat,dat,dat) 
fullDat$grp <- rep(c('Cumulative number of new mutations', 
         'Cumulative mean pairwise distance', 
         'Number of horses sampled'),each = nrow(dat)) 

ggplot(fullDat,aes(x = date)) + 
    facet_wrap(~grp,nrow = 3,scale = "free_y") + 
    geom_point(data = subset(fullDat,grp == 'Cumulative number of new mutations'), 
     aes(y = TotalMutObs)) + 
    stat_smooth(data = subset(fullDat,grp == 'Cumulative number of new mutations'), 
     aes(y = TotalMutObs),se = FALSE) + 
    geom_point(data = subset(fullDat,grp == 'Cumulative mean pairwise distance'), 
     aes(y = distance)) + 
    stat_smooth(data = subset(fullDat,grp == 'Cumulative mean pairwise distance'), 
     aes(y = distance),se = FALSE) + 
    geom_bar(data = subset(fullDat,grp == 'Number of horses sampled'), 
     aes(y = NbOfHorses),stat = "identity") + 
    labs(x = NULL,y = NULL) 

enter image description here

+0

這很美,非常感謝! – blJOg 2012-04-10 09:08:24

相關問題