tests/testthat/test_plotting.R

In bomeara/treevo: Using ABC to Understand Trait Evolution

test_that("plotABC_3D works", {

	if (requireNamespace("gpclib", quietly = TRUE) 
		& requireNamespace("rgl", quietly = TRUE
		)) {

		set.seed(1)		
		data(simRunExample)
		
		plotABC_3D(
			particleDataFrame = resultsBMExample[[1]]$particleDataFrame,
			parameter = 7,
			show.particles = "none",
			plot.parent = FALSE,
			realParam = FALSE,
			realParamValues = NA
			)	
			
		}
	})

test_that("plotPosteriors works", {

    set.seed(1)
    data(simRunExample)
    plotPosteriors(
        particleDataFrame = resultsBMExample[[1]]$particleDataFrame,
        priorsList = resultsBMExample[[1]]$priorList,
        realParam = TRUE,
        realParamValues = c(ancStateExample, genRateExample)
        )		
		
})

test_that("plotPosteriors correctly gives an error", {
    set.seed(1)
    data(simRunExample)
    expect_error(
        plotPosteriors(
            resultsBMExample[[1]]$particleDataFrame,
            resultsBMExample[[1]]$priorList)
        ,NA)
})

test_that("plotPriorPost works", {

  set.seed(1)

  data(simRunExample)
  
  plotPrior(
	priorFn = "exponential", 
	priorVariables = c(10)
	)

  plotPrior(
	priorFn = "normal", 
	priorVariables = c(1,2)
	)

  plotPrior(
	priorFn = "gamma", 
	priorVariables = c(2,0.2), 
	plotQuants = FALSE, 
	plotLegend = FALSE
	)
	
  priorKernal <- getUnivariatePriorCurve(
	priorFn = "normal",
    priorVariables = c(28, 2), 
	nPoints = 1e+05, 
	from = NULL,
    to = NULL, 
	alpha = 0.95
	)
	
  postKernal <- getUnivariatePosteriorCurve(
	acceptedValues = resultsBMExample[[1]]$particleDataFrame$starting_1,
    from = NULL, 
	to = NULL, 
	alpha = 0.95
	)
	
  priorKernal
  postKernal
  
  plotUnivariatePosteriorVsPrior(
	posteriorCurve = postKernal,
    priorCurve = priorKernal, 
	label = "parameter",
    trueValue = NULL
	)
})

test_that("parentOffspringPlots works", {
    set.seed(1)
  data(simRunExample)
  parentOffspringPlots(resultsBMExample[[1]]$particleDataFrame)
})