Raghavendran Shankar
1. The hepatic injury data set was described in the introductory chapter and
contains 281 unique compounds, each of which has been classified as causing no
liver damage, mild damage, or severe damage (Fig. 1.2). These compounds were
analyzed with 184 biological screens (i.e., experiments) to assess each
compound’s effect on a particular biologically relevant target in the body. The
larger the value of each of these predictors, the higher the activity of the
compound. In addition to biological screens, 192 chemical fingerprint predictors
were determined for these compounds. Each of these predictors represent a
substructure (i.e., an atom or combination of atoms within the compound) and
are either counts of the number of substructures or an indicator of presence or
absence of the particular substructure. The objective of this data set is to build a
predictive model for hepatic injury so that other compounds can be screened for
the likelihood of causing hepatic injury. Start R and use these commands to load
the data:
(a) Given the classification imbalance in hepatic injury status, describe how you
would create a training and testing set.
A: We use stratified random sampling to split the data to cope up with the imbalance in
hepatic injury status. Stratified random sampling is used to split the training and test data
in balance according to the hepatic status label (None, Mild, Severe) using
CreateDataPartition() method.
(b) Which classification statistic would you choose to optimize for this exercise and
why?
A: Accuracy is used as a classification statistic. Accuracy can be used to optimize as it
make good decisions to select optimal model for training and testing set. Accuracy tells
how good a classification model is functioning.
(c) Split the data into a training and a testing set, pre-process the data, and build models
described in this chapter for the biological predictors and separately for the chemical
fingerprint predictors. Which model has the best predictive ability for the biological
predictors and what is the optimal performance? Which model has the best predictive
ability for the chemical predictors and what is the optimal performance? Based on
these results, which set of predictors contains the most information about hepatic
toxicity?
A:
Biological Data:
GLM:
,PLSDA:
,
, LDA:
GLMNET: