Package modeval is designed to assist novice to intermediate analysts in choosing an optimal classification model, particularly for working with relatively small data sets. It provides cross-validated results comparing several different models at once using a consistent set of performance metrics, so users can hone in on the most promising approach rather than attempting single model fittings at a time. The package pre-defined 12 most common classification models, although users are free to select from the 200+ other options available in caret package. In total, these models have been proven to address a wide range of classification challenges. Pre-defined models are as follows:
There are many different modeling options in R and each has particular syntax for model training and/or prediction. At the same time, model fit structures are also different, so it is challenging to consolidate results and efficiently compare performance between models. The caret package provides a uniform interface for more than 200 regression and classification models; modeval was largely built on that platform, and is intended to be an easier, simpler entry into solving classification problems.
Presented with a classification problem, one of the first questions facing an analyst is, what kinds of models should I try? With hundreds of options available, it is daunting to work through enough models and confidently pick the one that best suits a particular dataset and classification objective. So, analysts tend to rely on models that are familiar and easy to use, which may not necessarily yield the best possible result. modeval automates model evaluation and enables comparison of performance with minimal user manipulation.
Analysts may also enter into classification model fitting with assumptions of data normality and class balance, which can each threaten their ability to achieve a good result. This package includes automated checks and recommendations for transforming and subsampling data to address issues of non-normality and class imbalance, respectively.
Following are some of the key features of modeval:
Below, several brief examples are included to demonstrate the workflow and general utility of modeval.
Here we draw a sample from a data set called "PimaIndiansDiabetes," which is embedded in R package mlbench and can also be downloaded from the University of California-Irvine's Machine Learning Repository. Included are 8 predictor variables referencing various health characteristics and one outcome class variable indicating either a positive or negative test for diabetes.
library(modeval) library(mlbench) data(PimaIndiansDiabetes) index <- sample(seq_len(nrow(PimaIndiansDiabetes)), 500) trainingSet <- PimaIndiansDiabetes[index, ] testSet <- PimaIndiansDiabetes[-index, ] x <- trainingSet[, -9] y <- trainingSet[, 9] x_test <- testSet[, -9] y_test <- testSet[, 9]
One of the first things an analyst will want to do is inspect the data set to understand the characteristics of all the variables, including the distributions of the predictor variables and the prevalence of each class outcome.
Two common challenges with classifcation problems are: (1) significant non-normality in one or more variables, and (2) class imbalance. modeval not only automates checking for both issues, but also guides the user to select a remedy, if appropriate.
In the example below, we use the
suggest_tranformation function to check for non-normality (skew and kurtosis), and view the results of applying three data transformations. (Options for addressing class imbalance are presented within function
add_model, specifically the sampling parameter.)
Four plots are generated, each displaying the skewness (x-axis) and kurtosis (y-axis) of the predictor variables. The first plot presents the untransformed data, while the others demonstrate the impact on skew and kurtosis of applying three well-established tranformation approaches. Each is labled with a tag that can be applied in the
add_model function described below, should the user wish to use one of the tranformations. The four plots:
Skew and kurtosis values of >2 and <-2 are generally considered problematic, and each plot helps the user visually inspect the extent to which variables are outside those values, and thus transformation may be useful. First, a shaded square between -2 and +2 is presented. And second, variable points are colored red, blue, or purple if skew, kurtosis, or both are outside the boundaries, respectively.
In the above example, we see that all three transformation options reduce skew and kurtosis but only tf2 (Yeo-Johnson) gets all variables inside the desired boundaries.
add_model, users can conduct model training and add each model fit to the summary list. The function will check that the model and data supports classification. If model supports class probability, then the best model is choosen based on AUC. If model doesnt support class probabilities, the best model is choosen based on accuracy. Ten-fold cross validation is used by default.
In the example below, we specify four models. Although modeval supports 12 models, the user can use additional models if they choose, as long as they are supported in caret. Here, one of those non-native models,
qrnn, is excluded from training because it only applicable to regression models, and the other,
rpartScore will only produce Accuracy and Kappa metrics as it doesn't support class probability prediction.
sSummary is a list of model fits and each model fit is named. The code below allows us to see the current contents of
sSummary as defined above.
sSummary <- list() # empty list where we store model fit results models = c("glm", "qrnn", "lda", "rpartScore") # create a character vector with function names sSummary <- add_model(sSummary, x, y, models)
## ## Model: rpartScore ## >> Accuracy and Kappa metrics are available. ## >> Best model is selected based on accuracy. ## ## ## Model: glm, lda ## >> ROC, Sens, Spec, Accuracy and Kappa metrics are available. ## >> Best model is selected based on AUC. ## ## ## Model: qrnn ## >> Model(s) not support classification problem. ##
suggest_transformation). Default = NULL.
Note that SMOTE and ROSE require the DMwR and ROSE packages, respectively.
Adding more models is easy. Simply run the
add_model function. Here are a few examples:
sSummary <- add_model(sSummary, x, y, c("knn", "nnet", "qda"), modelTag = "Nonlinear", tf="tf2") sSummary <- add_model(sSummary, x, y, c("rf", "rpart", "treebag"), modelTag = "TreeBased", sampling = "down") sSummary <- add_model(sSummary, x, y, c("svmLinear", "svmRadial", "svmPoly"), modelTag = "svmFamily")
In this case, we notice that each model fit is labeled "modelname_modelTag." The modelTag will be useful when we want to visualize results by the four main categories provided in this package (linear, nonlinear, tree-based, support vector machines).
Users can easily visualize performance with a single line of code and have the option of viewing average training time for a single tuning. Options with and without time are presented below.
suggest_auc(sSummary) suggest_auc(sSummary, time = TRUE)
If a user wishes to select specific models to compare, they can leverage modelTag. Once modelTag is added, it only includes the model fits that includes modelTag in the name.
suggest_auc(sSummary, time = TRUE, "TreeBased")
Users can also indicate model names to grep the models for comparison. The following code grabs any model fits that include "glm" or "Tree."
suggest_auc(sSummary, time = TRUE, "glm|Tree")
Users interested in visualizing accuracy and Kappa metrics can use the function
suggest_accuracy, which uses the same syntax and arguments as
suggest_auc above. As with
suggest_auc, time=TRUE will generate a training time chart. Here's an example:
With this function, users can get better understanding of caracteristics of each prediction model family. Users may also want to use modelTag to limit visualization to certain models. Example:
Some classification and regression models produce variable and importance indices based on their own algorithms. Weights and priorities are not always comparable across different model types, so average or median value from differnet models is not always useful and may be misleading. However, observing and comparing variable importance results can be very useful because we can check which variables are consistantly important across different models and which aren't. This helps the analyst develop a deeper level of understanding of the caracteristic of each model and the problem itself. Users may wish to use tag to extract the models for comparison.
suggest_variable(sSummary) suggest_variable(sSummary, "TreeBased")
To be able to use functions
suggest_gain, users will need to add prediction functions with the test dataset suing the
add_prob function. Using our sample data set, class probability predictions of a positive diabetes test ("pos") are generated using x_test and y_test.
sSummary <- add_prob(sSummary, x_test, y_test, "pos")
While AUC demonstrates overall performance, sometimes analysts will be more interested in a specific area. By observing density distribution by population and by probability cut-off, we can more clearly understand each model's predictive performance as well as similarities and differences between models.
suggest_probPop(sSummary, "pos") suggest_probPop(sSummary, "pos", modelTag = "Tree") suggest_probCut(sSummary, "pos") suggest_probCut(sSummary, "pos", modelTag = "Tree")
Gain and Lift charts are widely used for marketing purposes. They indicate the effectiveness of predictive models compared to the results obtained without the predictive model. With modeval, users can create Gain, Lift, Accumulated Event Percent, and Event Percent for each population bucket.
suggest_gain(sSummary, "pos", modelTag = "Tree") suggest_gain(sSummary, "pos", type = "Gain") suggest_gain(sSummary, "pos", type = "Lift") suggest_gain(sSummary, "pos", type = "PctAcc") suggest_gain(sSummary, "pos", type = "Pct") suggest_gain(sSummary, "pos", type = "Gain", modelTag = "Tree")
When single plots are created using the "type" argument, users can adjust the plot by adding ggplot2 syntax. In this example, we adjust the x-axis limit.
suggest_gain(sSummary, "pos", type = "Gain") + xlim(0, 0.5)
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.