R中帶有3個分類和1個連續變量的SE的點圖

Question

我正在嘗試生成一個點圖，該點圖在具有三個分類變量（均值，tre和性）的設計中具有單個度量值（len）的特征。

我制作了一個圖，該圖具有我想要的全部內容，可分為六個不同的子圖：

但我希望盡可能將它們全部放在一個圖上，最好使用ggplot。

這是我當前的R代碼：

ggplot(mydf, aes(x=factor(mea), y=len), group=sex) +
  geom_point() + geom_errorbar(limits, width=0.1) + facet_wrap(~ tre + sex)

以及來自dput的樣本數據：

structure(list(mea = structure(c(1L, 1L, 1L, 1L, 1L, 1L, 
2L, 2L, 2L, 2L, 2L, 2L), .Label = c("PO_P", "Melaniz"), class = "factor"), 
    tre = structure(c(1L, 1L, 2L, 2L, 3L, 3L, 1L, 1L, 2L, 
    2L, 3L, 3L), .Label = c("a", "b", "c"), class = "factor"), 
    Sex = structure(c(1L, 2L, 1L, 2L, 1L, 2L, 1L, 2L, 1L, 2L, 
    1L, 2L), .Label = c("Male", "Female"), class = "factor"), 
    N = c(26, 26, 25, 25, 27, 27, 14, 13, 12, 11, 14, 13), len = c(10.6615384615385, 
    10.5807692307692, 10.292, 10.6, 10.2851851851852, 10.6518518518519, 
    11.4785714285714, 11.7153846153846, 11.7083333333333, 11.5, 
    11.6214285714286, 11.8923076923077), sd = c(0.869057845290829, 
    0.779753412698774, 0.722218803410712, 0.654471797202395, 
    0.906686148609193, 0.8040141456708, 1.0123685642542, 0.805032249712347, 
    1.13654846981659, 0.822192191643779, 0.833139171519908, 0.739889111580849
    ), se = c(0.170436265829955, 0.152922225659293, 0.144443760682142, 
    0.130894359440479, 0.174491830656674, 0.154732594478434, 
    0.270566879755675, 0.223275773441538, 0.328093282497832, 
    0.247900273203854, 0.222665809666299, 0.205208317689404), 
    ci = c(0.351020060264102, 0.314949219318153, 0.298117269908016, 
    0.270152680174426, 0.358673094717481, 0.318057403068012, 
    0.584524206501098, 0.486476119728297, 0.722128445903482, 
    0.552356230143519, 0.481040236068982, 0.447110515336101)), .Names = c("mea", "tre", "sex", "N", "len", "sd", "se", "ci"), row.names = c(NA, 
-12L), class = "data.frame")

Answer 1

為了您可以使用顏色，形狀，大小，pointtype所有這些分類變量之間區分，等等。這是一個使用的顏色和點型為例sex和tre ，

library(ggplot2)
limits <- aes(ymax=mydf$len+mydf$se, ymin=mydf$len-mydf$se)

ggplot(mydf, aes(x=factor(mea), y=len, color=sex, pch=tre)) +
  geom_point(position=position_dodge(width=0.5)) +
  geom_errorbar(limits, position=position_dodge(width=0.5))