Nothing
library(smcfcs)
library(survival)
context("Error trap testing")
test_that("Checking outcome model check for logistic models", {
expect_error({
set.seed(1234)
n <- 100
x <- rnorm(n)
y <- as.factor(1 * (runif(n) < 0.5))
x[(runif(n) < 0.5)] <- NA
simData <- data.frame(x, y)
imps <- smcfcs(simData,
smtype = "logistic", smformula = "y~x",
method = c("norm", "")
)
})
})
test_that("Checking outcome model check for logistic models", {
expect_error({
set.seed(1234)
n <- 100
x <- rnorm(n)
y <- 1 + 1 * (runif(n) < 0.5)
x[(runif(n) < 0.5)] <- NA
simData <- data.frame(x, y)
imps <- smcfcs(simData,
smtype = "logistic", smformula = "y~x",
method = c("norm", "")
)
})
})
test_that("Checking outcome model check for logistic models", {
expect_output({
set.seed(1234)
n <- 100
x <- rnorm(n)
y <- 1 * (runif(n) < 0.5)
x[(runif(n) < 0.5)] <- NA
simData <- data.frame(x, y)
imps <- smcfcs(simData,
smtype = "logistic", smformula = "y~x",
method = c("norm", "")
)
})
})
test_that("Checking error trap for method statement 1", {
expect_error({
set.seed(1234)
n <- 100
x1 <- rnorm(n)
x2 <- rnorm(n)
y <- y <- x1 + x2 + rnorm(n)
x1[(runif(n) < 0.5)] <- NA
simData <- data.frame(x1, x2, y)
imps <- smcfcs(simData,
smtype = "lm", smformula = "y~x1+x2",
method = c("", "", "")
)
})
})
test_that("Checking error trap for method statement 2", {
expect_error({
set.seed(1234)
n <- 100
x1 <- rnorm(n)
x2 <- rnorm(n)
y <- y <- x1 + x2 + rnorm(n)
x1[(runif(n) < 0.5)] <- NA
simData <- data.frame(x1, x2, y)
imps <- smcfcs(simData,
smtype = "lm", smformula = "y~x1+x2",
method = c("", "norm", "")
)
})
})
test_that("Checking measurement error error checks 1", {
expect_error({
set.seed(1234)
n <- 100
x <- rnorm(n)
w1 <- x + rnorm(n)
w2 <- x + rnorm(n)
y <- y <- x + rnorm(n)
x <- rep(NA, n)
simData <- data.frame(x, w1, w2, y)
imps <- smcfcs(simData,
smtype = "lm", smformula = "y~x",
method = c("latnorm", "", "", "")
)
})
})
test_that("Checking measurement error error checks 2", {
expect_error({
set.seed(1234)
n <- 100
x <- rnorm(n)
w1 <- x + rnorm(n)
y <- y <- x + rnorm(n)
x <- rep(NA, n)
simData <- data.frame(x, w1, y)
errMat <- array(0, dim = c(3, 3))
imps <- smcfcs(simData,
smtype = "lm", smformula = "y~x",
method = c("latnorm", "", ""), errorProneMatrix = errMat
)
})
})
test_that("Checking measurement error error checks 3", {
expect_error({
set.seed(1234)
n <- 100
x <- rnorm(n)
w1 <- x + rnorm(n)
y <- y <- x + rnorm(n)
x <- rep(NA, n)
simData <- data.frame(x, w1, y)
errMat <- array(0, dim = c(3, 3))
errMat[1, 1] <- 1
imps <- smcfcs(simData,
smtype = "lm", smformula = "y~x",
method = c("latnorm", "", ""), errorProneMatrix = errMat
)
})
})
test_that("Checking measurement error error checks 4", {
expect_error({
set.seed(1234)
n <- 100
x <- rnorm(n)
w1 <- x + rnorm(n)
w2 <- x + rnorm(n)
y <- y <- x + rnorm(n)
x <- rep(NA, n)
simData <- data.frame(x, w1, w2, y)
errMat <- array(0, dim = c(4, 4))
errMat[1, 2] <- 1
errMat[1, 3] <- 1
errMat[4, 2] <- 1
imps <- smcfcs(simData,
smtype = "lm", smformula = "y~x",
method = c("latnorm", "", "", ""), errorProneMatrix = errMat
)
})
})
test_that("Checking measurement error error checks 5", {
expect_error({
set.seed(1234)
n <- 100
x <- rnorm(n)
w1 <- x + rnorm(n)
w2 <- x + rnorm(n)
y <- y <- x + rnorm(n)
x <- rep(NA, n)
simData <- data.frame(x, w1, w2, y)
errMat <- array(0, dim = c(4, 4))
errMat[1, 2] <- 1
errMat[1, 3] <- 1
errMat[4, 2] <- 2
imps <- smcfcs(simData,
smtype = "lm", smformula = "y~x",
method = c("latnorm", "", "", ""), errorProneMatrix = errMat
)
})
})
test_that("Checking measurement error error checks 6", {
expect_error({
set.seed(1234)
n <- 100
x <- rnorm(n)
w1 <- x + rnorm(n)
w2 <- x + rnorm(n)
y <- y <- x + rnorm(n)
x <- rep(NA, n)
simData <- data.frame(x, w1, w2, y)
errMat <- array(0, dim = c(5, 5))
errMat[1, 2] <- 1
errMat[1, 3] <- 1
imps <- smcfcs(simData,
smtype = "lm", smformula = "y~x",
method = c("latnorm", "", "", ""), errorProneMatrix = errMat
)
})
})
test_that("Checking measurement error error checks 7", {
expect_error({
set.seed(1234)
n <- 100
x <- rnorm(n)
w1 <- x + rnorm(n)
w2 <- x + rnorm(n)
y <- y <- x + rnorm(n)
x <- rep(NA, n)
z <- rnorm(n)
z[1:50] <- NA
simData <- data.frame(w1, w2, y, z)
errMat <- array(0, dim = c(4, 4))
errMat[3, 1] <- 1
errMat[3, 2] <- 1
imps <- smcfcs(simData,
smtype = "lm", smformula = "y~x",
method = c("", "", "", "norm"), errorProneMatrix = errMat
)
})
})
test_that("Cox imputation fails if event indicator is not coded right", {
expect_error({
set.seed(1234)
n <- 1000
z <- rnorm(n)
x <- z + rnorm(n)
t <- -log(runif(n)) / (1 * exp(x + z))
d <- 1 * (t < 10) + 1
t[d == 1] <- 10
x[(runif(n) < 0.5)] <- NA
simData <- data.frame(t, d, x, z)
imps <- smcfcs(simData,
smtype = "coxph", smformula = "Surv(t, d)~x+z",
method = c("", "", "norm", "")
)
})
})
test_that("Cox error check doesn't fail when events are 1 1 1 0 0 0", {
expect_error(
{
set.seed(1234)
n <- 100
z <- rnorm(n)
x <- z + rnorm(n)
t <- runif(n)
d <- c(rep(1, n / 2), rep(0, n / 2))
x[(runif(n) < 0.5)] <- NA
simData <- data.frame(t, d, x, z)
imps <- smcfcs(simData,
smtype = "coxph", smformula = "Surv(t, d)~x+z",
method = c("", "", "norm", "")
)
},
NA
)
})
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.