tests/testthat/test_null_model_logistic_regression.R

In UW-GAC/genesis2_tests: genesis assocation testing code rewrite

context("check null model logistic regression")

test_that("logistic", {
    dat <- .testNullInputs(binary=TRUE)
    
    nullmod <- fitNullMod(dat$y, dat$X, family="binomial", verbose=FALSE)

    glm.mod <- glm(dat$y ~ -1 + dat$X, family = "binomial")

    expect_equal(nullmod$family$family, "binomial")
    expect_false(nullmod$hetResid)
    expect_equivalent(nullmod$fitted.values, fitted(glm.mod))
    expect_false(nullmod$family$mixedmodel)
    expect_equivalent(nullmod$resid.marginal, resid(glm.mod, type = "response"))
    expect_true(all(nullmod$fixef == summary(glm.mod)$coef))
    expect_equivalent(nullmod$varComp, fitted(glm.mod)*(1-fitted(glm.mod)))
    expect_null(nullmod$varCompCov)
    expect_equivalent(nullmod$betaCov, vcov(glm.mod))
    expect_equivalent(nullmod$fitted.values, fitted(glm.mod))
    expect_equal(nullmod$logLik, as.numeric(logLik(glm.mod)))
    expect_equal(nullmod$AIC, AIC(glm.mod))
    expect_equivalent(nullmod$workingY, dat$y)
    expect_equivalent(nullmod$outcome, dat$y)
    expect_equivalent(nullmod$model.matrix, dat$X)
    expect_equivalent(diag(nullmod$cholSigmaInv), 1/sqrt(fitted(glm.mod)*(1-fitted(glm.mod))))
    expect_equal(nullmod$converged, glm.mod$converged)
    expect_null(nullmod$zeroFLAG)
    expect_equal(nullmod$RSS, 1)
})