tests/testthat/test-backcrossPedigree.R

In rohan-shah/mpMap2: Genetic Analysis of Multi-Parent Recombinant Inbred Lines

context("Test backcross pedigree")
test_that("Test that genetic data is as expected",
	{
		map <- qtl::sim.map(len = 500, n.mar = 11, anchor.tel=TRUE, include.x=FALSE, eq.spacing=TRUE)
		pedigree <- backcrossPedigree(1000)
		cross <- simulateMPCross(map = map, pedigree = pedigree, mapFunction = haldane)
		proportions <- table(cross@geneticData[[1]]@finals)/ length(cross@geneticData[[1]]@finals)
		#Expect that there are only two genotypes
		expect_identical(length(proportions), 2L)
		#Expect that proportions are close to half. 
		expect_lt(abs(proportions[1] - 0.5), 0.04)
		expect_lt(abs(proportions[2] - 0.5), 0.04)
	})