Multiple plots in one function

tidyverse
1

Hello there,

I am trying to create a function with multiple ggplots in it, so that all the plot get arranged in a single window. However, when I try to do it, all I get is the last plot in my code, or I get the following error message:

Error in ggplot + geom_line(data = a, aes(x = DEPTH, y = TEMPERATURE,  : 
  non-numeric argument to binary operator

Can anyone help me figure out what I'm doing wrong?

Here is my code:

filt <- function(x) {
  filter(wq, SITE_ID == x)
}

all.plots <- function(a,x){
  library(ggplot2)
  library(gridExtra)
  p1 <- ggplot() + geom_line(data=a, aes(x=DEPTH, y=DO, color=CAST)) + labs(title = "Dissolved Oxygen (mg/L)", subtitle = x) +
    coord_cartesian(xlim=NULL, ylim=c(0,15))
  p2 <- ggplot() + geom_line(data=a, aes(x=DEPTH, y=PH, color=CAST)) + labs(title = "pH", subtitle = x) +
    coord_cartesian(xlim=NULL, ylim=c(0,14))
  p3 <- ggplot + geom_line(data=a, aes(x=DEPTH, y=TEMPERATURE, color=CAST)) + labs(title = "Temperature (C)", subtitle = x) +
    coord_cartesian(xlim=NULL, ylim=c(0,30))
  p4 <- ggplot + geom_line(data=a, aes(x=DEPTH, y=CONDUCTIVITY, color=CAST)) + labs(title = "Conductivity (uS)", subtitle = x) +
    coord_cartesian(xlim=NULL, ylim=c(200,400))
  grid.arrange(p1, p2, p3, p4, ncol=2)
}
glba1013 <- filt(x="GLBA15-1013")
all.plots(a=glba1013, x="GLBA15-1013")

Also, here is the dataset from which I am pulling:

wq <- structure(list(PUBLICATION_DATE = c(42933, 42933, 42933, 42933, 
42933, 42933, 42933, 42933, 42933, 42933, 42933, 42933, 42933, 
42933, 42933, 42933, 42933, 42933, 42933, 42933, 42933, 42933, 
42933, 42933, 42933, 42933, 42933, 42933, 42933, 42933, 42933, 
42933, 42933, 42933, 42933, 42933, 42933, 42933, 42933, 42933, 
42933, 42933, 42933, 42933, 42933, 42933, 42933, 42933, 42933, 
42933, 42933, 42933, 42933, 42933, 42933, 42933, 42933, 42933, 
42933, 42933, 42933, 42933, 42933, 42933, 42933, 42933, 42933, 
42933, 42933, 42933, 42933, 42933, 42933, 42933, 42933, 42933, 
42933, 42933, 42933, 42933, 42933, 42933, 42933, 42933, 42933, 
42933, 42933, 42933, 42933, 42933, 42933, 42933, 42933, 42933, 
42933, 42933, 42933, 42933, 42933, 42933), UID = c(169770, 173491, 
169770, 173491, 169770, 173491, 169770, 173491, 169770, 173491, 
169770, 173491, 169770, 173491, 169770, 173491, 169770, 173491, 
169770, 173491, 169770, 173491, 169770, 173491, 169770, 173491, 
169770, 173491, 169770, 173491, 169770, 173491, 169770, 173491, 
169011, 169011, 169011, 169011, 169011, 169011, 169011, 169542, 
169542, 169542, 169542, 169542, 169542, 169542, 169542, 169542, 
169542, 169542, 1000017, 1000017, 1000017, 1000017, 1000017, 
1000017, 1000017, 1000017, 1000017, 1000028, 1000028, 1000028, 
1000028, 1000028, 1000028, 1000028, 1000016, 1000016, 1000016, 
1000016, 1000016, 1000016, 1000016, 1000016, 1000016, 1000016, 
1000016, 171524, 171524, 171524, 171524, 171524, 171524, 171524, 
171524, 171524, 171524, 171524, 171524, 171524, 176960, 176960, 
176960, 176960, 176960, 176960, 176960, 176960), SITE_ID = c("GLBA15-1013", 
"GLBA15-1013", "GLBA15-1013", "GLBA15-1013", "GLBA15-1013", "GLBA15-1013", 
"GLBA15-1013", "GLBA15-1013", "GLBA15-1013", "GLBA15-1013", "GLBA15-1013", 
"GLBA15-1013", "GLBA15-1013", "GLBA15-1013", "GLBA15-1013", "GLBA15-1013", 
"GLBA15-1013", "GLBA15-1013", "GLBA15-1013", "GLBA15-1013", "GLBA15-1013", 
"GLBA15-1013", "GLBA15-1013", "GLBA15-1013", "GLBA15-1013", "GLBA15-1013", 
"GLBA15-1013", "GLBA15-1013", "GLBA15-1013", "GLBA15-1013", "GLBA15-1013", 
"GLBA15-1013", "GLBA15-1013", "GLBA15-1013", "GLBA15-1034", "GLBA15-1034", 
"GLBA15-1034", "GLBA15-1034", "GLBA15-1034", "GLBA15-1034", "GLBA15-1034", 
"GLBA15-1052", "GLBA15-1052", "GLBA15-1052", "GLBA15-1052", "GLBA15-1052", 
"GLBA15-1052", "GLBA15-1052", "GLBA15-1052", "GLBA15-1052", "GLBA15-1052", 
"GLBA15-1052", "GLBA15-1056", "GLBA15-1056", "GLBA15-1056", "GLBA15-1056", 
"GLBA15-1056", "GLBA15-1056", "GLBA15-1056", "GLBA15-1056", "GLBA15-1056", 
"GLBA15-1064", "GLBA15-1064", "GLBA15-1064", "GLBA15-1064", "GLBA15-1064", 
"GLBA15-1064", "GLBA15-1064", "GLBA15-3207", "GLBA15-3207", "GLBA15-3207", 
"GLBA15-3207", "GLBA15-3207", "GLBA15-3207", "GLBA15-3207", "GLBA15-3207", 
"GLBA15-3207", "GLBA15-3207", "GLBA15-3207", "GLBA15-3219", "GLBA15-3219", 
"GLBA15-3219", "GLBA15-3219", "GLBA15-3219", "GLBA15-3219", "GLBA15-3219", 
"GLBA15-3219", "GLBA15-3219", "GLBA15-3219", "GLBA15-3219", "GLBA15-3219", 
"GLBA15-3219", "GLBA15-3223", "GLBA15-3223", "GLBA15-3223", "GLBA15-3223", 
"GLBA15-3223", "GLBA15-3223", "GLBA15-3223", "GLBA15-3223"), 
    DATE_COL = c(42179, 42207, 42179, 42207, 42179, 42207, 42179, 
    42207, 42179, 42207, 42179, 42207, 42179, 42207, 42179, 42207, 
    42179, 42207, 42179, 42207, 42179, 42207, 42179, 42207, 42179, 
    42207, 42179, 42207, 42179, 42207, 42179, 42207, 42179, 42207, 
    42172, 42172, 42172, 42172, 42172, 42172, 42172, 42177, 42177, 
    42177, 42177, 42177, 42177, 42177, 42177, 42177, 42177, 42177, 
    42166, 42166, 42166, 42166, 42166, 42166, 42166, 42166, 42166, 
    42173, 42173, 42173, 42173, 42173, 42173, 42173, 42165, 42165, 
    42165, 42165, 42165, 42165, 42165, 42165, 42165, 42165, 42165, 
    42193, 42193, 42193, 42193, 42193, 42193, 42193, 42193, 42193, 
    42193, 42193, 42193, 42193, 42233, 42233, 42233, 42233, 42233, 
    42233, 42233, 42233), VISIT_NO = c(1, 2, 1, 2, 1, 2, 1, 2, 
    1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 
    2, 1, 2, 1, 2, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 
    1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 
    1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 
    1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1), STUDY = c("NCCA 2015", 
    "NCCA 2015", "NCCA 2015", "NCCA 2015", "NCCA 2015", "NCCA 2015", 
    "NCCA 2015", "NCCA 2015", "NCCA 2015", "NCCA 2015", "NCCA 2015", 
    "NCCA 2015", "NCCA 2015", "NCCA 2015", "NCCA 2015", "NCCA 2015", 
    "NCCA 2015", "NCCA 2015", "NCCA 2015", "NCCA 2015", "NCCA 2015", 
    "NCCA 2015", "NCCA 2015", "NCCA 2015", "NCCA 2015", "NCCA 2015", 
    "NCCA 2015", "NCCA 2015", "NCCA 2015", "NCCA 2015", "NCCA 2015", 
    "NCCA 2015", "NCCA 2015", "NCCA 2015", "NCCA 2015", "NCCA 2015", 
    "NCCA 2015", "NCCA 2015", "NCCA 2015", "NCCA 2015", "NCCA 2015", 
    "NCCA 2015", "NCCA 2015", "NCCA 2015", "NCCA 2015", "NCCA 2015", 
    "NCCA 2015", "NCCA 2015", "NCCA 2015", "NCCA 2015", "NCCA 2015", 
    "NCCA 2015", "NCCA 2015", "NCCA 2015", "NCCA 2015", "NCCA 2015", 
    "NCCA 2015", "NCCA 2015", "NCCA 2015", "NCCA 2015", "NCCA 2015", 
    "NCCA 2015", "NCCA 2015", "NCCA 2015", "NCCA 2015", "NCCA 2015", 
    "NCCA 2015", "NCCA 2015", "NCCA 2015", "NCCA 2015", "NCCA 2015", 
    "NCCA 2015", "NCCA 2015", "NCCA 2015", "NCCA 2015", "NCCA 2015", 
    "NCCA 2015", "NCCA 2015", "NCCA 2015", "NCCA 2015", "NCCA 2015", 
    "NCCA 2015", "NCCA 2015", "NCCA 2015", "NCCA 2015", "NCCA 2015", 
    "NCCA 2015", "NCCA 2015", "NCCA 2015", "NCCA 2015", "NCCA 2015", 
    "NCCA 2015", "NCCA 2015", "NCCA 2015", "NCCA 2015", "NCCA 2015", 
    "NCCA 2015", "NCCA 2015", "NCCA 2015", "NCCA 2015"), PSTL_CODE = c("PA", 
    "PA", "PA", "PA", "PA", "PA", "PA", "PA", "PA", "PA", "PA", 
    "PA", "PA", "PA", "PA", "PA", "PA", "PA", "PA", "PA", "PA", 
    "PA", "PA", "PA", "PA", "PA", "PA", "PA", "PA", "PA", "PA", 
    "PA", "PA", "PA", "MI", "MI", "MI", "MI", "MI", "MI", "MI", 
    "PA", "PA", "PA", "PA", "PA", "PA", "PA", "PA", "PA", "PA", 
    "PA", "OH", "OH", "OH", "OH", "OH", "OH", "OH", "OH", "OH", 
    "OH", "OH", "OH", "OH", "OH", "OH", "OH", "OH", "OH", "OH", 
    "OH", "OH", "OH", "OH", "OH", "OH", "OH", "OH", "NY", "NY", 
    "NY", "NY", "NY", "NY", "NY", "NY", "NY", "NY", "NY", "NY", 
    "NY", "OH", "OH", "OH", "OH", "OH", "OH", "OH", "OH"), SAMPLE_TYPE = c("HYDRO", 
    "HYDRO", "HYDRO", "HYDRO", "HYDRO", "HYDRO", "HYDRO", "HYDRO", 
    "HYDRO", "HYDRO", "HYDRO", "HYDRO", "HYDRO", "HYDRO", "HYDRO", 
    "HYDRO", "HYDRO", "HYDRO", "HYDRO", "HYDRO", "HYDRO", "HYDRO", 
    "HYDRO", "HYDRO", "HYDRO", "HYDRO", "HYDRO", "HYDRO", "HYDRO", 
    "HYDRO", "HYDRO", "HYDRO", "HYDRO", "HYDRO", "HYDRO", "HYDRO", 
    "HYDRO", "HYDRO", "HYDRO", "HYDRO", "HYDRO", "HYDRO", "HYDRO", 
    "HYDRO", "HYDRO", "HYDRO", "HYDRO", "HYDRO", "HYDRO", "HYDRO", 
    "HYDRO", "HYDRO", "HYDRO", "HYDRO", "HYDRO", "HYDRO", "HYDRO", 
    "HYDRO", "HYDRO", "HYDRO", "HYDRO", "HYDRO", "HYDRO", "HYDRO", 
    "HYDRO", "HYDRO", "HYDRO", "HYDRO", "HYDRO", "HYDRO", "HYDRO", 
    "HYDRO", "HYDRO", "HYDRO", "HYDRO", "HYDRO", "HYDRO", "HYDRO", 
    "HYDRO", "HYDRO", "HYDRO", "HYDRO", "HYDRO", "HYDRO", "HYDRO", 
    "HYDRO", "HYDRO", "HYDRO", "HYDRO", "HYDRO", "HYDRO", "HYDRO", 
    "HYDRO", "HYDRO", "HYDRO", "HYDRO", "HYDRO", "HYDRO", "HYDRO", 
    "HYDRO"), LINE = c(1, 1, 2, 2, 3, 3, 4, 4, 5, 5, 6, 6, 7, 
    7, 8, 8, 16, 24, 15, 23, 14, 22, 13, 21, 12, 20, 11, 19, 
    10, 18, 9, 17, 0, 0, 1, 2, 3, 6, 5, 4, 0, 1, 2, 3, 4, 5, 
    21, 20, 19, 18, 17, 0, 1, 2, 3, 4, 8, 7, 6, 5, 0, 1, 2, 3, 
    6, 5, 4, 0, 1, 2, 3, 4, 5, 10, 9, 8, 7, 6, 0, 1, 2, 3, 4, 
    5, 6, 22, 21, 20, 19, 18, 17, 0, 1, 2, 3, 4, 5, 6, 7, 8), 
    CAST = c("DOWNCAST", "DOWNCAST", "DOWNCAST", "DOWNCAST", 
    "DOWNCAST", "DOWNCAST", "DOWNCAST", "DOWNCAST", "DOWNCAST", 
    "DOWNCAST", "DOWNCAST", "DOWNCAST", "DOWNCAST", "DOWNCAST", 
    "DOWNCAST", "DOWNCAST", "UPCAST", "UPCAST", "UPCAST", "UPCAST", 
    "UPCAST", "UPCAST", "UPCAST", "UPCAST", "UPCAST", "UPCAST", 
    "UPCAST", "UPCAST", "UPCAST", "UPCAST", "UPCAST", "UPCAST", 
    NA, NA, "DOWNCAST", "DOWNCAST", "DOWNCAST", "UPCAST", "UPCAST", 
    "UPCAST", NA, "DOWNCAST", "DOWNCAST", "DOWNCAST", "DOWNCAST", 
    "DOWNCAST", "UPCAST", "UPCAST", "UPCAST", "UPCAST", "UPCAST", 
    NA, "DOWNCAST", "DOWNCAST", "DOWNCAST", "DOWNCAST", "UPCAST", 
    "UPCAST", "UPCAST", "UPCAST", NA, "DOWNCAST", "DOWNCAST", 
    "DOWNCAST", "UPCAST", "UPCAST", "UPCAST", NA, "DOWNCAST", 
    "DOWNCAST", "DOWNCAST", "DOWNCAST", "DOWNCAST", "UPCAST", 
    "UPCAST", "UPCAST", "UPCAST", "UPCAST", NA, "DOWNCAST", "DOWNCAST", 
    "DOWNCAST", "DOWNCAST", "DOWNCAST", "DOWNCAST", "UPCAST", 
    "UPCAST", "UPCAST", "UPCAST", "UPCAST", "UPCAST", NA, "DOWNCAST", 
    "DOWNCAST", "DOWNCAST", "DOWNCAST", "DOWNCAST", "DOWNCAST", 
    "DOWNCAST", "DOWNCAST"), COL_LOC = c("Surface", "Surface", 
    "Mid-water", "Mid-water", "Mid-water", "Mid-water", "Mid-water", 
    "Mid-water", "Mid-water", "Mid-water", "Mid-water", "Mid-water", 
    "Mid-water", "Mid-water", "Bottom", "Bottom", "Surface", 
    "Surface", "Mid-water", "Mid-water", "Mid-water", "Mid-water", 
    "Mid-water", "Mid-water", "Mid-water", "Mid-water", "Mid-water", 
    "Mid-water", "Mid-water", "Mid-water", "Bottom", "Bottom", 
    NA, NA, "Surface", "Mid-water", "Bottom", "Surface", "Mid-water", 
    "Bottom", NA, "Surface", "Mid-water", "Mid-water", "Mid-water", 
    "Bottom", "Surface", "Mid-water", "Mid-water", "Mid-water", 
    "Bottom", NA, "Surface", "Mid-water", "Mid-water", "Bottom", 
    "Surface", "Mid-water", "Mid-water", "Bottom", NA, "Surface", 
    "Mid-water", "Bottom", "Surface", "Mid-water", "Bottom", 
    NA, "Surface", "Mid-water", "Mid-water", "Mid-water", "Bottom", 
    "Surface", "Mid-water", "Mid-water", "Mid-water", "Bottom", 
    NA, "Surface", "Mid-water", "Mid-water", NA, "Mid-water", 
    "Bottom", "Surface", "Mid-water", "Mid-water", "Mid-water", 
    "Mid-water", "Bottom", NA, "Surface", "Mid-water", "Mid-water", 
    "Mid-water", "Mid-water", "Mid-water", "Mid-water", "Mid-water"
    ), DEPTH = c(0.1, 0.1, 0.5, 0.5, 1, 1, 2, 2, 3, 3, 4, 4, 
    5, 5, 5.2, 5.3, 0.1, 0.1, 0.5, 0.5, 1, 1, 2, 2, 3, 3, 4, 
    4, 5, 5, 5.2, 5.3, NA, NA, 0.1, 0.5, 0.8, 0.1, 0.5, 0.8, 
    NA, 0.1, 0.5, 1, 2, 2.8, 0.1, 0.5, 1, 2, 2.8, NA, 0.1, 0.5, 
    1, 1.5, 0.1, 0.5, 1, 1.5, NA, 0.1, 0.5, 1, 0.1, 0.5, 1, NA, 
    0.1, 0.5, 1, 2, 2.3, 0.1, 0.5, 1, 2, 2.3, NA, 0.1, 0.5, 1, 
    2, 3, 3.6, 0.1, 0.5, 1, 2, 3, 3.6, NA, 0.1, 0.5, 1, 2, 3, 
    4, 5, 6), CONDUCTIVITY = c(330, 297.8, 330.2, 297.8, 329.8, 
    297.8, 329.9, 297.7, 328.7, 297.7, 328, 297.7, 327.3, 297.8, 
    324.4, 297.8, 325.7, 297.9, 325.9, 297.9, 325.8, 297.9, 325, 
    297.8, 324.8, 297.7, 325.2, 297.7, 323.2, 305.8, 324.3, 297.4, 
    NA, NA, 708, 715, 733, 725, 735, 737, NA, 303.3, 302.9, 302.8, 
    302.9, 301.1, 302.9, 302.9, 303, 302.6, 301.2, NA, 28.1, 
    27.8, 27.8, 27.8, 28.2, 28.2, 28, 27.9, NA, 396, 396, 396, 
    395, 395, 395, NA, 538, 538, 538, 537, 531, 539, 539, 539, 
    538, 422, NA, 278.5, 278.6, 278.7, 278.9, 278.9, 278.7, 278.3, 
    278.2, 278.8, 278.8, 278.7, 278.7, NA, 317.3, 317.3, 317.4, 
    317.3, 317.3, 317.1, 319.2, 312), DO = c(9.17, 9.99, 9.16, 
    10.05, 9.19, 10.06, 9.3, 10.02, 9.01, 10.04, 8.93, 9.99, 
    8.67, 9.6, 8.1, 5.44, 9.19, 10.09, 9.19, 10.09, 9.19, 10.06, 
    9.19, 10.01, 8.84, 9.96, 8.63, 9.79, 8.12, 8.5, 8.05, 5.4, 
    NA, NA, 7.3, 7.3, 7.35, 7.36, 7.4, 7.43, NA, 9.09, 9.21, 
    9.28, 9.32, 10.03, 9.16, 9.17, 9.26, 9.36, 10.01, NA, 8.7, 
    8.8, 8.8, 8.8, 8.8, 8.8, 8.8, 8.8, NA, 5.6, 5.3, 5.3, 5.2, 
    5.2, 5.2, NA, 10.3, 10.5, 10.4, 10.4, 10.2, 10.6, 10.5, 10.5, 
    10.5, 10.3, NA, 9.77, 9.68, 9.58, 9.44, 9.44, 9.52, 9.87, 
    9.89, 9.52, 9.42, 9.5, 9.51, NA, 9.48, 9.5, 9.5, 9.53, 9.52, 
    9.45, 7.84, 7.53), LIGHT_AMB = c(1258, 1661, 1008, 1602, 
    1355, 1603, 1663, 1621, 1745, 1698, 1280, 1716, 1316, 1701, 
    1500, 1647, 1789, 1631, 1763, 1662, 1509, 1670, 1356, 1669, 
    1482, 1647, 1330, 1625, 1298, 1605, 1043, 1610, NA, NA, 264.3, 
    280.8, 291, 334.3, 317, 307.2, NA, 1548, 1621, 1670, 1586, 
    1628, 1633, 1662, 1647, 1610, 1644, NA, 1727, 1732, 1750, 
    1775, 1746, 1765, 1771, 1800, NA, 1742, 2092, 1262, 1106, 
    934.1, 1129, NA, 348.2, 332.8, 337.3, 332.4, 340.1, 785.8, 
    419.9, 391.2, 377.2, 360.9, NA, 1345, 1364, 1429, 1405, 1450, 
    1452, 1496, 1473, 1463, 1453, 1462, 1455, NA, 1984, 1994, 
    1993, 1975, 1962, 1958, 1963, 1944), LIGHT_UW = c(259, 845.1, 
    209.8, 470.4, 186.2, 178.4, 197.8, 218.2, 173.3, 98.6, 68.8, 
    42, 28.5, 15.7, 23.4, 8.9, 210.8, 843.7, 205.9, 820, 177.8, 
    461.8, 101, 226.8, 61.8, 110.8, 58.5, 40.8, 33.2, 10.6, 23.5, 
    8.6, NA, NA, 111.2, 13.35, 4.85, 106.7, 12.89, 4.193, NA, 
    171.4, 140.4, 153.9, 134, 97.76, 195.5, 144.5, 137.2, 124.2, 
    100.9, NA, 1217, 795, 666.4, 394.9, 1163, 1088, 675, 312, 
    NA, 181.2, 150.5, 0.7, 177.4, 11.62, 0.5, NA, NA, NA, NA, 
    0, 0, NA, NA, NA, 0, 0, NA, 647.1, 560.9, 479, 342.1, 188.2, 
    203.1, 634.3, 414.2, 350.5, 269.7, 200.4, 175.8, NA, 1422, 
    1028, 673.6, 307.4, 167.3, 87.36, 50.59, 28.18), NCCA_FLAG = c(NA_character_, 
    NA_character_, NA_character_, NA_character_, NA_character_, 
    NA_character_, NA_character_, NA_character_, NA_character_, 
    NA_character_, NA_character_, NA_character_, NA_character_, 
    NA_character_, NA_character_, NA_character_, NA_character_, 
    NA_character_, NA_character_, NA_character_, NA_character_, 
    NA_character_, NA_character_, NA_character_, NA_character_, 
    NA_character_, NA_character_, NA_character_, NA_character_, 
    NA_character_, NA_character_, NA_character_, NA_character_, 
    NA_character_, NA_character_, NA_character_, NA_character_, 
    NA_character_, NA_character_, NA_character_, NA_character_, 
    NA_character_, NA_character_, NA_character_, NA_character_, 
    NA_character_, NA_character_, NA_character_, NA_character_, 
    NA_character_, NA_character_, NA_character_, NA_character_, 
    NA_character_, NA_character_, NA_character_, NA_character_, 
    NA_character_, NA_character_, NA_character_, NA_character_, 
    NA_character_, NA_character_, NA_character_, NA_character_, 
    NA_character_, NA_character_, NA_character_, NA_character_, 
    NA_character_, NA_character_, NA_character_, NA_character_, 
    NA_character_, NA_character_, NA_character_, NA_character_, 
    NA_character_, NA_character_, NA_character_, NA_character_, 
    NA_character_, NA_character_, NA_character_, NA_character_, 
    NA_character_, NA_character_, NA_character_, NA_character_, 
    NA_character_, NA_character_, NA_character_, NA_character_, 
    NA_character_, NA_character_, NA_character_, NA_character_, 
    NA_character_, NA_character_, NA_character_), PH = c(8.3, 
    8.7, 8.3, 8.8, 8.3, 8.8, 8.3, 8.7, 8.3, 8.8, 8.3, 8.8, 8.2, 
    8.7, 8.1, 7.9, 8.3, 8.8, 8.3, 8.8, 8.3, 8.8, 8.3, 8.8, 8.2, 
    8.8, 8.2, 8.7, 8.1, 8.4, 8.1, 7.8, NA, NA, 8, 8, 8, 8, 8, 
    8, NA, 8.2, 8.2, 8.2, 8.2, 8.3, 8.2, 8.2, 8.2, 8.2, 8.3, 
    NA, 8.7, 8.7, 8.7, 8.7, 8.7, 8.7, 8.7, 8.7, NA, 7.8, 7.7, 
    7.6, 7.5, 7.5, 7.5, NA, 8.5, 8.5, 8.5, 8.4, 8.3, 8.5, 8.5, 
    8.4, 8.4, 8, NA, 8.4, 8.3, 8.3, 8.2, 8.3, 8.3, 8.3, 8.3, 
    8.2, 8.2, 8.3, 8.3, NA, 8.5, 8.5, 8.5, 8.5, 8.5, 8.5, 8.2, 
    8.2), SALINITY = c(NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
    NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
    NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
    NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
    NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
    NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
    NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA
    ), STATION_DEPTH = c(NA, NA, NA, NA, NA, NA, NA, NA, NA, 
    NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
    NA, NA, NA, NA, NA, NA, NA, NA, 5.7, 5.8, NA, NA, NA, NA, 
    NA, NA, 1.2, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 3.3, 
    NA, NA, NA, NA, NA, NA, NA, NA, 2, NA, NA, NA, NA, NA, NA, 
    1.5, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 2.8, NA, NA, 
    NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 4.1, NA, NA, NA, 
    NA, NA, NA, NA, NA), TEMPERATURE = c(21.3, 23.6, 21.4, 23.6, 
    21.2, 23.6, 21.1, 23.6, 21, 23.5, 21, 23.5, 21, 23.3, 20.7, 
    20.3, 21.4, 23.6, 21.4, 23.6, 21.4, 23.6, 21.1, 23.6, 21, 
    23.5, 21, 23.5, 20.7, 21.8, 20.7, 20.3, NA, NA, 23.5, 23.5, 
    23.5, 23.5, 23.4, 23.4, NA, 20.9, 20.6, 19.9, 19.5, 19.1, 
    21, 20.8, 20.1, 19.5, 19.1, NA, 23, 23, 22.9, 22.9, 23, 23, 
    22.9, 22.9, NA, 22.6, 22.5, 22.4, 22.6, 22.4, 22.4, NA, 21.5, 
    21.5, 21.5, 21.5, 21.5, 21.6, 21.6, 21.5, 21, 21.5, NA, 21, 
    21, 21, 20.9, 20.9, 21, 21, 21, 20.9, 20.9, 20.9, 21, NA, 
    25.6, 25.6, 25.6, 25.6, 25.6, 25.5, 24.6, 24.3), NCCA_COMMENT = c(NA_character_, 
    NA_character_, NA_character_, NA_character_, NA_character_, 
    NA_character_, NA_character_, NA_character_, NA_character_, 
    NA_character_, NA_character_, NA_character_, NA_character_, 
    NA_character_, NA_character_, NA_character_, NA_character_, 
    NA_character_, NA_character_, NA_character_, NA_character_, 
    NA_character_, NA_character_, NA_character_, NA_character_, 
    NA_character_, NA_character_, NA_character_, NA_character_, 
    NA_character_, NA_character_, NA_character_, NA_character_, 
    NA_character_, NA_character_, NA_character_, NA_character_, 
    NA_character_, NA_character_, NA_character_, NA_character_, 
    NA_character_, NA_character_, NA_character_, NA_character_, 
    NA_character_, NA_character_, NA_character_, NA_character_, 
    NA_character_, NA_character_, NA_character_, NA_character_, 
    NA_character_, NA_character_, NA_character_, NA_character_, 
    NA_character_, NA_character_, NA_character_, NA_character_, 
    NA_character_, NA_character_, NA_character_, NA_character_, 
    NA_character_, NA_character_, NA_character_, NA_character_, 
    NA_character_, NA_character_, NA_character_, NA_character_, 
    NA_character_, NA_character_, NA_character_, NA_character_, 
    NA_character_, NA_character_, NA_character_, NA_character_, 
    NA_character_, NA_character_, NA_character_, NA_character_, 
    NA_character_, NA_character_, NA_character_, NA_character_, 
    NA_character_, NA_character_, NA_character_, NA_character_, 
    NA_character_, NA_character_, NA_character_, NA_character_, 
    NA_character_, NA_character_, NA_character_), WATERBODY = c("Lake_Erie", 
    "Lake_Erie", "Lake_Erie", "Lake_Erie", "Lake_Erie", "Lake_Erie", 
    "Lake_Erie", "Lake_Erie", "Lake_Erie", "Lake_Erie", "Lake_Erie", 
    "Lake_Erie", "Lake_Erie", "Lake_Erie", "Lake_Erie", "Lake_Erie", 
    "Lake_Erie", "Lake_Erie", "Lake_Erie", "Lake_Erie", "Lake_Erie", 
    "Lake_Erie", "Lake_Erie", "Lake_Erie", "Lake_Erie", "Lake_Erie", 
    "Lake_Erie", "Lake_Erie", "Lake_Erie", "Lake_Erie", "Lake_Erie", 
    "Lake_Erie", "Lake_Erie", "Lake_Erie", "Lake_Erie", "Lake_Erie", 
    "Lake_Erie", "Lake_Erie", "Lake_Erie", "Lake_Erie", "Lake_Erie", 
    "Lake_Erie", "Lake_Erie", "Lake_Erie", "Lake_Erie", "Lake_Erie", 
    "Lake_Erie", "Lake_Erie", "Lake_Erie", "Lake_Erie", "Lake_Erie", 
    "Lake_Erie", "Lake_Erie", "Lake_Erie", "Lake_Erie", "Lake_Erie", 
    "Lake_Erie", "Lake_Erie", "Lake_Erie", "Lake_Erie", "Lake_Erie", 
    "Lake_Erie", "Lake_Erie", "Lake_Erie", "Lake_Erie", "Lake_Erie", 
    "Lake_Erie", "Lake_Erie", "Lake_Erie", "Lake_Erie", "Lake_Erie", 
    "Lake_Erie", "Lake_Erie", "Lake_Erie", "Lake_Erie", "Lake_Erie", 
    "Lake_Erie", "Lake_Erie", "Lake_Erie", "Lake_Erie", "Lake_Erie", 
    "Lake_Erie", "Lake_Erie", "Lake_Erie", "Lake_Erie", "Lake_Erie", 
    "Lake_Erie", "Lake_Erie", "Lake_Erie", "Lake_Erie", "Lake_Erie", 
    "Lake_Erie", "Lake_Erie", "Lake_Erie", "Lake_Erie", "Lake_Erie", 
    "Lake_Erie", "Lake_Erie", "Lake_Erie", "Lake_Erie")), row.names = c(NA, 
-100L), class = c("tbl_df", "tbl", "data.frame"))

Any help would be greatly appreciated! Thanks so much!

2

You can just use this format:

ggplot() +
  geom_point(...) +
  geom_point(...) +
  ...

In short, just need to use ggplot() once.

1 Like
3

to put multiple ggplots into a single grob, use a package such as cowplot or patchwork.

4

Thank you so much for answering! Unfortunately that doesn't work, because I actually am trying to create multiple plots with different coordinate systems :frowning: I did try it though, wondering why only my last plot shows up...

5

I think you made a typo in the last two ggplot calls.
I also omitted the record where TEMPERATURE was NA .
Maybe you have to do something more clever than that to correct this.

library(ggplot2)
library(gridExtra)
library(dplyr)
#> 
#> Attaching package: 'dplyr'
#> The following object is masked from 'package:gridExtra':
#> 
#>     combine
#> The following objects are masked from 'package:stats':
#> 
#>     filter, lag
#> The following objects are masked from 'package:base':
#> 
#>     intersect, setdiff, setequal, union

wq <- structure(list(PUBLICATION_DATE = c(42933, 42933, 42933, 42933, 
# HanOostdijk : omitted most lines from your specification
.....
    "Lake_Erie", "Lake_Erie", "Lake_Erie", "Lake_Erie")), row.names = c(NA, 
-100L), class = c("tbl_df", "tbl", "data.frame"))

filt <- function(x) {
  filter(wq, SITE_ID == x)  %>%
  filter(!is.na(TEMPERATURE)) 
  # %>%
  # select(DEPTH,DO,CAST,PH,TEMPERATURE,CONDUCTIVITY)
}

all.plots <- function(a,x){
  p1 <- ggplot() + geom_line(data=a, aes(x=DEPTH, y=DO, color=CAST)) + labs(title = "Dissolved Oxygen (mg/L)", subtitle = x) +
    coord_cartesian(xlim=NULL, ylim=c(0,15))
  p2 <- ggplot() + geom_line(data=a, aes(x=DEPTH, y=PH, color=CAST)) + labs(title = "pH", subtitle = x) +
    coord_cartesian(xlim=NULL, ylim=c(0,14))
  p3 <- ggplot() + geom_line(data=a, aes(x=DEPTH, y=TEMPERATURE, color=CAST)) + labs(title = "Temperature (C)", subtitle = x) +
    coord_cartesian(xlim=NULL, ylim=c(0,30))
  p4 <- ggplot() + geom_line(data=a, aes(x=DEPTH, y=CONDUCTIVITY, color=CAST)) + labs(title = "Conductivity (uS)", subtitle = x) +
    coord_cartesian(xlim=NULL, ylim=c(200,400))
  grid.arrange(p1, p2, p3, p4, ncol=2)
}
glba1013 <- filt(x="GLBA15-1013")
all.plots(a=glba1013, x="GLBA15-1013")

Created on 2021-07-30 by the reprex package (v2.0.0)

image
image719×494 9.41 KB

1 Like
6

You're amazing! Thank you SO much - the worst part is I can't seem to find where my typos were! Too much time in R and staring at multiple lines of code. Thank you again, you're a life saver!

7

The typo was

ggplot +

In stead of

ggplot() +
8

This topic was automatically closed 7 days after the last reply. New replies are no longer allowed.

If you have a query related to it or one of the replies, start a new topic and refer back with a link.