@@ -3,8 +3,8 @@ Type: Package

 Title: A framework for cross-validated classification problems, with

        applications to differential variability and differential

        distribution testing

-Version: 3.1.24

-Date: 2022-10-17

+Version: 3.1.25

+Date: 2022-10-18

 Author: Dario Strbenac, Ellis Patrick, John Ormerod, Graham Mann, Jean Yang

 Maintainer: Dario Strbenac <dario.strbenac@sydney.edu.au>

 VignetteBuilder: knitr

@@ -10,8 +10,8 @@

 #' @param outcome A vector of class labels of class \code{\link{factor}} of the

 #' same length as the number of samples in \code{measurements} or a character vector of length 1 containing the

 #' column name in \code{measurements} if it is a \code{\link{DataFrame}}. Or a \code{\link{Surv}} object or a character vector of

-#' length 2 or 3 specifying the time and event columns in \code{measurements} for survival outcome.

-#' @param outcomeColumns If \code{measurements} is a \code{\link{MultiAssayExperiment}}, the column name(s) in \code{colData(measurements)} representing the outcome.

+#' length 2 or 3 specifying the time and event columns in \code{measurements} for survival outcome. If \code{measurements} is a

+#' \code{\link{MultiAssayExperiment}}, the column name(s) in \code{colData(measurements)} representing the outcome.

 #' @param outcomeTrain For the \code{train} function, either a factor vector of classes, a \code{\link{Surv}} object, or

 #' a character string, or vector of such strings, containing column name(s) of column(s)

 #' containing either classes or time and event information about survival.

@@ -201,7 +201,6 @@ Using an ordinary GLM instead.")

                   # The below loops over different combinations of assays and merges them together.

                   # This allows someone to answer which combinations of the assays might be most useful.

-

                   if(!is.list(assayCombinations) && assayCombinations == "all") assayCombinations <- do.call("c", sapply(seq_along(assayIDs), function(nChoose) combn(assayIDs, nChoose, simplify = FALSE)))

                   result <- sapply(assayCombinations, function(assayIndex){

@@ -237,7 +236,6 @@ Using an ordinary GLM instead.")

                       if(length(assayCombinations) == 0) stop("No assayCombinations with \"clinical\" data")

                   }

-

                   result <- sapply(assayCombinations, function(assayIndex){

                       CV(measurements = measurements[, mcols(measurements)[["assay"]] %in% assayIndex],

                          outcome = outcome, assayIDs = assayIndex,

@@ -299,7 +297,7 @@ Using an ordinary GLM instead.")

 # One or more omics data sets, possibly with clinical data.

 setMethod("crossValidate", "MultiAssayExperiment",

           function(measurements,

-                   outcomeColumns, 

+                   outcome, 

                    nFeatures = 20,

                    selectionMethod = "t-test",

                    selectionOptimisation = "Resubstitution",

@@ -312,7 +310,7 @@ setMethod("crossValidate", "MultiAssayExperiment",

                    nCores = 1,

                    characteristicsLabel = NULL, ...)

           {

-              measurementsAndOutcome <- prepareData(measurements, outcomeColumns, ...)

+              measurementsAndOutcome <- prepareData(measurements, outcome, ...)

               crossValidate(measurements = measurementsAndOutcome[["measurements"]],

                             outcome = measurementsAndOutcome[["outcome"]], 

@@ -691,6 +689,7 @@ generateMultiviewParams <- function(assayIDs,

                                           nFeatures = nFeatures,

                                           selectionMethod = selectionMethod,

                                           selectionOptimisation = "none",

+                                          performanceType = performanceType,

                                           classifier = classifier,

                                           multiViewMethod = "none")

@@ -1043,9 +1042,9 @@ train.list <- function(x, outcomeTrain, ...)

 #' @rdname crossValidate

 #' @method train MultiAssayExperiment

 #' @export

-train.MultiAssayExperiment <- function(x, outcomeColumns, ...)

+train.MultiAssayExperiment <- function(x, outcome, ...)

           {

-              prepArgs <- list(x, outcomeColumns)

+              prepArgs <- list(x, outcome)

               extraInputs <- list(...)

               prepExtras <- trainExtras <- numeric()

               if(length(extraInputs) > 0)

@@ -165,11 +165,13 @@ setMethod("prepareData", "MultiAssayExperiment",

   if(!all(useFeatures[, "assay"] %in% c(names(measurements), "clinical")))

     stop("Some assay names in first column of 'useFeatures' are not assay names in 'measurements' or \"clinical\".")

-  clinicalColumns <- colnames(MultiAssayExperiment::colData(measurements))

+  clinicalColumnsDataset <- colnames(MultiAssayExperiment::colData(measurements))

   if("clinical" %in% useFeatures[, "assay"])

   {

     clinicalRows <- useFeatures[, "assay"] == "clinical"      

     clinicalColumns <- useFeatures[clinicalRows, "feature"]

+    if(length(clinicalColumns) == 1 && clinicalColumns == "all")

+      clinicalColumns <- setdiff(clinicalColumnsDataset, outcomeColumns)

     useFeatures <- useFeatures[!clinicalRows, ]

   } else {

     clinicalColumns <- NULL

@@ -178,7 +180,6 @@ setMethod("prepareData", "MultiAssayExperiment",

   if(nrow(useFeatures) > 0)

   {

     measurements <- measurements[, , unique(useFeatures[, "assay"])]

-  

     # Get all desired measurements tables and clinical columns (other than the columns representing outcome).

     # These form the independent variables to be used for making predictions with.

     # Variable names will have names like RNA_BRAF for traceability.

@@ -13,7 +13,10 @@

 #' are features.

 #' @param outcome Either a factor vector of classes, a \code{\link{Surv}} object, or

 #' a character string, or vector of such strings, containing column name(s) of column(s)

-#' containing either classes or time and event information about survival.

+#' containing either classes or time and event information about survival. If

+#' \code{measurements} is a \code{MultiAssayExperiment}, the names of the column (class) or

+#' columns (survival) in the table extracted by \code{colData(data)} that contain(s) the samples'

+#' outcome to use for prediction.

 #' @param crossValParams An object of class \code{\link{CrossValParams}},

 #' specifying the kind of cross-validation to be done.

 #' @param modellingParams An object of class \code{\link{ModellingParams}},

@@ -26,9 +29,6 @@

 #' package.  Transformation, selection and prediction functions provided by

 #' this package will cause the characteristics to be automatically determined

 #' and this can be left blank.

-#' @param outcomeColumns If \code{measurementsTrain} is a \code{MultiAssayExperiment}, the

-#' names of the column (class) or columns (survival) in the table extracted by \code{colData(data)}

-#' that contain(s)s the samples' outcome to use for prediction.

 #' @param ... Variables not used by the \code{matrix} nor the \code{MultiAssayExperiment} method which

 #' are passed into and used by the \code{DataFrame} method or passed onwards to \code{\link{prepareData}}.

 #' @param verbose Default: 1. A number between 0 and 3 for the amount of

@@ -70,9 +70,11 @@ setMethod("runTests", c("matrix"), function(measurements, outcome, ...) # Matrix

 setMethod("runTests", "DataFrame", function(measurements, outcome, crossValParams = CrossValParams(), modellingParams = ModellingParams(),

            characteristics = S4Vectors::DataFrame(), ..., verbose = 1)

 {

-  # Get out the outcome if inside of data table.           

   if(is.null(rownames(measurements)))

-    stop("'measurements' DataFrame must have sample identifiers as its row names.")

+  {

+    warning("'measurements' DataFrame must have sample identifiers as its row names. Generating generic ones.")

+    rownames(measurements) <- paste("Sample", seq_len(nrow(measurements)))

+  }

   if(any(is.na(measurements)))

     stop("Some data elements are missing and classifiers don't work with missing data. Consider imputation or filtering.")            

@@ -93,7 +95,7 @@ input data. Autmomatically reducing to smaller number.")

   # Element names of the list returned by runTest, in order.

   resultTypes <- c("ranked", "selected", "models", "testSet", "predictions", "tune", "importance")

-  

+

   # Create all partitions of training and testing sets.

   samplesSplits <- .samplesSplits(crossValParams, outcome)

   splitsTestInfo <- .splitsTestInfo(crossValParams, samplesSplits)

@@ -112,6 +114,7 @@ input data. Autmomatically reducing to smaller number.")

       message("Processing sample set ", setNumber, '.')

     # crossValParams is needed at least for nested feature tuning.

+    

     runTest(measurements[trainingSamples, , drop = FALSE], outcome[trainingSamples],

             measurements[testSamples, , drop = FALSE], outcome[testSamples],

             crossValParams, modellingParams, characteristics, verbose,

@@ -180,9 +183,9 @@ input data. Autmomatically reducing to smaller number.")

 #' @import MultiAssayExperiment methods

 #' @export

 setMethod("runTests", c("MultiAssayExperiment"),

-          function(measurements, outcomeColumns, ...)

+          function(measurements, outcome, ...)

 {

-  prepArgs <- list(measurements, outcomeColumns)              

+  prepArgs <- list(measurements, outcome)              

   extraInputs <- list(...)

   prepExtras <- numeric()

   if(length(extraInputs) > 0)

@@ -37,7 +37,7 @@ crossValidate(measurements, outcome, ...)

 \S4method{crossValidate}{MultiAssayExperiment}(

   measurements,

-  outcomeColumns,

+  outcome,

   nFeatures = 20,

   selectionMethod = "t-test",

   selectionOptimisation = "Resubstitution",

@@ -119,7 +119,7 @@ crossValidate(measurements, outcome, ...)

 \method{train}{list}(x, outcomeTrain, ...)

-\method{train}{MultiAssayExperiment}(x, outcomeColumns, ...)

+\method{train}{MultiAssayExperiment}(x, outcome, ...)

 \method{predict}{trainedByClassifyR}(object, newData, ...)

 }

@@ -130,7 +130,8 @@ or a list of these objects containing the data.}

 \item{outcome}{A vector of class labels of class \code{\link{factor}} of the

 same length as the number of samples in \code{measurements} or a character vector of length 1 containing the

 column name in \code{measurements} if it is a \code{\link{DataFrame}}. Or a \code{\link{Surv}} object or a character vector of

-length 2 or 3 specifying the time and event columns in \code{measurements} for survival outcome.}

+length 2 or 3 specifying the time and event columns in \code{measurements} for survival outcome. If \code{measurements} is a

+\code{\link{MultiAssayExperiment}}, the column name(s) in \code{colData(measurements)} representing the outcome.}

 \item{...}{Parameters passed into \code{\link{prepareData}} which control subsetting and filtering of input data.}

@@ -162,8 +163,6 @@ with each element being a vector of assays to combine. Special value \code{"all"

 \item{characteristicsLabel}{A character specifying an additional label for the cross-validation run.}

-\item{outcomeColumns}{If \code{measurements} is a \code{\link{MultiAssayExperiment}}, the column name(s) in \code{colData(measurements)} representing the outcome.}

-

 \item{x}{Same as \code{measurements} but only training samples.}

 \item{outcomeTrain}{For the \code{train} function, either a factor vector of classes, a \code{\link{Surv}} object, or

@@ -19,7 +19,7 @@

   verbose = 1

 )

-\S4method{runTests}{MultiAssayExperiment}(measurements, outcomeColumns, ...)

+\S4method{runTests}{MultiAssayExperiment}(measurements, outcome, ...)

 }

 \arguments{

 \item{measurements}{Either a \code{\link{matrix}}, \code{\link{DataFrame}}

@@ -32,7 +32,10 @@ are passed into and used by the \code{DataFrame} method or passed onwards to \co

 \item{outcome}{Either a factor vector of classes, a \code{\link{Surv}} object, or

 a character string, or vector of such strings, containing column name(s) of column(s)

-containing either classes or time and event information about survival.}

+containing either classes or time and event information about survival. If

+\code{measurements} is a \code{MultiAssayExperiment}, the names of the column (class) or

+columns (survival) in the table extracted by \code{colData(data)} that contain(s) the samples'

+outcome to use for prediction.}

 \item{crossValParams}{An object of class \code{\link{CrossValParams}},

 specifying the kind of cross-validation to be done.}

@@ -52,10 +55,6 @@ and this can be left blank.}

 \item{verbose}{Default: 1. A number between 0 and 3 for the amount of

 progress messages to give.  A higher number will produce more messages as

 more lower-level functions print messages.}

-

-\item{outcomeColumns}{If \code{measurementsTrain} is a \code{MultiAssayExperiment}, the

-names of the column (class) or columns (survival) in the table extracted by \code{colData(data)}

-that contain(s)s the samples' outcome to use for prediction.}

 }

 \value{

 An object of class \code{\link{ClassifyResult}}.