Nothing
context("Data extraction and simulation")
library(missCompare)
data("clindata_miss")
small <- clindata_miss[1:80, 1:4]
small$age[1:15] <- NA
# expecting error if non numeric variables are in the dataframe
test_that("numeric errors in get_data()", {
expect_error(get_data(small))
})
cleaned <- clean(small)
y <- get_data(cleaned, matrixplot_sort = T)
simulated <- simulate(rownum = y$Rows, colnum =y$Columns, cormat=y$Corr_matrix)
res <- all_patterns(X_hat = simulated$Simulated_matrix,
MD_pattern = y$MD_Pattern,
NA_fraction = y$Fraction_missingness,
min_PDM = 3)
# checking no errors in functions
test_that("no errors in get_data()", {
expect_error(get_data(cleaned), NA)
})
test_that("no errors in simulate()", {
expect_error(simulate(rownum = y$Rows, colnum =y$Columns, cormat=y$Corr_matrix), NA)
})
# checking output
test_that("correlation matrix dimensions correct", {
expect_equal(dim(y$Corr_matrix), c(ncol(small),ncol(small)))
})
test_that("fraction of missingness vector length correct", {
expect_equal(length(y$Fraction_missingness_per_variable), ncol(small))
})
test_that("equal dimensions of simulated matrix and original dataframe", {
expect_equal(dim(simulated$Simulated_matrix), dim(cleaned))
})
test_that("equal dimensions of simulated matrix and MCAR matrix", {
expect_equal(dim(simulated$Simulated_matrix), dim(res$MCAR_matrix))
})
test_that("equal dimensions of simulated matrix and MNAR matrix", {
expect_equal(dim(simulated$Simulated_matrix), dim(res$MNAR_matrix))
})
test_that("equal dimensions of simulated matrix and MAR matrix", {
expect_equal(dim(simulated$Simulated_matrix), dim(res$MAR_matrix))
})
rm(list=ls())
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.