Nothing
## test-stratified.R
## started 2021-05-03
## stratification openCR >= 2.0.0
library(openCR)
Sys.setenv(RCPP_PARALLEL_BACKEND = "tinythread")
# create small working proximity dataset
suppressWarnings(smallCH <- join(subset(ovenCHp, sessions = 1:3, traps = 1:20, dropnullocc = FALSE)))
suppressWarnings(stratifiedch3 <- stratify (smallCH, covariate = 'Sex'))
suppressWarnings(ch2 <- subset(ovenCHp, sessions = 1:3, traps = 1:20, dropnullocc = FALSE))
stratifiedch2 <- stratify (ch2[1:3])
test_that("stratify generates suitable data for openCR.fit", {
## case 3: split joined data by covariate
expect_equal(summary(stratifiedch3, terse = TRUE),
matrix(c(29,21,15,20,29,39,13,20), nrow = 4),
check.attributes = FALSE)
## case 2: split joined data by covariate
expect_equal(summary(stratifiedch2, terse = TRUE),
matrix(c(9,18,11,20,10,21,11,20,10,21,13,20), nrow = 4),
check.attributes = FALSE)
})
test_that("correct stratified likelihood", {
args <- list(capthist = stratifiedch3, type = 'PLBf',
start = list(p = 0.4, phi = 0.497, f = 0.560),
model = p ~ stratum, details = list(LLonly = TRUE),
stratified = TRUE)
expect_equal(do.call(openCR.fit, args)[1], -246.589008072,
tolerance = 1e-4, check.attributes = FALSE)
})
## 2022-01-28 Heiko Hinneberg phidot butterphi
test_that("successful predict() with parameter b and varying sessions per stratum", {
# fitting too slow, so pre-fit 2022-01-28
# suppressWarnings(ch3 <- subset(ovenCHp, sessions = 1:2, traps = 1:20, dropnullocc = FALSE))
# ch2 <- reduce(ch2, by='all') # 3 sessions
# ch3 <- reduce(ch3, by='all') # 2 sessions
# stratifiedch4 <- stratify(list(a=ch2, b=ch3))
# args <- list(capthist = stratifiedch4, type = 'PLBb',
# start = list(p = 0.4, phi = 0.5, b = 0.1),
# model = p ~ stratum, details = list(LLonly = FALSE),
# stratified = TRUE, trace=T)
# stratifiedfit <- do.call(openCR.fit, args)
# saveRDS(stratifiedfit, file = 'd:/open populations/openCR/inst/exampledata/stratifiedfit.RDS')
stratifiedfit <- readRDS(system.file("exampledata", "stratifiedfit.RDS", package = "openCR"))
expect_silent(pred <- predict(stratifiedfit))
})
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.