tests/testthat/test-generateMissingSSRows.R

In ices-tools-dev/icesRDBES: Functions for the ICES Regional Database and Estimation System (RDBES)

capture.output({ ## suppresses printing of console output when running test()

  test_that("generateMissingSSRows does not add any SS rows if none are missing (H1)", {

     myH1RawObject <- importRDBESDataCSV(rdbesExtractPath = "./h1_v_20250211")

    # Only use a subset of the test data
    myH1RawObject <- filterRDBESDataObject(myH1RawObject, c("DEstratumName"), c("DE_stratum1_H1", "DE_stratum2_H1", "DE_stratum3_H1"))
    myH1RawObject <- findAndKillOrphans(myH1RawObject, verbose = FALSE)

    # Change all FO rows to just be Lan (to match SS which just has Lan)
    myH1RawObject[["FO"]]$FOcatReg <- "Lan"

    # Try to generate any missing SS rows - there shouldn't be any missing
    mySSAfter <- generateMissingSSRows(myH1RawObject,
                                           "ZW_1965_SpeciesList1",
                                           verbose = FALSE)
    SSbefore <- nrow(myH1RawObject[["SS"]])
    SSafter <- nrow(mySSAfter)

    expect_equal(SSafter,SSbefore)
  })

  test_that("generateMissingSSRows produces an error if there is no SL data", {

    myH1RawObject <- importRDBESDataCSV(rdbesExtractPath = "./h1_v_20250211")
    # Only use a subset of the test data
    myH1RawObject <- filterRDBESDataObject(myH1RawObject, c("DEstratumName"), c("DE_stratum1_H1", "DE_stratum2_H1", "DE_stratum3_H1"))
    myH1RawObject <- findAndKillOrphans(myH1RawObject, verbose = FALSE)

    tableToRemove <- "SL"
    # Get rid of data
    myH1RawObject <- myH1RawObject[!(sapply(myH1RawObject, is.null))]
    myDFs <- lapply(myH1RawObject, as.data.frame)
    myDFs[[tableToRemove]] <- NULL
    myH1RawObject2 <- suppressWarnings(createRDBESDataObject(myDFs, castToCorrectDataTypes = FALSE))

    expect_error(generateMissingSSRows(myH1RawObject2,
                                       "ZW_1965_SpeciesList",
                                       verbose = FALSE),
                 paste0(tableToRemove," data does not exist in the input data but it is required"))
  })

  test_that("generateMissingSSRows produces an error if there is no SS data", {

    myH1RawObject <- importRDBESDataCSV(rdbesExtractPath = "./h1_v_20250211")
    # Only use a subset of the test data
    myH1RawObject <- filterRDBESDataObject(myH1RawObject, c("DEstratumName"), c("DE_stratum1_H1", "DE_stratum2_H1", "DE_stratum3_H1"))
    myH1RawObject <- findAndKillOrphans(myH1RawObject, verbose = FALSE)

    tableToRemove <- "SS"
    # Get rid of data
    myH1RawObject <- myH1RawObject[!(sapply(myH1RawObject, is.null))]
    myDFs <- lapply(myH1RawObject, as.data.frame)
    myDFs[[tableToRemove]] <- NULL
    myH1RawObject2 <- suppressWarnings(createRDBESDataObject(myDFs, castToCorrectDataTypes = FALSE))

    expect_error(generateMissingSSRows(myH1RawObject2,
                                       "ZW_1965_SpeciesList",
                                       verbose = FALSE),
                 paste0(tableToRemove," data does not exist in the input data but it is required"))
  })

  test_that("generateMissingSSRows produces an error if there is no FO data", {

    myH1RawObject <- importRDBESDataCSV(rdbesExtractPath = "./h1_v_20250211")
    # Only use a subset of the test data
    myH1RawObject <- filterRDBESDataObject(myH1RawObject, c("DEstratumName"), c("DE_stratum1_H1", "DE_stratum2_H1", "DE_stratum3_H1"))
    myH1RawObject <- findAndKillOrphans(myH1RawObject, verbose = FALSE)

    tableToRemove <- "FO"
    # Get rid of data
    myH1RawObject <- myH1RawObject[!(sapply(myH1RawObject, is.null))]
    myDFs <- lapply(myH1RawObject, as.data.frame)
    myDFs[[tableToRemove]] <- NULL
    myH1RawObject2 <- suppressWarnings(createRDBESDataObject(myDFs, castToCorrectDataTypes = FALSE))

    expect_error(generateMissingSSRows(myH1RawObject2,
                                       "ZW_1965_SpeciesList",
                                       verbose = FALSE),
                 paste0(tableToRemove," data does not exist in the input data but it is required"))
  })

  test_that("generateMissingSSRows produces an error if a species list name that does not exist is provided", {

    myH1RawObject <- importRDBESDataCSV(rdbesExtractPath = "./h1_v_20250211")
    # Only use a subset of the test data
    myH1RawObject <- filterRDBESDataObject(myH1RawObject, c("DEstratumName"), c("DE_stratum1_H1", "DE_stratum2_H1", "DE_stratum3_H1"))
    myH1RawObject <- findAndKillOrphans(myH1RawObject, verbose = FALSE)

    expect_error(generateMissingSSRows(myH1RawObject,
                                       "Stupid name",
                                       verbose = FALSE),
                 "The requested species list name does not exist in the input data")
  })

  test_that("generateMissingSSRows produces an error if a species list name is provided that is not compatible with the Year/Country", {

    myH1RawObject <- importRDBESDataCSV(rdbesExtractPath = "./h1_v_20250211")
    # Only use a subset of the test data
    myH1RawObject <- filterRDBESDataObject(myH1RawObject, c("DEstratumName"), c("DE_stratum1_H1", "DE_stratum2_H1", "DE_stratum3_H1"))
    myH1RawObject <- findAndKillOrphans(myH1RawObject, verbose = FALSE)

    # Change the SL country - this shoudl cause a compatibility error
    myH1RawObject[["SL"]]$SLcou <- "TT"

    expect_error(generateMissingSSRows(myH1RawObject,
                                       "ZW_1965_SpeciesList1",
                                       verbose = FALSE),
                 "The requested species list is not compatible with the combination of SS country and year")
  })

  test_that("generateMissingSSRows runs correctly for FO Dis, when there are no missing SS Dis rows (H1)", {

    myH1RawObject <- importRDBESDataCSV(rdbesExtractPath = "./h1_v_20250211")
    # Only use a subset of the test data
    myH1RawObject <- filterRDBESDataObject(myH1RawObject, c("DEstratumName"), c("DE_stratum1_H1", "DE_stratum2_H1", "DE_stratum3_H1"))
    myH1RawObject <- findAndKillOrphans(myH1RawObject, verbose = FALSE)

    # Set FO to Dis
    myH1RawObject[["FO"]]$FOcatReg <- "Dis"
    # Ensure only SS Dis rows exists
    myH1RawObject[["SS"]]$SScatchFra <- "Dis"

    # Try to generate any missing SS rows
    mySSAfter <- generateMissingSSRows(myH1RawObject,
                                       "ZW_1965_SpeciesList1",
                                       verbose = FALSE)
    SSbefore <- nrow(myH1RawObject[["SS"]])
    catchFracBefore <- sort(unique(myH1RawObject[["SS"]]$SScatchFra))
    SSafter <- nrow(mySSAfter)
    catchFracAfter <- sort(unique(mySSAfter$SScatchFra))

    expect_equal(SSafter,SSbefore)
    expect_equal(catchFracAfter,catchFracBefore)

    # Take our new SS rows and add them back into our original object
    myH1RawObject[["SS"]] <- mySSAfter
    expect_silent(validateRDBESDataObject(myH1RawObject))
  })

  test_that("generateMissingSSRows runs correctly for FO Dis, generates missing SS Dis rows (H1)", {

    myH1RawObject <- importRDBESDataCSV(rdbesExtractPath = "./h1_v_20250211")
    # Only use a subset of the test data
    myH1RawObject <- filterRDBESDataObject(myH1RawObject, c("DEstratumName"), c("DE_stratum1_H1", "DE_stratum2_H1", "DE_stratum3_H1"))
    myH1RawObject <- findAndKillOrphans(myH1RawObject, verbose = FALSE)

    # Set FO to Dis
    myH1RawObject[["FO"]]$FOcatReg <- "Dis"
    # Ensure only SS Lan rows exists
    myH1RawObject[["SS"]]$SScatchFra <- "Lan"

    # Try to generate any missing SS rows
    # All FO rows now have catReg "Dis", but SS only has Lan rows - so we expect Dis
    # rows to be added to SS
    mySSAfter <- generateMissingSSRows(myH1RawObject,
                                       "ZW_1965_SpeciesList1",
                                       verbose = FALSE)
    SSbefore <- nrow(myH1RawObject[["SS"]])
    catchFracBefore <- sort(unique(myH1RawObject[["SS"]]$SScatchFra))
    SSafter <- nrow(mySSAfter)
    catchFracAfter <- sort(unique(mySSAfter$SScatchFra))

    # We expect the number of SS rows to double
    expect_equal(SSafter,2*SSbefore)
    # We expect to have Dis rows added to SS
    expect_equal(catchFracAfter,c("Dis",catchFracBefore))

    # Take our new SS rows and add them back into our original object
    myH1RawObject[["SS"]] <- mySSAfter
    expect_silent(validateRDBESDataObject(myH1RawObject))
  })

  test_that("generateMissingSSRows runs correctly for FO Lan, when there are no missing SS Lan rows (H1)", {

    myH1RawObject <- importRDBESDataCSV(rdbesExtractPath = "./h1_v_20250211")
    # Only use a subset of the test data
    myH1RawObject <- filterRDBESDataObject(myH1RawObject, c("DEstratumName"), c("DE_stratum1_H1", "DE_stratum2_H1", "DE_stratum3_H1"))
    myH1RawObject <- findAndKillOrphans(myH1RawObject, verbose = FALSE)

    # Set FO to Lan
    myH1RawObject[["FO"]]$FOcatReg <- "Lan"
    # Ensure only SS Lan rows exists
    myH1RawObject[["SS"]]$SScatchFra <- "Lan"

    # Try to generate any missing SS rows
    mySSAfter <- generateMissingSSRows(myH1RawObject,
                                       "ZW_1965_SpeciesList1",
                                       verbose = FALSE)
    SSbefore <- nrow(myH1RawObject[["SS"]])
    catchFracBefore <- sort(unique(myH1RawObject[["SS"]]$SScatchFra))
    SSafter <- nrow(mySSAfter)
    catchFracAfter <- sort(unique(mySSAfter$SScatchFra))

    expect_equal(SSafter,SSbefore)
    expect_equal(catchFracAfter,catchFracBefore)

    # Take our new SS rows and add them back into our original object
    myH1RawObject[["SS"]] <- mySSAfter
    expect_silent(validateRDBESDataObject(myH1RawObject))
  })

  test_that("generateMissingSSRows runs correctly for FO Lan, generates missing SS Lan rows (H1)", {

    myH1RawObject <- importRDBESDataCSV(rdbesExtractPath = "./h1_v_20250211")
    # Only use a subset of the test data
    myH1RawObject <- filterRDBESDataObject(myH1RawObject, c("DEstratumName"), c("DE_stratum1_H1", "DE_stratum2_H1", "DE_stratum3_H1"))
    myH1RawObject <- findAndKillOrphans(myH1RawObject, verbose = FALSE)

    # Set FO to Lan
    myH1RawObject[["FO"]]$FOcatReg <- "Lan"
    # Ensure only SS Dis rows exists
    myH1RawObject[["SS"]]$SScatchFra <- "Dis"

    # Try to generate any missing SS rows
    mySSAfter <- generateMissingSSRows(myH1RawObject,
                                       "ZW_1965_SpeciesList1",
                                       verbose = FALSE)
    SSbefore <- nrow(myH1RawObject[["SS"]])
    catchFracBefore <- sort(unique(myH1RawObject[["SS"]]$SScatchFra))
    SSafter <- nrow(mySSAfter)
    catchFracAfter <- sort(unique(mySSAfter$SScatchFra))

    # We expect the number of SS rows to double
    expect_equal(SSafter,2*SSbefore)
    # We expect to have Lan rows added to SS
    expect_equal(catchFracAfter,c(catchFracBefore,"Lan"))

    # Take our new SS rows and add them back into our original object
    myH1RawObject[["SS"]] <- mySSAfter
    expect_silent(validateRDBESDataObject(myH1RawObject))
  })

  test_that("generateMissingSSRows runs correctly for FO All, generates missing SS Lan rows (H1)", {

    myH1RawObject <- importRDBESDataCSV(rdbesExtractPath = "./h1_v_20250211")
    # Only use a subset of the test data
    myH1RawObject <- filterRDBESDataObject(myH1RawObject, c("DEstratumName"), c("DE_stratum1_H1", "DE_stratum2_H1", "DE_stratum3_H1"))
    myH1RawObject <- findAndKillOrphans(myH1RawObject, verbose = FALSE)

    # Set FO to Lan
    myH1RawObject[["FO"]]$FOcatReg <- "All"
    # Ensure only SS Dis rows exists
    myH1RawObject[["SS"]]$SScatchFra <- "Dis"

    # Try to generate any missing SS rows
    mySSAfter <- generateMissingSSRows(myH1RawObject,
                                       "ZW_1965_SpeciesList1",
                                       verbose = FALSE)
    SSbefore <- nrow(myH1RawObject[["SS"]])
    catchFracBefore <- sort(unique(myH1RawObject[["SS"]]$SScatchFra))
    SSafter <- nrow(mySSAfter)
    catchFracAfter <- sort(unique(mySSAfter$SScatchFra))

    # We expect the number of SS rows to double
    expect_equal(SSafter,2*SSbefore)
    # We expect to have Lan rows added to SS
    expect_equal(catchFracAfter,c(catchFracBefore,"Lan"))

    # Take our new SS rows and add them back into our original object
    myH1RawObject[["SS"]] <- mySSAfter
    expect_silent(validateRDBESDataObject(myH1RawObject))
  })

  test_that("generateMissingSSRows runs correctly for FO All, generates missing SS Dis rows (H1)", {

    myH1RawObject <- importRDBESDataCSV(rdbesExtractPath = "./h1_v_20250211")
    # Only use a subset of the test data
    myH1RawObject <- filterRDBESDataObject(myH1RawObject, c("DEstratumName"), c("DE_stratum1_H1", "DE_stratum2_H1", "DE_stratum3_H1"))
    myH1RawObject <- findAndKillOrphans(myH1RawObject, verbose = FALSE)

    # Set FO to Lan
    myH1RawObject[["FO"]]$FOcatReg <- "All"
    # Ensure only SS Dis rows exists
    myH1RawObject[["SS"]]$SScatchFra <- "Lan"

    # Try to generate any missing SS rows
    mySSAfter <- generateMissingSSRows(myH1RawObject,
                                       "ZW_1965_SpeciesList1",
                                       verbose = FALSE)
    SSbefore <- nrow(myH1RawObject[["SS"]])
    catchFracBefore <- sort(unique(myH1RawObject[["SS"]]$SScatchFra))
    SSafter <- nrow(mySSAfter)
    catchFracAfter <- sort(unique(mySSAfter$SScatchFra))

    # We expect the number of SS rows to double
    expect_equal(SSafter,2*SSbefore)
    # We expect to have Lan rows added to SS
    expect_equal(catchFracAfter,c("Dis",catchFracBefore))

    # Take our new SS rows and add them back into our original object
    myH1RawObject[["SS"]] <- mySSAfter
    expect_silent(validateRDBESDataObject(myH1RawObject))
  })

  test_that("generateMissingSSRows runs correctly for FO All, generates missing SS Dis rows (H1).  Add extra SS col (strict validation = FALSE)", {

    myH1RawObject <- importRDBESDataCSV(rdbesExtractPath = "./h1_v_20250211")
    # Only use a subset of the test data
    myH1RawObject <- filterRDBESDataObject(myH1RawObject, c("DEstratumName"), c("DE_stratum1_H1", "DE_stratum2_H1", "DE_stratum3_H1"))
    myH1RawObject <- findAndKillOrphans(myH1RawObject, verbose = FALSE)

    # Set FO to Lan
    myH1RawObject[["FO"]]$FOcatReg <- "All"
    # Ensure only SS Dis rows exists
    myH1RawObject[["SS"]]$SScatchFra <- "Lan"

    # Add a col to SS
    myH1RawObject[["SS"]]$SSextracCol <- "Test"

    # Try to generate any missing SS rows (strict = FALSE)
    mySSAfter <- generateMissingSSRows(myH1RawObject,
                                       "ZW_1965_SpeciesList1",
                                       verbose = FALSE,
                                       strict = FALSE)
    SSbefore <- nrow(myH1RawObject[["SS"]])
    catchFracBefore <- sort(unique(myH1RawObject[["SS"]]$SScatchFra))
    SSafter <- nrow(mySSAfter)
    catchFracAfter <- sort(unique(mySSAfter$SScatchFra))

    # We expect the number of SS rows to double
    expect_equal(SSafter,2*SSbefore)
    # We expect to have Lan rows added to SS
    expect_equal(catchFracAfter,c("Dis",catchFracBefore))

    # Take our new SS rows and add them back into our original object
    myH1RawObject[["SS"]] <- mySSAfter
    # strict = FALSE validation
    expect_silent(validateRDBESDataObject(myH1RawObject, strict = FALSE))
  })

  test_that("generateMissingSSRows runs correctly for FO All, generates missing SS Dis rows (H1).  Add extra FO col (strict validation = FALSE)", {

    myH1RawObject <- importRDBESDataCSV(rdbesExtractPath = "./h1_v_20250211")
    # Only use a subset of the test data
    myH1RawObject <- filterRDBESDataObject(myH1RawObject, c("DEstratumName"), c("DE_stratum1_H1", "DE_stratum2_H1", "DE_stratum3_H1"))
    myH1RawObject <- findAndKillOrphans(myH1RawObject, verbose = FALSE)

    # Set FO to Lan
    myH1RawObject[["FO"]]$FOcatReg <- "All"
    # Ensure only SS Dis rows exists
    myH1RawObject[["SS"]]$SScatchFra <- "Lan"

    # Add a col to FO
    myH1RawObject[["FO"]]$FOextracCol <- "Test"

    # Try to generate any missing SS rows (strict = FALSE)
    mySSAfter <- generateMissingSSRows(myH1RawObject,
                                       "ZW_1965_SpeciesList1",
                                       verbose = FALSE,
                                       strict = FALSE)
    SSbefore <- nrow(myH1RawObject[["SS"]])
    catchFracBefore <- sort(unique(myH1RawObject[["SS"]]$SScatchFra))
    SSafter <- nrow(mySSAfter)
    catchFracAfter <- sort(unique(mySSAfter$SScatchFra))

    # We expect the number of SS rows to double
    expect_equal(SSafter,2*SSbefore)
    # We expect to have Lan rows added to SS
    expect_equal(catchFracAfter,c("Dis",catchFracBefore))

    # Take our new SS rows and add them back into our original object
    myH1RawObject[["SS"]] <- mySSAfter
    # strict = FALSE validation
    expect_silent(validateRDBESDataObject(myH1RawObject, strict = FALSE))
  })

  test_that("generateMissingSSRows runs correctly for FO All, generates missing SS Dis rows (H1).  Remove non-essential FO col (strict validation = FALSE)", {

    myH1RawObject <- importRDBESDataCSV(rdbesExtractPath = "./h1_v_20250211")
    # Only use a subset of the test data
    myH1RawObject <- filterRDBESDataObject(myH1RawObject, c("DEstratumName"), c("DE_stratum1_H1", "DE_stratum2_H1", "DE_stratum3_H1"))
    myH1RawObject <- findAndKillOrphans(myH1RawObject, verbose = FALSE)

    # Set FO to Lan
    myH1RawObject[["FO"]]$FOcatReg <- "All"
    # Ensure only SS Dis rows exists
    myH1RawObject[["SS"]]$SScatchFra <- "Lan"

    # Remove FOval col
    myH1RawObject[["FO"]]$FOval <- NULL

    # Try to generate any missing SS rows (strict = FALSE)
    mySSAfter <- generateMissingSSRows(myH1RawObject,
                                       "ZW_1965_SpeciesList1",
                                       verbose = FALSE,
                                       strict = FALSE)
    SSbefore <- nrow(myH1RawObject[["SS"]])
    catchFracBefore <- sort(unique(myH1RawObject[["SS"]]$SScatchFra))
    SSafter <- nrow(mySSAfter)
    catchFracAfter <- sort(unique(mySSAfter$SScatchFra))

    # We expect the number of SS rows to double
    expect_equal(SSafter,2*SSbefore)
    # We expect to have Lan rows added to SS
    expect_equal(catchFracAfter,c("Dis",catchFracBefore))

    # Take our new SS rows and add them back into our original object
    myH1RawObject[["SS"]] <- mySSAfter
    # strict = FALSE validation
    expect_silent(validateRDBESDataObject(myH1RawObject, strict = FALSE))
  })


  test_that("generateMissingSSRows runs correctly for FO All, generates missing SS Dis rows (H1).  Add extra SL col (strict validation = FALSE)", {

    myH1RawObject <- importRDBESDataCSV(rdbesExtractPath = "./h1_v_20250211")
    # Only use a subset of the test data
    myH1RawObject <- filterRDBESDataObject(myH1RawObject, c("DEstratumName"), c("DE_stratum1_H1", "DE_stratum2_H1", "DE_stratum3_H1"))
    myH1RawObject <- findAndKillOrphans(myH1RawObject, verbose = FALSE)

    # Set FO to Lan
    myH1RawObject[["FO"]]$FOcatReg <- "All"
    # Ensure only SS Dis rows exists
    myH1RawObject[["SS"]]$SScatchFra <- "Lan"

    # Add a col to SL
    myH1RawObject[["SL"]]$SLextracCol <- "Test"

    # Try to generate any missing SS rows (strict = FALSE)
    mySSAfter <- generateMissingSSRows(myH1RawObject,
                                       "ZW_1965_SpeciesList1",
                                       verbose = FALSE,
                                       strict = FALSE)
    SSbefore <- nrow(myH1RawObject[["SS"]])
    catchFracBefore <- sort(unique(myH1RawObject[["SS"]]$SScatchFra))
    SSafter <- nrow(mySSAfter)
    catchFracAfter <- sort(unique(mySSAfter$SScatchFra))

    # We expect the number of SS rows to double
    expect_equal(SSafter,2*SSbefore)
    # We expect to have Lan rows added to SS
    expect_equal(catchFracAfter,c("Dis",catchFracBefore))

    # Take our new SS rows and add them back into our original object
    myH1RawObject[["SS"]] <- mySSAfter
    # strict = FALSE validation
    expect_silent(validateRDBESDataObject(myH1RawObject, strict = FALSE))
  })

  test_that("generateMissingSSRows runs correctly for FO All, when there are no missing SS Catch rows (H1)", {

    myH1RawObject <- importRDBESDataCSV(rdbesExtractPath = "./h1_v_20250211")
    # Only use a subset of the test data
    myH1RawObject <- filterRDBESDataObject(myH1RawObject, c("DEstratumName"), c("DE_stratum1_H1", "DE_stratum2_H1", "DE_stratum3_H1"))
    myH1RawObject <- findAndKillOrphans(myH1RawObject, verbose = FALSE)

    # Set FO to Lan
    myH1RawObject[["FO"]]$FOcatReg <- "All"
    # Ensure only SS Dis rows exists
    myH1RawObject[["SS"]]$SScatchFra <- "Catch"

    # Try to generate any missing SS rows
    mySSAfter <- generateMissingSSRows(myH1RawObject,
                                       "ZW_1965_SpeciesList1",
                                       verbose = FALSE)
    SSbefore <- nrow(myH1RawObject[["SS"]])
    catchFracBefore <- sort(unique(myH1RawObject[["SS"]]$SScatchFra))
    SSafter <- nrow(mySSAfter)
    catchFracAfter <- sort(unique(mySSAfter$SScatchFra))

    expect_equal(SSafter,SSbefore)
    expect_equal(catchFracAfter,catchFracBefore)

    # Take our new SS rows and add them back into our original object
    myH1RawObject[["SS"]] <- mySSAfter
    expect_silent(validateRDBESDataObject(myH1RawObject))
  })


  test_that("generateMissingSSRows can be used to create a validRDBESDataObject", {

    myH1RawObject <- importRDBESDataCSV(rdbesExtractPath = "./h1_v_20250211")
    # Only use a subset of the test data
    myH1RawObject <- filterRDBESDataObject(myH1RawObject, c("DEstratumName"), c("DE_stratum1_H1", "DE_stratum2_H1", "DE_stratum3_H1"))
    myH1RawObject <- findAndKillOrphans(myH1RawObject, verbose = FALSE)

    # Set FO to Lan
    myH1RawObject[["FO"]]$FOcatReg <- "All"
    # Ensure only SS Dis rows exists
    myH1RawObject[["SS"]]$SScatchFra <- "Lan"

    # Try to generate any missing SS rows
    mySSAfter <- generateMissingSSRows(myH1RawObject,
                                       "ZW_1965_SpeciesList1",
                                       verbose = FALSE)
    SSbefore <- nrow(myH1RawObject[["SS"]])
    catchFracBefore <- sort(unique(myH1RawObject[["SS"]]$SScatchFra))
    SSafter <- nrow(mySSAfter)
    catchFracAfter <- sort(unique(mySSAfter$SScatchFra))

    # We expect the number of SS rows to double
    expect_equal(SSafter,2*SSbefore)
    # We expect to have Lan rows added to SS
    expect_equal(catchFracAfter,c("Dis",catchFracBefore))

    # Take our new SS rows and add them back into our original object
    myH1RawObject[["SS"]] <- mySSAfter
    expect_silent(validateRDBESDataObject(myH1RawObject))
  })


}) ## end capture.output