inst/tinytest/test_writeIBDs.R

### Test writeIBDs function.

## Define file locations.
SxMloc <- system.file("extdata/SxM", "SxM_geno.txt", package = "statgenIBD")
SxMmap <- system.file("extdata/SxM", "SxM_map.txt", package = "statgenIBD")

## Compute IBDs
SxMIBD <- calcIBD(popType = "DH", markerFile = SxMloc, mapFile = SxMmap)

## Check that input checks are working correctly.
expect_error(writeIBDs(1),
             "IBDprob should be an object of class IBDprob")
expect_error(writeIBDs(SxMIBD, outFile = 1),
             "outFile should be a single character string ending in .txt")
expect_error(writeIBDs(SxMIBD, outFile = "a/b.txt"),
             "No permission to write to")

## Define tempt output file.
ibdOut <- tempfile(fileext = ".txt")

expect_error(writeIBDs(SxMIBD, outFile = ibdOut, decimals = "a"),
             "decimals should be a single numerical value")
expect_error(writeIBDs(SxMIBD, outFile = ibdOut, minProb = 0.75),
             "minProb should be a numerical value between 0 and")

## There is very little that can actually be tested.
## Just checking silent execution.
expect_silent(writeIBDs(SxMIBD, outFile = ibdOut))
expect_silent(writeIBDs(SxMIBD, outFile = ibdOut, decimals = 4, minProb = 0.3))
expect_silent(writeIBDs(SxMIBD, outFile = ibdOut, compress = TRUE))

Try the statgenIBD package in your browser

Any scripts or data that you put into this service are public.

statgenIBD documentation built on May 29, 2024, 1:32 a.m.