Error: undefined columns selected when subsetting data based on max value

HABSaremyjam · August 13, 2021, 2:44pm

Hello!

I am trying to subset a dataset by separating them between the maximum value in the first column. One dataframe will be from the lowest to highest value (dat.down), and the other will be all the other values (dat.bottle). The code works with data.down, but for dat.bottle, I get the following error message:

Error in `[.data.frame`(dat.3, deepest:length(dat.3[, 1]), c("Depth (m)",  : 
  undefined columns selected

Can anyone help me figure out what's wrong? Here is my dataset and code:

dat.3 <- structure(list(DepFM = c(2, 2.25, 2.5, 2.75, 3, 3.25, 3.5, 3.75, 
4, 4.25, 4.5, 4.75, 5, 5.25, 5.5, 5.75, 6, 6.25, 6.5, 6.75), 
    T090C = c(24.1904, 24.1865, -9.99e-29, 24.1825, -9.99e-29, 
    24.1758, -9.99e-29, -9.99e-29, 24.1553, 24.1451, -9.99e-29, 
    -9.99e-29, 24.1367, -9.99e-29, -9.99e-29, 24.1013, 24.1111, 
    -9.99e-29, -9.99e-29, -9.99e-29), `Sbeox0Mg/L` = c(8.2004, 
    8.1958, -9.99e-29, 8.2201, -9.99e-29, 8.2201, -9.99e-29, 
    -9.99e-29, 8.223, 8.2315, -9.99e-29, -9.99e-29, 8.2377, -9.99e-29, 
    -9.99e-29, 8.2575, 8.2631, -9.99e-29, -9.99e-29, -9.99e-29
    ), `C0uS/cm` = c(268.275978, 268.263418, -9.99e-29, 268.281074, 
    -9.99e-29, 268.292471, -9.99e-29, -9.99e-29, 268.177553, 
    268.108738, -9.99e-29, -9.99e-29, 268.065145, -9.99e-29, 
    -9.99e-29, 267.888035, 267.942992, -9.99e-29, -9.99e-29, 
    -9.99e-29), Specc = c(272.659, 272.669, -9.99e-29, 272.709, 
    -9.99e-29, 272.757, -9.99e-29, -9.99e-29, 272.754, 272.741, 
    -9.99e-29, -9.99e-29, 272.743, -9.99e-29, -9.99e-29, 272.759, 
    272.76, -9.99e-29, -9.99e-29, -9.99e-29), Bat = c(0.1928, 
    0.1875, -9.99e-29, 0.1894, -9.99e-29, 0.1903, -9.99e-29, 
    -9.99e-29, 0.1964, 0.1952, -9.99e-29, -9.99e-29, 0.1989, 
    -9.99e-29, -9.99e-29, 0.2039, 0.2084, -9.99e-29, -9.99e-29, 
    -9.99e-29), Xmiss = c(95.2938, 95.4214, -9.99e-29, 95.3744, 
    -9.99e-29, 95.3548, -9.99e-29, -9.99e-29, 95.2078, 95.2365, 
    -9.99e-29, -9.99e-29, 95.149, -9.99e-29, -9.99e-29, 95.0313, 
    94.9223, -9.99e-29, -9.99e-29, -9.99e-29), Spar = c(2814.1, 
    2812.8, -9.99e-29, 2837.8, -9.99e-29, 2851.1, -9.99e-29, 
    -9.99e-29, 2856.6, 2858.4, -9.99e-29, -9.99e-29, 2854.8, 
    -9.99e-29, -9.99e-29, 2851.1, 2863, -9.99e-29, -9.99e-29, 
    -9.99e-29), Par = c(380.04, 455.91, -9.99e-29, 310.51, -9.99e-29, 
    433.91, -9.99e-29, -9.99e-29, 382.26, 271.45, -9.99e-29, 
    -9.99e-29, 264.25, -9.99e-29, -9.99e-29, 235.4, 228.16, -9.99e-29, 
    -9.99e-29, -9.99e-29), Cpar = c(13.505, 16.208, -9.99e-29, 
    10.947, -9.99e-29, 15.219, -9.99e-29, -9.99e-29, 13.381, 
    9.4966, -9.99e-29, -9.99e-29, 9.2564, -9.99e-29, -9.99e-29, 
    8.2562, 7.9691, -9.99e-29, -9.99e-29, -9.99e-29), Density00 = c(997.3594, 
    997.3615, -9.99e-29, 997.3649, -9.99e-29, 997.3688, -9.99e-29, 
    -9.99e-29, 997.3769, 997.3806, -9.99e-29, -9.99e-29, 997.3872, 
    -9.99e-29, -9.99e-29, 997.3977, 997.3977, -9.99e-29, -9.99e-29, 
    -9.99e-29), TimeS = c(60.642, 62.087, -9.99e-29, 82, -9.99e-29, 
    119, -9.99e-29, -9.99e-29, 122, 124.929, -9.99e-29, -9.99e-29, 
    126, -9.99e-29, -9.99e-29, 129, 130.854, -9.99e-29, -9.99e-29, 
    -9.99e-29), PrDM = c(1.961, 2.207, -9.99e-29, 2.66, -9.99e-29, 
    3.078, -9.99e-29, -9.99e-29, 3.804, 4.167, -9.99e-29, -9.99e-29, 
    4.796, -9.99e-29, -9.99e-29, 5.541, 5.884, -9.99e-29, -9.99e-29, 
    -9.99e-29), Nbf = c(0, 0, -9.99e-29, 0, -9.99e-29, 0, -9.99e-29, 
    -9.99e-29, 0, 0, -9.99e-29, -9.99e-29, 0, -9.99e-29, -9.99e-29, 
    0, 0, -9.99e-29, -9.99e-29, -9.99e-29), AltM = c(23.34, 26.1, 
    -9.99e-29, 35.1, -9.99e-29, 24.74, -9.99e-29, -9.99e-29, 
    39.41, 22.08, -9.99e-29, -9.99e-29, 21.86, -9.99e-29, -9.99e-29, 
    20.37, 19.34, -9.99e-29, -9.99e-29, -9.99e-29), WetStar = c(-0.1429, 
    -0.1435, -9.99e-29, -0.1444, -9.99e-29, -0.1416, -9.99e-29, 
    -9.99e-29, -0.1406, -0.1419, -9.99e-29, -9.99e-29, -0.1333, 
    -9.99e-29, -9.99e-29, -0.1411, -0.1351, -9.99e-29, -9.99e-29, 
    -9.99e-29), Ph = c(8.621, 8.634, -9.99e-29, 8.621, -9.99e-29, 
    8.622, -9.99e-29, -9.99e-29, 8.628, 8.625, -9.99e-29, -9.99e-29, 
    8.63, -9.99e-29, -9.99e-29, 8.629, 8.629, -9.99e-29, -9.99e-29, 
    -9.99e-29), Pumps = c(1, 1, -9.99e-29, 1, -9.99e-29, 1, -9.99e-29, 
    -9.99e-29, 1, 1, -9.99e-29, -9.99e-29, 1, -9.99e-29, -9.99e-29, 
    1, 1, -9.99e-29, -9.99e-29, -9.99e-29), Latitude = c(41.93497, 
    41.93496, -9.99e-29, 41.93495, -9.99e-29, 41.93488, -9.99e-29, 
    -9.99e-29, 41.93488, 41.93487, -9.99e-29, -9.99e-29, 41.93486, 
    -9.99e-29, -9.99e-29, 41.93486, 41.93486, -9.99e-29, -9.99e-29, 
    -9.99e-29), Longitude = c(-81.47777, -81.47776, -9.99e-29, 
    -81.47765, -9.99e-29, -81.47746, -9.99e-29, -9.99e-29, -81.47744, 
    -81.47744, -9.99e-29, -9.99e-29, -81.47744, -9.99e-29, -9.99e-29, 
    -81.47741, -81.4774, -9.99e-29, -9.99e-29, -9.99e-29), Sbeox0PS = c(97.828, 
    97.808, -9.99e-29, 97.94, -9.99e-29, 97.902, -9.99e-29, -9.99e-29, 
    98.077, 97.982, -9.99e-29, -9.99e-29, 98.101, -9.99e-29, 
    -9.99e-29, 98.438, 98.555, -9.99e-29, -9.99e-29, -9.99e-29
    ), Nbin = c(1, 1, 0, 3, 0, 1, 0, 0, 1, 1, 0, 0, 1, 0, 0, 
    1, 2, 0, 0, 0), Flag = c(0, 0, -9.99e-29, 0, -9.99e-29, 0, 
    -9.99e-29, -9.99e-29, 0, 0, -9.99e-29, -9.99e-29, 0, -9.99e-29, 
    -9.99e-29, 0, 0, -9.99e-29, -9.99e-29, -9.99e-29)), row.names = c(NA, 
20L), class = "data.frame")
> dput(head(dat.3,100))
structure(list(DepFM = c(2, 2.25, 2.5, 2.75, 3, 3.25, 3.5, 3.75, 
4, 4.25, 4.5, 4.75, 5, 5.25, 5.5, 5.75, 6, 6.25, 6.5, 6.75, 7, 
7.25, 7.5, 7.75, 8, 8.25, 8.5, 8.75, 9, 9.25, 9.5, 9.75, 10, 
10.25, 10.5, 10.75, 11, 11.25, 11.5, 11.75, 12, 12.25, 12.5, 
12.75, 13, 13.25, 13.5, 13.75, 14, 14.25, 14.5, 14.75, 15, 15.25, 
15.5, 15.75, 16, 16.25, 16.5, 16.75, 17, 17.25, 17.5, 17.75, 
18, 18.25, 18.5, 18.75, 19, 19.25, 19.5, 19.75, 20, 20.25, 20.5, 
20.75, 21, 21.25, 21.5, 21.75, 22, 22.25, 22.5, 22.75, 23, 23.25, 
23, 22.75, 22.5, 22.25, 22, 21.75, 21.5, 21.25, 21, 20.75, 20.5, 
20.25, 20, 19.75), T090C = c(24.1904, 24.1865, -9.99e-29, 24.1825, 
-9.99e-29, 24.1758, -9.99e-29, -9.99e-29, 24.1553, 24.1451, -9.99e-29, 
-9.99e-29, 24.1367, -9.99e-29, -9.99e-29, 24.1013, 24.1111, -9.99e-29, 
-9.99e-29, -9.99e-29, 24.0919, 24.0937, -9.99e-29, 24.0938, -9.99e-29, 
-9.99e-29, 24.0898, 24.0946, -9.99e-29, -9.99e-29, -9.99e-29, 
24.0801, 24.0776, 24.0769, -9.99e-29, -9.99e-29, 24.0775, 24.0771, 
-9.99e-29, -9.99e-29, -9.99e-29, 24.0755, -9.99e-29, 24.073, 
-9.99e-29, 24.0726, 24.0691, -9.99e-29, 24.0698, 24.0514, -9.99e-29, 
24.0435, 24.0424, 24.0291, -9.99e-29, 23.9533, 23.9418, -9.99e-29, 
-9.99e-29, 20.7539, -9.99e-29, -9.99e-29, 16.4656, 16.1806, -9.99e-29, 
-9.99e-29, 12.7271, -9.99e-29, 13.4159, -9.99e-29, 12.2821, -9.99e-29, 
-9.99e-29, -9.99e-29, 12.1374, -9.99e-29, -9.99e-29, -9.99e-29, 
11.9609, 12.1614, 12.2216, -9.99e-29, 11.8538, 11.7177, 11.7268, 
11.7339, -9.99e-29, 11.748, -9.99e-29, 11.7496, -9.99e-29, -9.99e-29, 
11.7485, -9.99e-29, 11.7507, -9.99e-29, 11.7534, -9.99e-29, 11.7645, 
-9.99e-29), `Sbeox0Mg/L` = c(8.2004, 8.1958, -9.99e-29, 8.2201, 
-9.99e-29, 8.2201, -9.99e-29, -9.99e-29, 8.223, 8.2315, -9.99e-29, 
-9.99e-29, 8.2377, -9.99e-29, -9.99e-29, 8.2575, 8.2631, -9.99e-29, 
-9.99e-29, -9.99e-29, 8.2633, 8.2581, -9.99e-29, 8.2568, -9.99e-29, 
-9.99e-29, 8.2592, 8.2685, -9.99e-29, -9.99e-29, -9.99e-29, 8.2673, 
8.2631, 8.2587, -9.99e-29, -9.99e-29, 8.2596, 8.2703, -9.99e-29, 
-9.99e-29, -9.99e-29, 8.262, -9.99e-29, 8.2587, -9.99e-29, 8.2632, 
8.2784, -9.99e-29, 8.2859, 8.2945, -9.99e-29, 8.3058, 8.3295, 
8.3561, -9.99e-29, 8.4079, 8.5271, -9.99e-29, -9.99e-29, 8.8425, 
-9.99e-29, -9.99e-29, 9.0432, 9.6955, -9.99e-29, -9.99e-29, 9.5009, 
-9.99e-29, 9.3609, -9.99e-29, 8.6813, -9.99e-29, -9.99e-29, -9.99e-29, 
7.8796, -9.99e-29, -9.99e-29, -9.99e-29, 7.1098, 6.8814, 6.4711, 
-9.99e-29, 6.2685, 5.6928, 5.7317, 5.7267, -9.99e-29, 5.7271, 
-9.99e-29, 5.723, -9.99e-29, -9.99e-29, 5.7233, -9.99e-29, 5.7256, 
-9.99e-29, 5.7283, -9.99e-29, 5.7347, -9.99e-29), `C0uS/cm` = c(268.275978, 
268.263418, -9.99e-29, 268.281074, -9.99e-29, 268.292471, -9.99e-29, 
-9.99e-29, 268.177553, 268.108738, -9.99e-29, -9.99e-29, 268.065145, 
-9.99e-29, -9.99e-29, 267.888035, 267.942992, -9.99e-29, -9.99e-29, 
-9.99e-29, 267.890368, 267.897179, -9.99e-29, 267.896721, -9.99e-29, 
-9.99e-29, 267.882745, 267.918966, -9.99e-29, -9.99e-29, -9.99e-29, 
267.936035, 267.94692, 267.94416, -9.99e-29, -9.99e-29, 267.943993, 
267.979819, -9.99e-29, -9.99e-29, -9.99e-29, 267.966752, -9.99e-29, 
267.981438, -9.99e-29, 267.985339, 267.984076, -9.99e-29, 267.98454, 
268.008614, -9.99e-29, 268.026186, 267.987947, 268.037955, -9.99e-29, 
267.906022, 267.835207, -9.99e-29, -9.99e-29, 251.359333, -9.99e-29, 
-9.99e-29, 228.470706, 225.8807, -9.99e-29, -9.99e-29, 209.472704, 
-9.99e-29, 213.522534, -9.99e-29, 207.789976, -9.99e-29, -9.99e-29, 
-9.99e-29, 207.218023, -9.99e-29, -9.99e-29, -9.99e-29, 206.299225, 
207.197947, 207.416024, -9.99e-29, 205.551725, 204.651953, 204.663459, 
204.674921, -9.99e-29, 204.746665, -9.99e-29, 204.748845, -9.99e-29, 
-9.99e-29, 204.738945, -9.99e-29, 204.746327, -9.99e-29, 204.761559, 
-9.99e-29, 204.824514, -9.99e-29), Specc = c(272.659, 272.669, 
-9.99e-29, 272.709, -9.99e-29, 272.757, -9.99e-29, -9.99e-29, 
272.754, 272.741, -9.99e-29, -9.99e-29, 272.743, -9.99e-29, -9.99e-29, 
272.759, 272.76, -9.99e-29, -9.99e-29, -9.99e-29, 272.813, 272.811, 
-9.99e-29, 272.809, -9.99e-29, -9.99e-29, 272.818, 272.828, -9.99e-29, 
-9.99e-29, -9.99e-29, 272.926, 272.95, 272.952, -9.99e-29, -9.99e-29, 
272.948, 272.986, -9.99e-29, -9.99e-29, -9.99e-29, 272.983, -9.99e-29, 
273.012, -9.99e-29, 273.018, 273.035, -9.99e-29, 273.032, 273.159, 
-9.99e-29, 273.222, 273.189, 273.315, -9.99e-29, 273.6, 273.593, 
-9.99e-29, -9.99e-29, 274.681, -9.99e-29, -9.99e-29, 275.434, 
274.335, -9.99e-29, -9.99e-29, 277.56, -9.99e-29, 277.913, -9.99e-29, 
278.695, -9.99e-29, -9.99e-29, -9.99e-29, 278.966, -9.99e-29, 
-9.99e-29, -9.99e-29, 279.05, 278.759, 278.601, -9.99e-29, 278.854, 
278.662, 278.608, 278.57, -9.99e-29, 278.561, -9.99e-29, 278.551, 
-9.99e-29, -9.99e-29, 278.546, -9.99e-29, 278.54, -9.99e-29, 
278.539, -9.99e-29, 278.541, -9.99e-29), Bat = c(0.1928, 0.1875, 
-9.99e-29, 0.1894, -9.99e-29, 0.1903, -9.99e-29, -9.99e-29, 0.1964, 
0.1952, -9.99e-29, -9.99e-29, 0.1989, -9.99e-29, -9.99e-29, 0.2039, 
0.2084, -9.99e-29, -9.99e-29, -9.99e-29, 0.2125, 0.2125, -9.99e-29, 
0.2125, -9.99e-29, -9.99e-29, 0.2138, 0.2129, -9.99e-29, -9.99e-29, 
-9.99e-29, 0.2138, 0.2089, 0.2097, -9.99e-29, -9.99e-29, 0.2094, 
0.2105, -9.99e-29, -9.99e-29, -9.99e-29, 0.2101, -9.99e-29, 0.2125, 
-9.99e-29, 0.2138, 0.2113, -9.99e-29, 0.21, 0.2126, -9.99e-29, 
0.2088, 0.2088, 0.2063, -9.99e-29, 0.2051, 0.199, -9.99e-29, 
-9.99e-29, 0.1915, -9.99e-29, -9.99e-29, 0.2212, 0.2461, -9.99e-29, 
-9.99e-29, 0.2487, -9.99e-29, 0.27, -9.99e-29, 0.2927, -9.99e-29, 
-9.99e-29, -9.99e-29, 0.4321, -9.99e-29, -9.99e-29, -9.99e-29, 
0.5059, 0.5001, 0.5043, -9.99e-29, 0.5539, 0.6179, 0.6151, 0.6137, 
-9.99e-29, 0.632, -9.99e-29, 0.6059, -9.99e-29, -9.99e-29, 0.6028, 
-9.99e-29, 0.5936, -9.99e-29, 0.5985, -9.99e-29, 0.5773, -9.99e-29
), Xmiss = c(95.2938, 95.4214, -9.99e-29, 95.3744, -9.99e-29, 
95.3548, -9.99e-29, -9.99e-29, 95.2078, 95.2365, -9.99e-29, -9.99e-29, 
95.149, -9.99e-29, -9.99e-29, 95.0313, 94.9223, -9.99e-29, -9.99e-29, 
-9.99e-29, 94.8254, 94.8254, -9.99e-29, 94.8254, -9.99e-29, -9.99e-29, 
94.796, 94.8163, -9.99e-29, -9.99e-29, -9.99e-29, 94.796, 94.9115, 
94.8924, -9.99e-29, -9.99e-29, 94.899, 94.8731, -9.99e-29, -9.99e-29, 
-9.99e-29, 94.8843, -9.99e-29, 94.8254, -9.99e-29, 94.796, 94.8545, 
-9.99e-29, 94.8843, 94.8228, -9.99e-29, 94.9137, 94.9137, 94.9729, 
-9.99e-29, 95.0019, 95.146, -9.99e-29, -9.99e-29, 95.3254, -9.99e-29, 
-9.99e-29, 94.6196, 94.0324, -9.99e-29, -9.99e-29, 93.9725, -9.99e-29, 
93.4725, -9.99e-29, 92.9431, -9.99e-29, -9.99e-29, -9.99e-29, 
89.7668, -9.99e-29, -9.99e-29, -9.99e-29, 88.1197, 88.2489, 88.1575, 
-9.99e-29, 87.0707, 85.6828, 85.7462, 85.7764, -9.99e-29, 85.3845, 
-9.99e-29, 85.9433, -9.99e-29, -9.99e-29, 86.0119, -9.99e-29, 
86.208, -9.99e-29, 86.1051, -9.99e-29, 86.5609, -9.99e-29), Spar = c(2814.1, 
2812.8, -9.99e-29, 2837.8, -9.99e-29, 2851.1, -9.99e-29, -9.99e-29, 
2856.6, 2858.4, -9.99e-29, -9.99e-29, 2854.8, -9.99e-29, -9.99e-29, 
2851.1, 2863, -9.99e-29, -9.99e-29, -9.99e-29, 2843.9, 2850.1, 
-9.99e-29, 2858.4, -9.99e-29, -9.99e-29, 2796.5, 2875.9, -9.99e-29, 
-9.99e-29, -9.99e-29, 2800.2, 2819.8, 2840, -9.99e-29, -9.99e-29, 
2823.8, 2862.3, -9.99e-29, -9.99e-29, -9.99e-29, 2808.4, -9.99e-29, 
2812.9, -9.99e-29, 2822, 2849.1, -9.99e-29, 2849.3, 2870.3, -9.99e-29, 
2849.3, 2854, 2860.3, -9.99e-29, 2858.4, 2877.6, -9.99e-29, -9.99e-29, 
2873, -9.99e-29, -9.99e-29, 2863.9, 2889.1, -9.99e-29, -9.99e-29, 
2840.2, -9.99e-29, 2820.2, -9.99e-29, 2812.9, -9.99e-29, -9.99e-29, 
-9.99e-29, 2856.6, -9.99e-29, -9.99e-29, -9.99e-29, 2873, 2867.9, 
2866.7, -9.99e-29, 2780.2, 1711, 1401.4, 1747, -9.99e-29, 1997.3, 
-9.99e-29, 2231.8, -9.99e-29, -9.99e-29, 2546.5, -9.99e-29, 2512.5, 
-9.99e-29, 2518.9, -9.99e-29, 2609, -9.99e-29), Par = c(380.04, 
455.91, -9.99e-29, 310.51, -9.99e-29, 433.91, -9.99e-29, -9.99e-29, 
382.26, 271.45, -9.99e-29, -9.99e-29, 264.25, -9.99e-29, -9.99e-29, 
235.4, 228.16, -9.99e-29, -9.99e-29, -9.99e-29, 207.92, 201.19, 
-9.99e-29, 171.12, -9.99e-29, -9.99e-29, 152.84, 165.13, -9.99e-29, 
-9.99e-29, -9.99e-29, 152.41, 147.42, 133.76, -9.99e-29, -9.99e-29, 
114.88, 130.93, -9.99e-29, -9.99e-29, -9.99e-29, 106.87, -9.99e-29, 
95.825, -9.99e-29, 93.147, 88.788, -9.99e-29, 81.987, 78.888, 
-9.99e-29, 76.154, 70.602, 68.534, -9.99e-29, 64.214, 62.06, 
-9.99e-29, -9.99e-29, 58.29, -9.99e-29, -9.99e-29, 51.419, 50.969, 
-9.99e-29, -9.99e-29, 41.498, -9.99e-29, 37.97, -9.99e-29, 37.003, 
-9.99e-29, -9.99e-29, -9.99e-29, 28.241, -9.99e-29, -9.99e-29, 
-9.99e-29, 19.928, 18.151, 16.049, -9.99e-29, 13.084, 6.0438, 
4.9621, 5.9757, -9.99e-29, 6.6927, -9.99e-29, 9.683, -9.99e-29, 
-9.99e-29, 13.536, -9.99e-29, 14.725, -9.99e-29, 18.552, -9.99e-29, 
21.56, -9.99e-29), Cpar = c(13.505, 16.208, -9.99e-29, 10.947, 
-9.99e-29, 15.219, -9.99e-29, -9.99e-29, 13.381, 9.4966, -9.99e-29, 
-9.99e-29, 9.2564, -9.99e-29, -9.99e-29, 8.2562, 7.9691, -9.99e-29, 
-9.99e-29, -9.99e-29, 7.3111, 7.0606, -9.99e-29, 5.9865, -9.99e-29, 
-9.99e-29, 5.4653, 5.744, -9.99e-29, -9.99e-29, -9.99e-29, 5.4428, 
5.228, 4.7111, -9.99e-29, -9.99e-29, 4.0679, 4.5773, -9.99e-29, 
-9.99e-29, -9.99e-29, 3.8057, -9.99e-29, 3.4066, -9.99e-29, 3.3007, 
3.1163, -9.99e-29, 2.8774, 2.7484, -9.99e-29, 2.6727, 2.4738, 
2.396, -9.99e-29, 2.2465, 2.1569, -9.99e-29, -9.99e-29, 2.0289, 
-9.99e-29, -9.99e-29, 1.7954, 1.7641, -9.99e-29, -9.99e-29, 1.4611, 
-9.99e-29, 1.3464, -9.99e-29, 1.3155, -9.99e-29, -9.99e-29, -9.99e-29, 
0.98905, -9.99e-29, -9.99e-29, -9.99e-29, 0.69363, 0.63286, 0.5598, 
-9.99e-29, 0.47096, 0.36482, 0.3568, 0.34768, -9.99e-29, 0.33509, 
-9.99e-29, 0.43613, -9.99e-29, -9.99e-29, 0.53161, -9.99e-29, 
0.58608, -9.99e-29, 0.73649, -9.99e-29, 0.82637, -9.99e-29), 
    Density00 = c(997.3594, 997.3615, -9.99e-29, 997.3649, -9.99e-29, 
    997.3688, -9.99e-29, -9.99e-29, 997.3769, 997.3806, -9.99e-29, 
    -9.99e-29, 997.3872, -9.99e-29, -9.99e-29, 997.3977, 997.3977, 
    -9.99e-29, -9.99e-29, -9.99e-29, 997.406, 997.4071, -9.99e-29, 
    997.4091, -9.99e-29, -9.99e-29, 997.4135, 997.4136, -9.99e-29, 
    -9.99e-29, -9.99e-29, 997.4217, 997.4233, 997.4246, -9.99e-29, 
    -9.99e-29, 997.4281, 997.4292, -9.99e-29, -9.99e-29, -9.99e-29, 
    997.4335, -9.99e-29, 997.437, -9.99e-29, 997.4385, 997.441, 
    -9.99e-29, 997.4423, 997.4482, -9.99e-29, 997.4536, 997.452, 
    997.4604, -9.99e-29, 997.4766, 997.4806, -9.99e-29, -9.99e-29, 
    998.1793, -9.99e-29, -9.99e-29, 999.1069, 999.1694, -9.99e-29, 
    -9.99e-29, 999.5808, -9.99e-29, 999.5551, -9.99e-29, 999.6948, 
    -9.99e-29, -9.99e-29, -9.99e-29, 999.6776, -9.99e-29, -9.99e-29, 
    -9.99e-29, 999.7082, 999.6851, 999.6761, -9.99e-29, 999.725, 
    999.7391, 999.74, 999.7402, -9.99e-29, 999.7357, -9.99e-29, 
    999.7337, -9.99e-29, -9.99e-29, 999.7307, -9.99e-29, 999.7282, 
    -9.99e-29, 999.7254, -9.99e-29, 999.7212, -9.99e-29), TimeS = c(60.642, 
    62.087, -9.99e-29, 82, -9.99e-29, 119, -9.99e-29, -9.99e-29, 
    122, 124.929, -9.99e-29, -9.99e-29, 126, -9.99e-29, -9.99e-29, 
    129, 130.854, -9.99e-29, -9.99e-29, -9.99e-29, 133, 134.466, 
    -9.99e-29, 137, -9.99e-29, -9.99e-29, 138, 140.767, -9.99e-29, 
    -9.99e-29, -9.99e-29, 142, 142.981, 145.177, -9.99e-29, -9.99e-29, 
    147.5, 151.905, -9.99e-29, -9.99e-29, -9.99e-29, 151.5, -9.99e-29, 
    155, -9.99e-29, 156, 157.989, -9.99e-29, 159, 160.043, -9.99e-29, 
    161, 162.259, 163.006, -9.99e-29, 164, 165.13, -9.99e-29, 
    -9.99e-29, 167, -9.99e-29, -9.99e-29, 168, 172.607, -9.99e-29, 
    -9.99e-29, 172, -9.99e-29, 173, -9.99e-29, 176, -9.99e-29, 
    -9.99e-29, -9.99e-29, 178.5, -9.99e-29, -9.99e-29, -9.99e-29, 
    181, 182.199, 184.91, -9.99e-29, 187.667, 211.399, 228.277, 
    243.663, -9.99e-29, 247, -9.99e-29, 269.5, -9.99e-29, -9.99e-29, 
    290, -9.99e-29, 292, -9.99e-29, 293.5, -9.99e-29, 296, -9.99e-29
    ), PrDM = c(1.961, 2.207, -9.99e-29, 2.66, -9.99e-29, 3.078, 
    -9.99e-29, -9.99e-29, 3.804, 4.167, -9.99e-29, -9.99e-29, 
    4.796, -9.99e-29, -9.99e-29, 5.541, 5.884, -9.99e-29, -9.99e-29, 
    -9.99e-29, 6.833, 7.11, -9.99e-29, 7.512, -9.99e-29, -9.99e-29, 
    8.365, 8.581, -9.99e-29, -9.99e-29, -9.99e-29, 9.604, 9.806, 
    10.052, -9.99e-29, -9.99e-29, 10.869, 11.033, -9.99e-29, 
    -9.99e-29, -9.99e-29, 11.965, -9.99e-29, 12.608, -9.99e-29, 
    12.888, 13.239, -9.99e-29, 13.62, 13.975, -9.99e-29, 14.4, 
    14.71, 14.955, -9.99e-29, 15.513, 15.691, -9.99e-29, -9.99e-29, 
    16.332, -9.99e-29, -9.99e-29, 17.131, 17.407, -9.99e-29, 
    -9.99e-29, 18.244, -9.99e-29, 18.583, -9.99e-29, 19.123, 
    -9.99e-29, -9.99e-29, -9.99e-29, 20.155, -9.99e-29, -9.99e-29, 
    -9.99e-29, 21.008, 21.33, 21.575, -9.99e-29, 22.073, 22.31, 
    22.555, 22.8, -9.99e-29, 22.194, -9.99e-29, 21.835, -9.99e-29, 
    -9.99e-29, 21.119, -9.99e-29, 20.652, -9.99e-29, 20.118, 
    -9.99e-29, 19.525, -9.99e-29), Nbf = c(0, 0, -9.99e-29, 0, 
    -9.99e-29, 0, -9.99e-29, -9.99e-29, 0, 0, -9.99e-29, -9.99e-29, 
    0, -9.99e-29, -9.99e-29, 0, 0, -9.99e-29, -9.99e-29, -9.99e-29, 
    0, 0, -9.99e-29, 0, -9.99e-29, -9.99e-29, 0, 0, -9.99e-29, 
    -9.99e-29, -9.99e-29, 0, 0, 0, -9.99e-29, -9.99e-29, 0, 0, 
    -9.99e-29, -9.99e-29, -9.99e-29, 0, -9.99e-29, 0, -9.99e-29, 
    0, 0, -9.99e-29, 0, 0, -9.99e-29, 0, 0, 0, -9.99e-29, 0, 
    0, -9.99e-29, -9.99e-29, 0, -9.99e-29, -9.99e-29, 0, 0, -9.99e-29, 
    -9.99e-29, 0, -9.99e-29, 0, -9.99e-29, 0, -9.99e-29, -9.99e-29, 
    -9.99e-29, 0, -9.99e-29, -9.99e-29, -9.99e-29, 0, 0, 0, -9.99e-29, 
    0, 1, 2, 2, -9.99e-29, 3, -9.99e-29, 4, -9.99e-29, -9.99e-29, 
    4, -9.99e-29, 4, -9.99e-29, 4, -9.99e-29, 4, -9.99e-29), 
    AltM = c(23.34, 26.1, -9.99e-29, 35.1, -9.99e-29, 24.74, 
    -9.99e-29, -9.99e-29, 39.41, 22.08, -9.99e-29, -9.99e-29, 
    21.86, -9.99e-29, -9.99e-29, 20.37, 19.34, -9.99e-29, -9.99e-29, 
    -9.99e-29, 19.19, 18.05, -9.99e-29, 17.05, -9.99e-29, -9.99e-29, 
    18.07, 15.69, -9.99e-29, -9.99e-29, -9.99e-29, 17.22, 16.66, 
    15.11, -9.99e-29, -9.99e-29, 14.25, 12.72, -9.99e-29, -9.99e-29, 
    -9.99e-29, 13.84, -9.99e-29, 12.72, -9.99e-29, 11.92, 11.02, 
    -9.99e-29, 10.77, 10.29, -9.99e-29, 9.91, 9.52, 9.28, -9.99e-29, 
    9.01, 8.64, -9.99e-29, -9.99e-29, 8.13, -9.99e-29, -9.99e-29, 
    62.98, -20.98, -9.99e-29, -9.99e-29, 98.24, -9.99e-29, 67.59, 
    -9.99e-29, 24.64, -9.99e-29, -9.99e-29, -9.99e-29, 98.69, 
    -9.99e-29, -9.99e-29, -9.99e-29, 98.66, 47.21, -4.51, -9.99e-29, 
    1.74, 1.2, 6.18, 90.41, -9.99e-29, 66.76, -9.99e-29, 1.41, 
    -9.99e-29, -9.99e-29, 2.02, -9.99e-29, 2.27, -9.99e-29, 25.75, 
    -9.99e-29, 98.53, -9.99e-29), WetStar = c(-0.1429, -0.1435, 
    -9.99e-29, -0.1444, -9.99e-29, -0.1416, -9.99e-29, -9.99e-29, 
    -0.1406, -0.1419, -9.99e-29, -9.99e-29, -0.1333, -9.99e-29, 
    -9.99e-29, -0.1411, -0.1351, -9.99e-29, -9.99e-29, -9.99e-29, 
    -0.1396, -0.1335, -9.99e-29, -0.1413, -9.99e-29, -9.99e-29, 
    -0.1414, -0.1318, -9.99e-29, -9.99e-29, -9.99e-29, -0.1401, 
    -0.1343, -0.1399, -9.99e-29, -9.99e-29, -0.1417, 0.1843, 
    -9.99e-29, -9.99e-29, -9.99e-29, 0.0182, -9.99e-29, -0.1408, 
    -9.99e-29, -0.1414, -0.1403, -9.99e-29, -0.1432, -0.1402, 
    -9.99e-29, -0.1444, -0.1236, -0.1538, -9.99e-29, -0.1182, 
    -0.1213, -9.99e-29, -9.99e-29, 0.2326, -9.99e-29, -9.99e-29, 
    0.0746, 1.231, -9.99e-29, -9.99e-29, 0.6825, -9.99e-29, -0.1204, 
    -9.99e-29, 0.9729, -9.99e-29, -9.99e-29, -9.99e-29, 0.938, 
    -9.99e-29, -9.99e-29, -9.99e-29, 1.4622, 1.5999, 1.8269, 
    -9.99e-29, 1.2641, 1.9561, 1.8266, 1.7412, -9.99e-29, 1.7433, 
    -9.99e-29, 1.7738, -9.99e-29, -9.99e-29, 1.614, -9.99e-29, 
    1.6857, -9.99e-29, 1.5912, -9.99e-29, 1.6504, -9.99e-29), 
    Ph = c(8.621, 8.634, -9.99e-29, 8.621, -9.99e-29, 8.622, 
    -9.99e-29, -9.99e-29, 8.628, 8.625, -9.99e-29, -9.99e-29, 
    8.63, -9.99e-29, -9.99e-29, 8.629, 8.629, -9.99e-29, -9.99e-29, 
    -9.99e-29, 8.63, 8.63, -9.99e-29, 8.627, -9.99e-29, -9.99e-29, 
    8.629, 8.628, -9.99e-29, -9.99e-29, -9.99e-29, 8.63, 8.63, 
    8.63, -9.99e-29, -9.99e-29, 8.63, 8.63, -9.99e-29, -9.99e-29, 
    -9.99e-29, 8.63, -9.99e-29, 8.63, -9.99e-29, 8.63, 8.63, 
    -9.99e-29, 8.63, 8.63, -9.99e-29, 8.63, 8.63, 8.629, -9.99e-29, 
    8.627, 8.626, -9.99e-29, -9.99e-29, 8.637, -9.99e-29, -9.99e-29, 
    8.633, 8.627, -9.99e-29, -9.99e-29, 8.598, -9.99e-29, 8.553, 
    -9.99e-29, 8.572, -9.99e-29, -9.99e-29, -9.99e-29, 8.49, 
    -9.99e-29, -9.99e-29, -9.99e-29, 8.453, 8.436, 8.429, -9.99e-29, 
    8.402, 8.262, 8.209, 8.149, -9.99e-29, 8.139, -9.99e-29, 
    8.104, -9.99e-29, -9.99e-29, 8.065, -9.99e-29, 8.06, -9.99e-29, 
    8.056, -9.99e-29, 8.035, -9.99e-29), Pumps = c(1, 1, -9.99e-29, 
    1, -9.99e-29, 1, -9.99e-29, -9.99e-29, 1, 1, -9.99e-29, -9.99e-29, 
    1, -9.99e-29, -9.99e-29, 1, 1, -9.99e-29, -9.99e-29, -9.99e-29, 
    1, 1, -9.99e-29, 1, -9.99e-29, -9.99e-29, 1, 1, -9.99e-29, 
    -9.99e-29, -9.99e-29, 1, 1, 1, -9.99e-29, -9.99e-29, 1, 1, 
    -9.99e-29, -9.99e-29, -9.99e-29, 1, -9.99e-29, 1, -9.99e-29, 
    1, 1, -9.99e-29, 1, 1, -9.99e-29, 1, 1, 1, -9.99e-29, 1, 
    1, -9.99e-29, -9.99e-29, 1, -9.99e-29, -9.99e-29, 1, 1, -9.99e-29, 
    -9.99e-29, 1, -9.99e-29, 1, -9.99e-29, 1, -9.99e-29, -9.99e-29, 
    -9.99e-29, 1, -9.99e-29, -9.99e-29, -9.99e-29, 1, 1, 1, -9.99e-29, 
    1, 1, 1, 1, -9.99e-29, 1, -9.99e-29, 1, -9.99e-29, -9.99e-29, 
    1, -9.99e-29, 1, -9.99e-29, 1, -9.99e-29, 1, -9.99e-29), 
    Latitude = c(41.93497, 41.93496, -9.99e-29, 41.93495, -9.99e-29, 
    41.93488, -9.99e-29, -9.99e-29, 41.93488, 41.93487, -9.99e-29, 
    -9.99e-29, 41.93486, -9.99e-29, -9.99e-29, 41.93486, 41.93486, 
    -9.99e-29, -9.99e-29, -9.99e-29, 41.93486, 41.93486, -9.99e-29, 
    41.93484, -9.99e-29, -9.99e-29, 41.93484, 41.93484, -9.99e-29, 
    -9.99e-29, -9.99e-29, 41.93484, 41.93484, 41.93483, -9.99e-29, 
    -9.99e-29, 41.93482, 41.93482, -9.99e-29, -9.99e-29, -9.99e-29, 
    41.93482, -9.99e-29, 41.9348, -9.99e-29, 41.9348, 41.9348, 
    -9.99e-29, 41.9348, 41.9348, -9.99e-29, 41.9348, 41.9348, 
    41.9348, -9.99e-29, 41.93478, 41.93478, -9.99e-29, -9.99e-29, 
    41.93478, -9.99e-29, -9.99e-29, 41.93478, 41.93478, -9.99e-29, 
    -9.99e-29, 41.93478, -9.99e-29, 41.93478, -9.99e-29, 41.93476, 
    -9.99e-29, -9.99e-29, -9.99e-29, 41.93476, -9.99e-29, -9.99e-29, 
    -9.99e-29, 41.93476, 41.93476, 41.93476, -9.99e-29, 41.93475, 
    41.93471, 41.9347, 41.9347, -9.99e-29, 41.9347, -9.99e-29, 
    41.93468, -9.99e-29, -9.99e-29, 41.93466, -9.99e-29, 41.93466, 
    -9.99e-29, 41.93466, -9.99e-29, 41.93466, -9.99e-29), Longitude = c(-81.47777, 
    -81.47776, -9.99e-29, -81.47765, -9.99e-29, -81.47746, -9.99e-29, 
    -9.99e-29, -81.47744, -81.47744, -9.99e-29, -9.99e-29, -81.47744, 
    -9.99e-29, -9.99e-29, -81.47741, -81.4774, -9.99e-29, -9.99e-29, 
    -9.99e-29, -81.4774, -81.47738, -9.99e-29, -81.47738, -9.99e-29, 
    -9.99e-29, -81.47738, -81.47736, -9.99e-29, -9.99e-29, -9.99e-29, 
    -81.47735, -81.47735, -81.47734, -9.99e-29, -9.99e-29, -81.47732, 
    -81.47732, -9.99e-29, -9.99e-29, -9.99e-29, -81.47732, -9.99e-29, 
    -81.47728, -9.99e-29, -81.47728, -81.47728, -9.99e-29, -81.47728, 
    -81.47728, -9.99e-29, -81.47726, -81.47726, -81.47725, -9.99e-29, 
    -81.47724, -81.47724, -9.99e-29, -9.99e-29, -81.47724, -9.99e-29, 
    -9.99e-29, -81.47722, -81.47722, -9.99e-29, -9.99e-29, -81.47722, 
    -9.99e-29, -81.47722, -9.99e-29, -81.4772, -9.99e-29, -9.99e-29, 
    -9.99e-29, -81.47718, -9.99e-29, -9.99e-29, -9.99e-29, -81.47718, 
    -81.47716, -81.47716, -9.99e-29, -81.47714, -81.47703, -81.47696, 
    -81.47692, -9.99e-29, -81.47692, -9.99e-29, -81.47682, -9.99e-29, 
    -9.99e-29, -81.47674, -9.99e-29, -81.47674, -9.99e-29, -81.47672, 
    -9.99e-29, -81.4767, -9.99e-29), Sbeox0PS = c(97.828, 97.808, 
    -9.99e-29, 97.94, -9.99e-29, 97.902, -9.99e-29, -9.99e-29, 
    98.077, 97.982, -9.99e-29, -9.99e-29, 98.101, -9.99e-29, 
    -9.99e-29, 98.438, 98.555, -9.99e-29, -9.99e-29, -9.99e-29, 
    98.501, 98.256, -9.99e-29, 98.266, -9.99e-29, -9.99e-29, 
    98.302, 98.623, -9.99e-29, -9.99e-29, -9.99e-29, 98.595, 
    98.282, 98.224, -9.99e-29, -9.99e-29, 98.266, 98.894, -9.99e-29, 
    -9.99e-29, -9.99e-29, 98.345, -9.99e-29, 98.147, -9.99e-29, 
    98.195, 98.603, -9.99e-29, 98.595, 98.55, -9.99e-29, 98.494, 
    98.504, 98.425, -9.99e-29, 98.435, 98.376, -9.99e-29, -9.99e-29, 
    98.526, -9.99e-29, -9.99e-29, 97.096, 104.394, -9.99e-29, 
    -9.99e-29, 97.311, -9.99e-29, 95.362, -9.99e-29, 85.347, 
    -9.99e-29, -9.99e-29, -9.99e-29, 74.218, -9.99e-29, -9.99e-29, 
    -9.99e-29, 65.012, 63.288, 59.134, -9.99e-29, 57.85, 52.591, 
    52.901, 52.866, -9.99e-29, 52.882, -9.99e-29, 52.827, -9.99e-29, 
    -9.99e-29, 52.842, -9.99e-29, 52.843, -9.99e-29, 52.88, -9.99e-29, 
    52.939, -9.99e-29), Nbin = c(1, 1, 0, 3, 0, 1, 0, 0, 1, 1, 
    0, 0, 1, 0, 0, 1, 2, 0, 0, 0, 1, 2, 0, 1, 0, 0, 1, 2, 0, 
    0, 0, 1, 1, 1, 0, 0, 2, 1, 0, 0, 0, 2, 0, 1, 0, 1, 1, 0, 
    1, 1, 0, 1, 1, 1, 0, 1, 2, 0, 0, 1, 0, 0, 1, 1, 0, 0, 1, 
    0, 1, 0, 1, 0, 0, 0, 2, 0, 0, 0, 1, 2, 2, 0, 3, 8, 10, 3, 
    0, 1, 0, 6, 0, 0, 3, 0, 1, 0, 2, 0, 1, 0), Flag = c(0, 0, 
    -9.99e-29, 0, -9.99e-29, 0, -9.99e-29, -9.99e-29, 0, 0, -9.99e-29, 
    -9.99e-29, 0, -9.99e-29, -9.99e-29, 0, 0, -9.99e-29, -9.99e-29, 
    -9.99e-29, 0, 0, -9.99e-29, 0, -9.99e-29, -9.99e-29, 0, 0, 
    -9.99e-29, -9.99e-29, -9.99e-29, 0, 0, 0, -9.99e-29, -9.99e-29, 
    0, 0, -9.99e-29, -9.99e-29, -9.99e-29, 0, -9.99e-29, 0, -9.99e-29, 
    0, 0, -9.99e-29, 0, 0, -9.99e-29, 0, 0, 0, -9.99e-29, 0, 
    0, -9.99e-29, -9.99e-29, 0, -9.99e-29, -9.99e-29, 0, 0, -9.99e-29, 
    -9.99e-29, 0, -9.99e-29, 0, -9.99e-29, 0, -9.99e-29, -9.99e-29, 
    -9.99e-29, 0, -9.99e-29, -9.99e-29, -9.99e-29, 0, 0, 0, -9.99e-29, 
    0, 0, 0, 0, -9.99e-29, 0, -9.99e-29, 0, -9.99e-29, -9.99e-29, 
    0, -9.99e-29, 0, -9.99e-29, 0, -9.99e-29, 0, -9.99e-29)), row.names = c(NA, 
100L), class = "data.frame")

# Determine which row is the deepest. This is where the downcast stops and the upcast starts.
  deepest <- which.max(dat.3[,1])
  
  # Subset a new data frame of only the downcast
  dat.down <- dat.3[1:deepest, ]
  
  # Subset a new data frame of only the depth and bottle number from the upcast.
  dat.bottle <- dat.3[deepest:length(dat.3[ ,1]), c("Depth (m)","Nbf")]

Thank you so much for your assistance!

dlsweet · August 13, 2021, 3:12pm

Hi, I took a look at your code and created your data sets. You have the right idea, but it looks like you don't have a column named "Depth (m)".

Could it be the column "DepFM"?

I got your code to work great by changing that column name.

HABSaremyjam · August 13, 2021, 3:56pm

Gah - of course it was that easy. Been staring at R too long, Thanks so much for being a second pair of eyes!

dlsweet · August 13, 2021, 4:08pm

Of course! We've all been there!

system · August 20, 2021, 4:09pm

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.