context("test_coefs")
library(survey)
library(svrepmisc)
library(nnet)
library(quantreg)
library(MASS)
library(crch)
# use survey package example
# replicate weights - jackknife (this is slower)
data(api)
dstrat<-svydesign(id=~1,strata=~stype, weights=~pw,
data=apistrat, fpc=~fpc)
jkstrat<-as.svrepdesign(dstrat)
test_that("coefs work for multinomial", {
svy <- as.vector(svymultinom(factor(stype)~factor(sch.wide)+api00+api99+meals,jkstrat))
reg <- coef(nnet::multinom(factor(stype)~factor(sch.wide)+api00+api99+meals,apistrat,weights=apistrat$pw))
reg <- as.vector(svrepmisc:::mat2vec(reg))
expect_equal(svy,reg)
})
test_that("coefs work for quantile", {
svy <- as.vector(svyrq(api00~api99+factor(sch.wide)+meals,jkstrat))
reg <- coef(quantreg::rq(api00~api99+factor(sch.wide)+meals,data=apistrat,weights=apistrat$pw))
reg <- as.vector(reg)
expect_equal(svy,reg)
})
test_that("coefs work for quantile smoother", {
svy <- as.vector(svyrqss(api00~api99+factor(sch.wide)+meals,jkstrat))
reg <- coef(quantreg::rqss(api00~api99+factor(sch.wide)+meals,data=apistrat,weights=apistrat$pw))
reg <- as.vector(reg)
expect_equal(svy,reg)
})
test_that("coefs work for negative binomial", {
svy <- as.vector(svynb(meals~api99+factor(sch.wide)+api00,jkstrat))
reg <- coef(MASS::glm.nb(meals~api99+factor(sch.wide)+api00,data=apistrat,weights=apistrat$pw))
reg <- as.vector(reg)
expect_equal(svy,reg)
})
test_that("coefs work for truncated reg", {
svy <- as.vector(svytrch(acs.core~api00+api99+factor(sch.wide)+meals,jkstrat,left=15))
reg <- coef(crch::trch(acs.core~api00+api99+factor(sch.wide)+meals,data=apistrat,weights=apistrat$pw,left=15))
reg <- as.vector(reg)
expect_equal(svy,reg)
})
test_that("coefs work for censored reg", {
svy <- as.vector(svycrch(acs.core~api00+api99+factor(sch.wide)+meals,jkstrat,left=15))
reg <- coef(crch::crch(acs.core~api00+api99+factor(sch.wide)+meals,data=apistrat,weights=apistrat$pw,left=15))
reg <- as.vector(reg)
expect_equal(svy,reg)
})
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.