context("GoodnessOfFitPlot")
zformula <- formula("Overall ~ Fees + Interest + Phone + Branch + Online + ATM")
data(bank, package = "flipExampleData")
bank <- bank[sample(nrow(bank), 300), ] # random sample of 300 rows to improve perfomance
sb <- bank$ID > 100
attr(sb, "label") <- "ID greater than 100"
wgt <- bank$ID
wgt[is.na(wgt)] = 0
attr(wgt, "label") <- "ID"
attr(bank$Overall, "label") <- "Overall satisfaction"
attr(bank$Fees, "label") <- "Fees paid"
attr(bank$Online, "label") <- "Online banking"
for(missing in c("Multiple imputation", "Imputation (replace missing values with estimates)",
"Exclude cases with missing data", "Dummy variable adjustment"))
for (type in c( "Linear","Poisson", "Quasi-Poisson", "Binary Logit", "Ordered Logit", "NBD", "Multinomial Logit"))
test_that(paste("GoodnessOfFitPlot:", missing, type),
{
z = suppressWarnings(z <- Regression(Overall ~ Fees + Interest + Phone + Branch + Online + ATM, missing = missing, data = bank, type = type))
expect_error(suppressWarnings(GoodnessOfFitPlot(z)), NA)
z = suppressWarnings(z <- Regression(Overall ~ Fees + Interest + Phone + Branch + Online + ATM, missing = missing, data = bank, subset = sb, type = type))
expect_error(suppressWarnings(GoodnessOfFitPlot(z)), NA)
z = suppressWarnings(z <- Regression(Overall ~ Fees + Interest + Phone + Branch + Online + ATM, weights = wgt, missing = missing, data = bank, type = type))
expect_error(suppressWarnings(GoodnessOfFitPlot(z)), NA)
z = suppressWarnings(z <- Regression(Overall ~ Fees + Interest + Phone + Branch + Online + ATM, weights = wgt, missing = missing, data = bank, subset = sb, type = type))
expect_error(suppressWarnings(GoodnessOfFitPlot(z)), NA)
})
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.