SignatureQuantitative {TBSignatureProfiler} | R Documentation

## Use logistic regression and bootstrap LOOCV to evaluate signatures.

### Description

This function takes as input a `data.frame` with genetic expression count data, and uses a bootstrapped leave-one-out cross validation procedure with logistic regression to allow for numeric and graphical comparison across any number of genetic signatures.

### Usage

    SignatureQuantitative(
      df.input,
      targetVec.num,
      signature.list = NULL,
      signature.name.vec = NULL,
      num.boot = 100,
      pb.show = TRUE
    )

### Arguments

`df.input` | a `data.frame` of gene expression count data. Required.
`targetVec.num` | a numeric binary vector of the response variable. The vector should be the same number of rows as `df`. Required.
`signature.list` | a `list` of signatures to run with their associated genes. This list should be in the same format as `TBsignatures`, included in the TBSignatureProfiler package. If `signature.list = NULL`, the default set of signatures `TBsignatures` list is used. For details, run `?TBsignatures`.
`signature.name.vec` | A vector specifying the names of the signatures to be compared. This should be the same length as `signature.list`. If `signature.name.vec = NULL`, the default set of signatures `TBsignatures` list is used.
`num.boot` | an integer specifying the number of bootstrap iterations.
`pb.show` | logical. If `TRUE` then a progress bar for the bootstrapping procedure will be displayed as output. The default is `TRUE`.
`name` | a character string giving a name for the outputted boxplot of bootstrapped AUCs. The default is `"Quantitative Evaluation of Signatures via Bootstrapped AUCs"`.

### Value

the AUC, sensitivity and specificity

### Examples

    inputTest <- matrix(rnorm(1000), 100, 20,
                        dimnames = list(paste0("gene", seq.int(1, 100)),
                                        paste0("sample", seq.int(1, 20))))
    inputTest <- as.data.frame(inputTest)
    targetVec <- sample(c(0,1), replace = TRUE, size = 20)
    signature.list <- list(sig1 = c("gene1", "gene2", "gene3"),
                           sig2 = c("gene4", "gene5", "gene6"))
    signature.name.vec <- c("sig1", "sig2")
    num.boot <- 5
    SignatureQuantitative(inputTest, targetVec.num = targetVec,
                          signature.list = signature.list,
                          signature.name.vec = signature.name.vec,
                          num.boot = num.boot)

[Package TBSignatureProfiler version 1.4.11 Index]