Create an array of ROC plots to compare signatures.
Usage
signatureROCplot(
inputData,
annotationData,
signatureColNames,
annotationColName,
scale = FALSE,
choose_colors = c("cornflowerblue", "gray24"),
name = "Signatures",
nrow = NULL,
ncol = NULL
)
Arguments
- inputData
an input data object. It should either be of the class
SummarizedExperiment
and contain the profiled signature data and annotation data as columns in thecolData
, or alternatively be of the classesdata.frame
ormatrix
and contain only the gene expression data. Required.- annotationData
a
data.frame
ormatrix
of annotation data, with one column. Only required ifinputData
is adata.frame
ormatrix
of signature data.- signatureColNames
a
vector
of the column names ofinputData
that contain the signature data. IfinputData
is aSummarizedExperiment
object, these are the column names of the objectcolData
.- annotationColName
a character string naming the column name in the
colData
that contains the annotation data to be used in making the boxplot. Only required if inputData is aSummarizedExperiment
object.- scale
logical. Setting
scale = TRUE
scales the signature data. The default isFALSE
.- choose_colors
a
vector
of length 2 defining the colors to be used in the ROC plots. The default isc("cornflowerblue", "gray24")
.- name
a character string giving the title of the boxplot. The default is
"Signatures"
.- nrow
integer giving the number of rows in the resulting array.
- ncol
integer giving the number of columns in the resulting array.
Examples
# Run signature profiling
choose_sigs <- subset(TBsignatures,
!(names(TBsignatures) %in% c("Lee_4", "Roe_OD_4")))[c(1,2)]
prof_indian <- runTBsigProfiler(TB_indian, useAssay = "logcounts",
algorithm = "ssGSEA",
signatures = choose_sigs,
parallel.sz = 1)
#> Parameter update_genes is TRUE. Gene names will be updated.
#> Running ssGSEA
#> Warning: 1 genes with constant expression values throuhgout the samples.
#> Estimating ssGSEA scores for 2 gene sets.
#> [1] "Calculating ranks..."
#> [1] "Calculating absolute values from ranks..."
#>
|
| | 0%
|
|== | 2%
|
|=== | 5%
|
|===== | 7%
|
|====== | 9%
|
|======== | 11%
|
|========== | 14%
|
|=========== | 16%
|
|============= | 18%
|
|============== | 20%
|
|================ | 23%
|
|================== | 25%
|
|=================== | 27%
|
|===================== | 30%
|
|====================== | 32%
|
|======================== | 34%
|
|========================= | 36%
|
|=========================== | 39%
|
|============================= | 41%
|
|============================== | 43%
|
|================================ | 45%
|
|================================= | 48%
|
|=================================== | 50%
|
|===================================== | 52%
|
|====================================== | 55%
|
|======================================== | 57%
|
|========================================= | 59%
|
|=========================================== | 61%
|
|============================================= | 64%
|
|============================================== | 66%
|
|================================================ | 68%
|
|================================================= | 70%
|
|=================================================== | 73%
|
|==================================================== | 75%
|
|====================================================== | 77%
|
|======================================================== | 80%
|
|========================================================= | 82%
|
|=========================================================== | 84%
|
|============================================================ | 86%
|
|============================================================== | 89%
|
|================================================================ | 91%
|
|================================================================= | 93%
|
|=================================================================== | 95%
|
|==================================================================== | 98%
|
|======================================================================| 100%
#>
#> [1] "Normalizing..."
# Create ROC plots
signatureROCplot(prof_indian, signatureColNames = names(choose_sigs),
annotationColName = "label")