R ggplot2 / ezPlot：對每個分位數的協變量繪制3x3 RM ANOVA設計並帶有誤差線

Question

首先，我的數據：

http://www.pastebin.ca/2599202 （我希望這不太方便，因為我無法創建合適的示例數據）

我基本上需要的是像我在這里所做的那樣的情節。

我為我的3個實驗組繪制了針對ias （連續因變量）的重復測量因子time （x軸，3個水平）。 我做了4次（對於我的特征量度MIHT每個分位數， MIHT ， miht.binned ）。

我必須承認我不是真正的R專業人士，而ggplot2手冊對我來說實在是ggplot2過分了。 我創建了情節ezPlot （從ezANOVA ），只有成功地做了一下布局調整ggplot2 ：

PlotIAS = ezPlot(                               
      data = MyData
      , dv = .(ias)
      , wid = .(id)
      , between = .(GROUP, miht.binned)
      , within = .(time)
      , x = .(time)
      , split = .(GROUP)
      , col   = .(miht.binned)
      , x_lab = 'time of measurement'
      , y_lab = 'IAS Score (Mean)'
      #, do_bars = FALSE
      , type = 3
    )


    PlotIAS = PlotIAS +  
      theme(
        panel.grid.major.y = element_line(colour = "gray80", size = NULL, linetype = NULL,  
                                          lineend = NULL)
        ,panel.grid.minor.y = element_line(colour = "gray90", size = NULL, linetype = NULL,
                                           lineend = NULL)
        ,panel.grid.major.x = element_blank()           
        ,panel.grid.minor.x = element_blank()
        ,legend.background = element_rect(fill = NULL, colour = "black") 

        ,panel.background = element_rect(fill = "white", colour = "white", size = NULL, 
                                         linetype = NULL)
      )

print(PlotIAS)

我沒有找到有關ezPlot創建的這些錯誤欄的任何信息。 它們對於每個點似乎都是相同的，並且可以使用bar_size =任意重新bar_size =長度。 我只需要擁有SE或CI的錯誤欄。 我不知道是否可以在基於ezPlot的代碼中添加這些代碼（以及如何添加），或者是否必須為此創建一個完整的新ggplot對象（這對我來說是一個過多的負擔...）。 非常感謝您的幫助。

Answer 1

我認為這很接近您想要的：

ggplot(MyData, aes(x=time, y=ias, colour=GROUP, group=GROUP, 
                   linetype=GROUP, shape=GROUP)) + 
  facet_grid(~miht.binned) +
  stat_summary(fun.data="mean_cl_boot", geom="errorbar", conf.int=90) + 
    #alternatives:
    #stat_summary(fun.data="mean_cl_normal", geom="errorbar") + 
    #stat_summary(fun.data="mean_sdl", geom="errorbar") + 
  stat_summary(fun.y="mean", geom="point", size=2) + 
  stat_summary(fun.y="mean", geom="line") + 
  theme_bw()

有關這些錯誤欄來自何處的更多信息，請參見?mean_cl_boot等。 conf.int是CI級別。 同樣，將所有三個顏色，線型和形狀都映射到GROUP似乎有點過分。 您可能沒有線型和形狀。

讓我補充一點

ggplot(MyData, aes(x=time, y=ias, fill=GROUP)) + facet_grid(~miht.binned) +
  geom_boxplot() + theme_bw()

實際上可能是一個更易於閱讀的圖（沒有交叉線/重疊線），同時保留了更多數據特征（最小/最大，離群值）。