Sélection de la direction des barres d'erreur multiples dans un tracé de ligne

En utilisant ggplot2, comment pouvons-nous choisir de manière sélective la direction des barres d'erreur? J'ai un graphique à 4 lignes dont les barres d'erreur se chevauchent et je voudrais choisir manuellement celles qui ne devraient apparaître que SE ou SD positives ou négatives.Sélection de la direction des barres d'erreur multiples dans un tracé de ligne

Voici ce que j'ai jusqu'à présent. Je premier tracé sans barres d'erreur:

require("ggplot2") 
ggplot(data.long, 
     aes(y = values, x = within_2, group = interaction(Between, within_1), linetype = within_1)) + 
    stat_summary(geom="line", fun.y="mean") + 
    stat_summary(geom="point", fun.y="mean", aes(shape = Between)) + 
    scale_shape_manual(values = c(0, 17)) + 
    labs(x = "xlab", y = "% change")

Alors je résume les données pour obtenir SE et SD:

require("Hmisc") 
data_summary <- summarySE(data.long, measurevar="values", groupvars=c("within_1","within_2", "Between"))

Terrain avec barres d'erreur:

ggplot(data_summary, 
     aes(y = values, x = within_2, group = interaction(Between, within_1), linetype = within_1)) + 
    stat_summary(geom="line", fun.y="mean") + 
    stat_summary(geom="point", fun.y="mean", aes(shape = Between)) + 
    scale_shape_manual(values = c(0, 17)) + 
    labs(x = "xlab", y = "% change") + 
    geom_errorbar(aes(ymin = values-se, ymax = values+se), 
        width=.2)

Comment tracer uniquement des barres d'erreur supérieures ou inférieures sur différentes lignes pour améliorer les visuels de tracé? Here there's a solution qui semble être sur le chemin, mais je suis incapable de le reproduire dans mon intrigue 4 lignes.

Voici les données:

data.long <- structure(list(Between = structure(c(2L, 2L, 2L, 2L, 2L, 2L, 
2L, 2L, 2L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 
2L, 2L, 2L, 2L, 2L, 2L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 
2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 
1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 1L, 1L, 1L, 1L, 
1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 1L, 
1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 
2L, 2L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 
2L, 2L, 2L, 2L, 2L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 
2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 
1L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 1L, 1L, 1L, 1L, 1L, 
1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 1L, 1L, 
1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 
2L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 
2L, 2L, 2L, 2L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L), .Label = c("CT", 
"TT"), class = "factor"), within_1 = structure(c(2L, 2L, 2L, 
2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 
2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 
2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 
2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 
2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 
2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 
2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 1L, 
1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 
1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 
1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 
1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 
1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 
1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 
1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 
1L), .Label = c("N", "H"), class = "factor"), within_2 = c("1", 
"1", "1", "1", "1", "1", "1", "1", "1", "1", "1", "1", "1", "1", 
"1", "1", "1", "1", "1", "2", "2", "2", "2", "2", "2", "2", "2", 
"2", "2", "2", "2", "2", "2", "2", "2", "2", "2", "2", "3", "3", 
"3", "3", "3", "3", "3", "3", "3", "3", "3", "3", "3", "3", "3", 
"3", "3", "3", "3", "4", "4", "4", "4", "4", "4", "4", "4", "4", 
"4", "4", "4", "4", "4", "4", "4", "4", "4", "4", "5", "5", "5", 
"5", "5", "5", "5", "5", "5", "5", "5", "5", "5", "5", "5", "5", 
"5", "5", "5", "6", "6", "6", "6", "6", "6", "6", "6", "6", "6", 
"6", "6", "6", "6", "6", "6", "6", "6", "6", "1", "1", "1", "1", 
"1", "1", "1", "1", "1", "1", "1", "1", "1", "1", "1", "1", "1", 
"1", "1", "2", "2", "2", "2", "2", "2", "2", "2", "2", "2", "2", 
"2", "2", "2", "2", "2", "2", "2", "2", "3", "3", "3", "3", "3", 
"3", "3", "3", "3", "3", "3", "3", "3", "3", "3", "3", "3", "3", 
"3", "4", "4", "4", "4", "4", "4", "4", "4", "4", "4", "4", "4", 
"4", "4", "4", "4", "4", "4", "4", "5", "5", "5", "5", "5", "5", 
"5", "5", "5", "5", "5", "5", "5", "5", "5", "5", "5", "5", "5", 
"6", "6", "6", "6", "6", "6", "6", "6", "6", "6", "6", "6", "6", 
"6", "6", "6", "6", "6", "6"), values = c(-1.4210854715202e-14, 
0, 0, 0, 0, 0, 0, 0, 0, 0, 1.4210854715202e-14, 0, 0, -1.4210854715202e-14, 
0, 0, 0, 1.4210854715202e-14, 0, -10.6057583104295, 3.82366147137421, 
4.70566944216206, 0.187724131379127, -2.03535330384851, 2.2524979765733, 
-2.9727148766444, 2.6083885402083, -4.61360697651398, -6.66382102505303, 
-1.89755148846209, 1.93353333867543, -3.11983762126631, -1.4345951465109, 
1.24952019949389, -3.11650324404344, -6.24550613366193, 1.00162236282669, 
3.03857695742765, -13.5070034150553, 3.77730864186645, 2.93181434857114, 
1.24485299776214, -8.79147436639039, 1.12367346753739, -4.94643186841199, 
-1.27354690713364, -9.08208771404033, -2.29587905024601, -0.886178182878069, 
4.50975449028302, -1.15072253151024, 0.418168345464125, 1.73306387514543, 
-3.40039886806825, -0.619776296430885, -2.99138061415866, 0.00261916388753036, 
-13.1610198532011, 2.36046192675914, 3.38443628501334, -1.67355568218787, 
-12.1136956835012, 1.89806525735736, -5.61536574685965, -4.07184931448117, 
-10.0281915596217, -8.25975424526955, -2.23420107114791, 6.16533079297051, 
-3.48741869496304, -5.34701895515406, 0.442597528682782, -4.8531422721731, 
-6.71971984650138, 0.55633231142032, 4.90437903753043, -18.1078379702318, 
1.62222015009712, 2.13720100979083, -1.17309255429642, -12.78396026415, 
0.729862841243275, -8.85206491114279, -8.31362171333841, -14.7101795669535, 
-1.79048494987613, -1.49496076749276, 3.41185281739396, -0.838416866619227, 
-2.9855687917721, 0.981470329055682, -0.468618335177922, 1.9761737856014, 
-3.11961635987447, 0.0122101591331756, -18.3764354546495, -0.839562634101327, 
-0.909215798953383, -3.04660579184322, -15.527856021636, -0.211986565635385, 
-10.2270437446076, -8.11032597620989, -17.0205715878628, -6.69427611524303, 
-3.74938133414415, 3.32631128248181, -8.32452015678412, -7.4699177381065, 
2.44417772472276, -6.43767201678928, -6.74659357477346, 0.344765610464776, 
3.07591251998821, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 
0, 0, 0, 0, -1.68055460693036, 0.561182122946207, 6.04376427020478, 
0.52606038309024, -1.92041213359749, 2.84769535274846, -4.09729336005844, 
-1.38630866031168, -4.40117348808879, -4.80696362592231, 0.524936720222485, 
-0.225035184695145, -0.720841177736872, 0.378127221604487, -2.17560352140093, 
-4.09583471719559, -5.77744921646007, 0.610060458577692, 4.7775865643377, 
-5.0208861691041, -0.881834167779218, 2.26361403828709, 1.81046574431237, 
-5.14112669366628, 3.73205231339362, -5.08511489358583, -3.82150725680953, 
-9.70438917132918, 0.33075871579895, -0.867931790168782, 1.19904914200561, 
-1.29825049887336, -0.113067560520705, 0.400223130726687, -2.18141887602395, 
-1.79457229366905, -2.16359725691149, 4.00029213506218, -9.21970490935017, 
-1.61596728535645, 2.18558304627001, 1.29045467896481, -8.54718013614946, 
2.05156900518132, -7.72721674043164, -9.38186074302632, -10.2589110749438, 
-5.64651239627536, -1.51893967360546, 2.6898530761015, -3.21933963350499, 
-0.553289096780574, -0.718499298125963, -7.41824549357261, -5.36249002000503, 
-2.33002523562874, 8.12201643395576, -11.4677701393646, -6.97071057624974, 
-0.551803400141694, 2.85443044484461, -10.4245256849277, 0.514338249427183, 
-9.02993382055641, -10.5768230419275, -14.1344173532463, 0.11489370036935, 
0.00285183138591094, 4.31327888241417, -1.27799113923298, -0.914744528454023, 
0.127682993945172, -4.87142828001647, -1.01903492042325, -0.83963875573221, 
0.896505058481708, -14.8096831429394, -5.90071739612186, 1.63305204156767, 
3.08190814542317, -14.00192066133, 0.641853802476064, -11.3190314697168, 
-13.2537662623402, -18.806303898898, -7.0410921988907, -2.65318926213971, 
7.46825055830635, -5.51894689174026, -1.20587792674647, 0.974813882949832, 
-5.8574589151346, -5.92200551366662, -2.7927135341117, 2.24223786198247 
)), row.names = c(NA, -228L), .Names = c("Between", "within_1", 
"within_2", "values"), class = "data.frame")

Source

2017-05-26 AJMA

Il me semble que vous devrez le faire manuellement. Résumez la moyenne et le SD dans un data.frame puis calculez ymax et ymin avec eux. Ensuite, changez la colonne ymax à la valeur de la moyenne pour supprimer la barre d'erreur supérieure, ou viceversa pour ceux où vous voulez tracer uniquement la barre d'erreur inférieure. Ensuite, alimentez geom_errorbar avec ces colonnes. – Osdorp

peut-être cela peut vous aider:

data_summary <- Rmisc::summarySE(data.long, measurevar="values", 
           groupvars=c("within_1","within_2", "Between")) 


#Based on the solution you found. 
data_summary$min <- data_summary$values - ((data_summary$Between=='CT' & data_summary$within_1=='H')| 
              (data_summary$Between=='TT'& data_summary$within_1=='H'))* 
              data_summary$se 

data_summary$max <- data_summary$values + ((data_summary$Between=='CT' & data_summary$within_1=='N')| 
              (data_summary$Between=='TT'& data_summary$within_1=='N'))* 
              data_summary$se 

ggplot(data_summary, 
     aes(y = values, x = within_2, group = interaction(Between, within_1), linetype = within_1)) + 
    stat_summary(geom="line", fun.y="mean") + 
    stat_summary(geom="point", fun.y="mean", aes(shape = Between)) + 
    scale_shape_manual(values = c(0, 17)) + 
    labs(x = "xlab", y = "% change",title="New columns:min and max") + 
    geom_errorbar(aes(ymin = min, ymax = max), 
       width=.2)

Vous pouvez toujours obtenir le même résultat en utilisant deux geom_errorbar() subsetting pour chaque geom_bar les données appropriées:

ggplot(data_summary, 
    aes(y = values, x = within_2, group = interaction(Between, within_1), linetype = within_1)) + 
stat_summary(geom="line", fun.y="mean") + 
stat_summary(geom="point", fun.y="mean", aes(shape = Between)) + 
scale_shape_manual(values = c(0, 17)) + 
labs(x = "xlab", y = "% change",title="2 geom_errorbar") + 
geom_errorbar(data=with(data_summary,data_summary[which(Between=='CT'& within_1=='H' | 
                  Between=='TT'& within_1=='H'),]), 
      aes(ymin = values-se, ymax = values), 
      width=.2)+ 
geom_errorbar(data=with(data_summary,data_summary[which(Between=='CT'& within_1=='N' | 
                  Between=='TT'& within_1=='N'),]) 
      ,aes(ymin = values, ymax = values+se), 
      width=.2)

capuchons sur une face

Une solution possible pourrait consister à mettre width = 0 dans les options de geom_errorbar et génèrent les bouchons en utilisant geom_segment.

ggplot(data_summary, 
    aes(y = values, x = within_2, group = interaction(Between, within_1), linetype = within_1)) + 

stat_summary(geom="line", fun.y="mean") + 
stat_summary(geom="point", fun.y="mean", aes(shape = Between)) + 
scale_shape_manual(values = c(0, 17)) + 
labs(x = "xlab", y = "% change",title="2 geom_errorbar - One-side caps") + 

#errorbar without caps 

geom_errorbar(data=with(data_summary,data_summary[which(Between=='CT'& within_1=='H' | 
                  Between=='TT'& within_1=='H'),]), 
      aes(ymin = values-se, ymax = values),width=0)+ 
geom_errorbar(data=with(data_summary,data_summary[which(Between=='CT'& within_1=='N' | 
                  Between=='TT'& within_1=='N'),]), 
      aes(ymin = values, ymax = values+se),width=0)+ 

#geom_segment for caps 

geom_segment(data=with(data_summary,data_summary[which(Between=='CT'& within_1=='H' | 
                  Between=='TT'& within_1=='H'),]), 
       aes(y=values-se,yend=values-se,x= as.numeric(within_2)-0.1,xend= as.numeric(within_2)+0.1))+ 

    geom_segment(data=with(data_summary,data_summary[which(Between=='CT'& within_1=='N' | 
                  Between=='TT'& within_1=='N'),]), 
       aes(y=values+se,yend=values+se,x= as.numeric(within_2)-0.1,xend= as.numeric(within_2)+0.1))

Hope this helps!

Source

2017-05-26 10:45:38

Merci qui fait le travail. Une idée de la façon dont les limites des barres d'erreur peuvent être supprimées d'un côté? – AJMA

La définition de geom_errorbar (width = 0) fonctionne, mais les majuscules sont supprimées des deux bords. Comment pourrait-il être fait seulement sur un bord? – AJMA

@AJMA, je ne sais pas si c'est possible, mais j'ai édité ma réponse avec une solution de contournement possible qui utilise geom_segment. –

Sélection de la direction des barres d'erreur multiples dans un tracé de ligne

Répondre

Questions connexes