context("test_2_impute_mat: impute_mat is ordered as original data")
X <- matrix(c(5.5,10.6,3.1,8.7,5.1,10.2,9.8,4.4,4.9), ncol=1)
Y <- matrix(c(102,105,120,130,100,80,94,108,96), ncol=1)
W <- matrix(c(1, 1, 1, 3, 2, 3, 2, 1, 2), ncol=1)
t2_in_imputemat <- multiMatch(Y,W,X,trimming = 0,match_on = "multinom")
tests_data <- quickLookup("test_toy_output.Rdata")
load(tests_data)
my_tolerance <- 1e-3
test_that("match on GPS with one X and no trimming returns same output", {
## test the estimates
expect_equal(
(t2_in_imputemat$results)$Estimate,
(baseline_tests2$results)$Estimate,
tolerance = my_tolerance
)
## test the AI 2006 variance estimates
expect_equal(
(t2_in_imputemat$results)$Variance,
(baseline_tests2$results)$Variance,
tolerance = my_tolerance
)
## test the AI 2016 variance estimates
expect_equal(
(t2_in_imputemat$results)$VarianceAI2016,
(baseline_tests2$results)$VarianceAI2012,
tolerance = my_tolerance
)
## test the impute matrix
expect_equal(
(t2_in_imputemat$impute_mat),
(baseline_tests2$impute_mat),
tolerance = my_tolerance,
check.attributes = FALSE
)
})
test_that(
"impute_mat behaves well",
{
baseline_imputes <- matrix(NA,ncol=3, nrow= length(W))
for (ii in 1:nrow(baseline_imputes)){
baseline_imputes[ii,W[ii]] <- Y[ii]
}
new_imputes <- (t2_in_imputemat$impute_mat)
new_imputes[is.na(baseline_imputes)] <- NA
expect_equal(
new_imputes[1:3,],
baseline_imputes[1:3,],
tolerance = my_tolerance,
check.attributes = FALSE
)
expect_equal(
new_imputes ,
baseline_imputes ,
tolerance = my_tolerance,
check.attributes = FALSE
)
}
)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.