How to combine multiple plot into one

Question

i am trying to plot mean, standard deviation, and standard error of mean of my data

Data table look like

freq mean-minP mean-correctedminP sd-minP sd-correctedminP se-minP se-correctedminP
0.04 0.002253326 0.2570858 0.00281535 0.2299418 5.365941e-06 0.0004382595  
0.06 0.002215208 0.2555567 0.002724357 0.2261158 4.193085e-06 0.0003480172    
0.08 0.002229885 0.2580928 0.002688148 0.2257434 4.788916e-06 0.0004021602    
0.10 0.002259921 0.2614918 0.002692339 0.2255382 5.17569e-06 0.0004335695    
0.12 0.002252171 0.2615271 0.00266768 0.2239421 5.341592e-06 0.0004484073    
0.14 0.002279967 0.2620723 0.002773566 0.226358 5.726879e-06 0.0004673856    
0.16 0.00227563 0.2634223 0.002696895 0.2248959 5.812342e-06 0.000484695    
0.18 0.002325724 0.2680895 0.002727271 0.226968 6.08316e-06 0.0005062506    
0.20 0.002338423 0.2698884 0.002716278 0.2267581 6.19585e-06 0.0005172369

Code is

library(ggplot2)

x = read.table("mean-sd-se.txt", header =T)

as.data.frame.matrix(x) 

class(x)

jpeg('mean-minP.jpg')

ggplot(x, aes(freq, mean.minP)) +
  geom_point() +
  geom_errorbar(aes(ymin = mean.minP - se.minP, ymax = mean.minP + se.minP))

dev.off()

jpeg('mean-correctedminP.jpg')

ggplot(x, aes(freq, mean.correctedminP )) +
  geom_point() +
  geom_errorbar(aes(ymin = mean.correctedminP - se.correctedminP, ymax = mean.correctedminP + se.correctedminP))

dev.off()

Third graph table is

freq ben-minP sd se
0.04 2.272638e-05 2.446807e-05 2.446807e-06
0.06 2.913836e-05 2.863156e-05 2.863156e-06
0.08 3.554738e-05 3.428234e-05 3.428234e-06
0.10 2.748583e-05 2.544388e-05 2.544388e-06
0.12 2.699006e-05 2.334158e-05 2.334158e-06
0.14 2.688554e-05 2.607751e-05 2.607751e-06
0.16 2.791795e-05 2.749683e-05 2.749683e-06
0.18 3.163485e-05 3.127313e-05 3.127313e-06
0.20 2.905028e-05 3.089178e-05 3.089178e-06

I want to put all these plots into one with different colors

Possible duplicate of [Side-by-side plots with ggplot2](https://stackoverflow.com/questions/1249548/side-by-side-plots-with-ggplot2) — pieca, Jun 04 '18 at 15:35
Are you trying to overlay all three plots? That would be a hell of a lot more confusing even with different colors and transparency. Not to mention the scales would be different. — acylam, Jun 04 '18 at 15:54
@useR i actually want to do that but i would appreciate suggestion. i created side by side using the help mentioned — star, Jun 04 '18 at 15:56

moodymudskipper · Answer 1 · 2018-06-04T16:16:07.087

There you go, the trickiest part is shaping the data, I use tidyverse for this.

library(tidyverse)
tidy_df <- df %>%
  select(freq,mean.minP,mean.correctedminP,se.minP,se.correctedminP) %>%
  gather(,,-1) %>%
  separate(key,into = c("fun","var")) %>%
  spread(fun,value)

tidy_df
#    freq           var        mean           se
# 1  0.04 correctedminP 0.257085800 4.382595e-04
# 2  0.04          minP 0.002253326 5.365941e-06
# 3  0.06 correctedminP 0.255556700 3.480172e-04
# 4  0.06          minP 0.002215208 4.193085e-06
# 5  0.08 correctedminP 0.258092800 4.021602e-04
# 6  0.08          minP 0.002229885 4.788916e-06
# 7  0.10 correctedminP 0.261491800 4.335695e-04
# 8  0.10          minP 0.002259921 5.175690e-06
# 9  0.12 correctedminP 0.261527100 4.484073e-04
# 10 0.12          minP 0.002252171 5.341592e-06
# 11 0.14 correctedminP 0.262072300 4.673856e-04
# 12 0.14          minP 0.002279967 5.726879e-06
# 13 0.16 correctedminP 0.263422300 4.846950e-04
# 14 0.16          minP 0.002275630 5.812342e-06
# 15 0.18 correctedminP 0.268089500 5.062506e-04
# 16 0.18          minP 0.002325724 6.083160e-06
# 17 0.20 correctedminP 0.269888400 5.172369e-04
# 18 0.20          minP 0.002338423 6.195850e-06

ggplot(tidy_df, aes(freq, mean,color=factor(var))) +
  geom_point() +
  geom_errorbar(aes(ymin = mean - se, ymax = mean + se))

Sorry I can't upload images from this computer. Is that what you wanted ? Values are very different so it's not so readable.

If you want to put 2 scales on one graph it's not at all straightforward with ggplot2 (we can rescale the data though).

data

df <- read.table(text=
"freq mean-minP mean-correctedminP sd-minP sd-correctedminP se-minP se-correctedminP
0.04 0.002253326 0.2570858 0.00281535 0.2299418 5.365941e-06 0.0004382595  
0.06 0.002215208 0.2555567 0.002724357 0.2261158 4.193085e-06 0.0003480172    
0.08 0.002229885 0.2580928 0.002688148 0.2257434 4.788916e-06 0.0004021602    
0.10 0.002259921 0.2614918 0.002692339 0.2255382 5.17569e-06 0.0004335695    
0.12 0.002252171 0.2615271 0.00266768 0.2239421 5.341592e-06 0.0004484073    
0.14 0.002279967 0.2620723 0.002773566 0.226358 5.726879e-06 0.0004673856    
0.16 0.00227563 0.2634223 0.002696895 0.2248959 5.812342e-06 0.000484695    
0.18 0.002325724 0.2680895 0.002727271 0.226968 6.08316e-06 0.0005062506    
0.20 0.002338423 0.2698884 0.002716278 0.2267581 6.19585e-06 0.0005172369",
h=T)

How to combine multiple plot into one

1 Answers1