scGPS introduction

Quan Nguyen and Michael Thompson

2020-10-27

1. Installation instruction

# To install scGPS from github (Depending on the configuration of the local
# computer or HPC, possible custom C++ compilation may be required - see
# installation trouble-shootings below)
devtools::install_github("IMB-Computational-Genomics-Lab/scGPS")

# for C++ compilation trouble-shooting, manual download and installation can be
# done from github

git clone https://github.com/IMB-Computational-Genomics-Lab/scGPS

# then check in scGPS/src if any of the precompiled (e.g.  those with *.so and
# *.o) files exist and delete them before recompiling

# then with the scGPS as the R working directory, manually install and load
# using devtools functionality
# Install the package
devtools::install()
#load the package to the workspace 
library(scGPS)

2. A simple workflow of the scGPS:

The purpose of this workflow is to solve the following task:

2.1 Create scGPS objects

2.2 Run prediction

2.3 Summarise results

3. A complete workflow of the scGPS:

The purpose of this workflow is to solve the following task:

3.1 Identify clusters in a dataset using CORE

(skip this step if clusters are known)

3.2 Identify clusters in a dataset using SCORE (Stable Clustering at Optimal REsolution)

(skip this step if clusters are known)

(SCORE aims to get stable subpopulation results by introducing bagging aggregation and bootstrapping to the CORE algorithm)

3.3 Visualise all cluster results in all iterations

3.4 Compare clustering results with other dimensional reduction methods (e.g., tSNE)

3.5 Find gene markers and annotate clusters

#load gene list (this can be any lists of user-selected genes)
genes <-training_gene_sample
genes <-genes$Merged_unique

#the gene list can also be objectively identified by differential expression
#analysis cluster information is requied for find_markers. Here, we use
#CORE results.

#colData(mixedpop2)[,1] <- unlist(SCORE_test$Cluster[SCORE_test$optimal_index])

suppressMessages(library(locfit))
suppressMessages(library(DESeq))

DEgenes <- find_markers(expression_matrix=assay(mixedpop2),
                            cluster = colData(mixedpop2)[,1],
                            selected_cluster=unique(colData(mixedpop2)[,1]))

#the output contains dataframes for each cluster.
#the data frame contains all genes, sorted by p-values
names(DEgenes)
#> [1] NA                      NA                      NA                     
#> [4] NA                      "DE_Subpop1vsRemaining" "DE_Subpop2vsRemaining"
#> [7] "DE_Subpop3vsRemaining" "DE_Subpop4vsRemaining"

#you can annotate the identified clusters
DEgeneList_1vsOthers <- DEgenes$DE_Subpop1vsRemaining$id

#users need to check the format of the gene input to make sure they are
#consistent to the gene names in the expression matrix

#the following command saves the file "PathwayEnrichment.xlsx" to the
#working dir
#use 500 top DE genes
suppressMessages(library(DOSE))
suppressMessages(library(ReactomePA))
suppressMessages(library(clusterProfiler))
genes500 <- as.factor(DEgeneList_1vsOthers[seq_len(500)])
enrichment_test <- annotate_clusters(genes, pvalueCutoff=0.05, gene_symbol=TRUE)

#the enrichment outputs can be displayed by running
clusterProfiler::dotplot(enrichment_test, showCategory=10, font.size = 6)

4. Relationship between clusters within one sample or between two samples

The purpose of this workflow is to solve the following task:

4.1 Start the scGPS prediction to find relationship between clusters

4.2 Display summary results for the prediction

#get the number of rows for the summary matrix
row_cluster <-length(unique(colData(mixedpop2)[,1]))

#summary results LDA to to show the percent of cells classified as cells
#belonging by LDA classifier
summary_prediction_lda(LSOLDA_dat=LSOLDA_dat, nPredSubpop = row_cluster )
#>                 V1               V2                                names
#> 1                0 68.9839572192513 LDA for subpop 1 in target mixedpop2
#> 2 37.1428571428571 45.7142857142857 LDA for subpop 2 in target mixedpop2
#> 3                0 48.8721804511278 LDA for subpop 3 in target mixedpop2
#> 4               15             47.5 LDA for subpop 4 in target mixedpop2

#summary results Lasso to show the percent of cells classified as cells
#belonging by Lasso classifier
summary_prediction_lasso(LSOLDA_dat=LSOLDA_dat, nPredSubpop = row_cluster)
#>                 V1               V2                                      names
#> 1 44.9197860962567 73.7967914438503 ElasticNet for subpop1 in target mixedpop2
#> 2 88.5714285714286               85 ElasticNet for subpop2 in target mixedpop2
#> 3 52.6315789473684 26.3157894736842 ElasticNet for subpop3 in target mixedpop2
#> 4             82.5               60 ElasticNet for subpop4 in target mixedpop2

# summary maximum deviance explained by the model during the model training
summary_deviance(object = LSOLDA_dat)
#> $allDeviance
#> [1] "38.97" "67.06"
#> 
#> $DeviMax
#>           dat_DE$Dfd          Deviance           DEgenes
#> 1                  0             67.06    genes_cluster1
#> 2                  1             67.06    genes_cluster1
#> 3                  2             67.06    genes_cluster1
#> 4                  3             67.06    genes_cluster1
#> 5                  5             67.06    genes_cluster1
#> 6                  7             67.06    genes_cluster1
#> 7                  8             67.06    genes_cluster1
#> 8                 12             67.06    genes_cluster1
#> 9                 14             67.06    genes_cluster1
#> 10                15             67.06    genes_cluster1
#> 11                20             67.06    genes_cluster1
#> 12                22             67.06    genes_cluster1
#> 13                25             67.06    genes_cluster1
#> 14                29             67.06    genes_cluster1
#> 15                30             67.06    genes_cluster1
#> 16                33             67.06    genes_cluster1
#> 17                35             67.06    genes_cluster1
#> 18                39             67.06    genes_cluster1
#> 19                41             67.06    genes_cluster1
#> 20 remaining DEgenes remaining DEgenes remaining DEgenes
#> 
#> $LassoGenesMax
#> NULL

# summary accuracy to check the model accuracy in the leave-out test set
summary_accuracy(object = LSOLDA_dat)
#> [1] 65.62500 70.98214

4.3 Plot the relationship between clusters in one sample

Here we look at one example use case to find relationship between clusters within one sample or between two sample

#run prediction for 3 clusters
cluster_mixedpop1 <- colData(mixedpop1)[,1]
cluster_mixedpop2 <- colData(mixedpop2)[,1]
#cluster_mixedpop2 <- as.numeric(as.vector(colData(mixedpop2)[,1]))

c_selectID <- 1
#top 200 gene markers distinguishing cluster 1
genes = DEgenes$DE_Subpop1vsRemaining$id[1:200]

LSOLDA_dat1 <- bootstrap_prediction(nboots = 2, mixedpop1 = mixedpop2,
                        mixedpop2 = mixedpop2, genes=genes, c_selectID, listData =list(),
                        cluster_mixedpop1 = cluster_mixedpop2,
                        cluster_mixedpop2 = cluster_mixedpop2)

c_selectID <- 2
genes = DEgenes$DE_Subpop2vsRemaining$id[1:200]

LSOLDA_dat2 <- bootstrap_prediction(nboots = 2,mixedpop1 = mixedpop2,
                        mixedpop2 = mixedpop2, genes=genes, c_selectID, listData =list(),
                        cluster_mixedpop1 = cluster_mixedpop2,
                        cluster_mixedpop2 = cluster_mixedpop2)

c_selectID <- 3
genes = DEgenes$DE_Subpop3vsRemaining$id[1:200]
LSOLDA_dat3 <- bootstrap_prediction(nboots = 2,mixedpop1 = mixedpop2,
                        mixedpop2 = mixedpop2, genes=genes, c_selectID, listData =list(),
                        cluster_mixedpop1 = cluster_mixedpop2,
                        cluster_mixedpop2 = cluster_mixedpop2)

c_selectID <- 4
genes = DEgenes$DE_Subpop4vsRemaining$id[1:200]
LSOLDA_dat4 <- bootstrap_prediction(nboots = 2,mixedpop1 = mixedpop2,
                        mixedpop2 = mixedpop2, genes=genes, c_selectID, listData =list(),
                        cluster_mixedpop1 = cluster_mixedpop2,
                        cluster_mixedpop2 = cluster_mixedpop2)


#prepare table input for sankey plot

LASSO_C1S2  <- reformat_LASSO(c_selectID=1, mp_selectID = 2,
                             LSOLDA_dat=LSOLDA_dat1,
                             nPredSubpop = length(unique(colData(mixedpop2)[,1])),
                             Nodes_group ="#7570b3")

LASSO_C2S2  <- reformat_LASSO(c_selectID=2, mp_selectID =2,
                             LSOLDA_dat=LSOLDA_dat2,
                             nPredSubpop = length(unique(colData(mixedpop2)[,1])),
                             Nodes_group ="#1b9e77")

LASSO_C3S2  <- reformat_LASSO(c_selectID=3, mp_selectID =2,
                             LSOLDA_dat=LSOLDA_dat3,
                             nPredSubpop = length(unique(colData(mixedpop2)[,1])),
                             Nodes_group ="#e7298a")

LASSO_C4S2  <- reformat_LASSO(c_selectID=4, mp_selectID =2,
                             LSOLDA_dat=LSOLDA_dat4,
                             nPredSubpop = length(unique(colData(mixedpop2)[,1])),
                             Nodes_group ="#00FFFF")

combined <- rbind(LASSO_C1S2,LASSO_C2S2,LASSO_C3S2, LASSO_C4S2 )
combined <- combined[is.na(combined$Value) != TRUE,]

nboots = 2
#links: source, target, value
#source: node, nodegroup
combined_D3obj <-list(Nodes=combined[,(nboots+3):(nboots+4)],
                     Links=combined[,c((nboots+2):(nboots+1),ncol(combined))])

library(networkD3)

Node_source <- as.vector(sort(unique(combined_D3obj$Links$Source)))
Node_target <- as.vector(sort(unique(combined_D3obj$Links$Target)))
Node_all <-unique(c(Node_source, Node_target))

#assign IDs for Source (start from 0)
Source <-combined_D3obj$Links$Source
Target <- combined_D3obj$Links$Target

for(i in 1:length(Node_all)){
   Source[Source==Node_all[i]] <-i-1
   Target[Target==Node_all[i]] <-i-1
}
# 
combined_D3obj$Links$Source <- as.numeric(Source)
combined_D3obj$Links$Target <- as.numeric(Target)
combined_D3obj$Links$LinkColor <- combined$NodeGroup

#prepare node info
node_df <-data.frame(Node=Node_all)
node_df$id <-as.numeric(c(0, 1:(length(Node_all)-1)))

suppressMessages(library(dplyr))
Color <- combined %>% count(Node, color=NodeGroup) %>% select(2)
node_df$color <- Color$color

suppressMessages(library(networkD3))
p1<-sankeyNetwork(Links =combined_D3obj$Links, Nodes = node_df,
                 Value = "Value", NodeGroup ="color", LinkGroup = "LinkColor", NodeID="Node",
                 Source="Source", Target="Target", fontSize = 22 )
p1

4.3 Plot the relationship between clusters in two samples

Here we look at one example use case to find relationship between clusters within one sample or between two sample

#run prediction for 3 clusters
cluster_mixedpop1 <- colData(mixedpop1)[,1]
cluster_mixedpop2 <- colData(mixedpop2)[,1]
row_cluster <-length(unique(colData(mixedpop2)[,1]))

c_selectID <- 1
#top 200 gene markers distinguishing cluster 1
genes = DEgenes$DE_Subpop1vsRemaining$id[1:200]
LSOLDA_dat1 <- bootstrap_prediction(nboots = 2, mixedpop1 = mixedpop1,
                        mixedpop2 = mixedpop2, genes=genes, c_selectID, listData =list(),
                        cluster_mixedpop1 = cluster_mixedpop1,
                        cluster_mixedpop2 = cluster_mixedpop2)


c_selectID <- 2
genes = DEgenes$DE_Subpop2vsRemaining$id[1:200]
LSOLDA_dat2 <- bootstrap_prediction(nboots = 2,mixedpop1 = mixedpop1,
                        mixedpop2 = mixedpop2, genes=genes, c_selectID, listData =list(),
                        cluster_mixedpop1 = cluster_mixedpop1,
                        cluster_mixedpop2 = cluster_mixedpop2)

c_selectID <- 3
genes = DEgenes$DE_Subpop3vsRemaining$id[1:200]
LSOLDA_dat3 <- bootstrap_prediction(nboots = 2,mixedpop1 = mixedpop1,
                        mixedpop2 = mixedpop2, genes=genes, c_selectID, listData =list(),
                        cluster_mixedpop1 = cluster_mixedpop1,
                        cluster_mixedpop2 = cluster_mixedpop2)

#prepare table input for sankey plot

LASSO_C1S1  <- reformat_LASSO(c_selectID=1, mp_selectID = 1,
                             LSOLDA_dat=LSOLDA_dat1, nPredSubpop = row_cluster, Nodes_group = "#7570b3")

LASSO_C2S1  <- reformat_LASSO(c_selectID=2, mp_selectID = 1,
                             LSOLDA_dat=LSOLDA_dat2, nPredSubpop = row_cluster, Nodes_group = "#1b9e77")

LASSO_C3S1  <- reformat_LASSO(c_selectID=3, mp_selectID = 1,
                             LSOLDA_dat=LSOLDA_dat3, nPredSubpop = row_cluster, Nodes_group = "#e7298a")


combined <- rbind(LASSO_C1S1,LASSO_C2S1,LASSO_C3S1)

nboots = 2
#links: source, target, value
#source: node, nodegroup
combined_D3obj <-list(Nodes=combined[,(nboots+3):(nboots+4)],
                     Links=combined[,c((nboots+2):(nboots+1),ncol(combined))])
combined <- combined[is.na(combined$Value) != TRUE,]


library(networkD3)

Node_source <- as.vector(sort(unique(combined_D3obj$Links$Source)))
Node_target <- as.vector(sort(unique(combined_D3obj$Links$Target)))
Node_all <-unique(c(Node_source, Node_target))

#assign IDs for Source (start from 0)
Source <-combined_D3obj$Links$Source
Target <- combined_D3obj$Links$Target

for(i in 1:length(Node_all)){
   Source[Source==Node_all[i]] <-i-1
   Target[Target==Node_all[i]] <-i-1
}

combined_D3obj$Links$Source <- as.numeric(Source)
combined_D3obj$Links$Target <- as.numeric(Target)
combined_D3obj$Links$LinkColor <- combined$NodeGroup

#prepare node info
node_df <-data.frame(Node=Node_all)
node_df$id <-as.numeric(c(0, 1:(length(Node_all)-1)))

suppressMessages(library(dplyr))
n <- length(unique(node_df$Node))
getPalette = colorRampPalette(RColorBrewer::brewer.pal(9, "Set1"))
Color = getPalette(n)
node_df$color <- Color
suppressMessages(library(networkD3))
p1<-sankeyNetwork(Links =combined_D3obj$Links, Nodes = node_df,
                 Value = "Value", NodeGroup ="color", LinkGroup = "LinkColor",
                 NodeID="Node", Source="Source", Target="Target", fontSize = 22)
p1
devtools::session_info()
#> ─ Session info ───────────────────────────────────────────────────────────────
#>  setting  value                       
#>  version  R version 4.0.3 (2020-10-10)
#>  os       Ubuntu 18.04.5 LTS          
#>  system   x86_64, linux-gnu           
#>  ui       X11                         
#>  language (EN)                        
#>  collate  C                           
#>  ctype    en_US.UTF-8                 
#>  tz       America/New_York            
#>  date     2020-10-27                  
#> 
#> ─ Packages ───────────────────────────────────────────────────────────────────
#>  package              * version    date       lib source        
#>  annotate               1.68.0     2020-10-27 [2] Bioconductor  
#>  AnnotationDbi        * 1.52.0     2020-10-27 [2] Bioconductor  
#>  assertthat             0.2.1      2019-03-21 [2] CRAN (R 4.0.3)
#>  backports              1.1.10     2020-09-15 [2] CRAN (R 4.0.3)
#>  Biobase              * 2.50.0     2020-10-27 [2] Bioconductor  
#>  BiocGenerics         * 0.36.0     2020-10-27 [2] Bioconductor  
#>  BiocManager            1.30.10    2019-11-16 [2] CRAN (R 4.0.3)
#>  BiocParallel           1.24.0     2020-10-27 [2] Bioconductor  
#>  bit                    4.0.4      2020-08-04 [2] CRAN (R 4.0.3)
#>  bit64                  4.0.5      2020-08-30 [2] CRAN (R 4.0.3)
#>  bitops                 1.0-6      2013-08-17 [2] CRAN (R 4.0.3)
#>  blob                   1.2.1      2020-01-20 [2] CRAN (R 4.0.3)
#>  callr                  3.5.1      2020-10-13 [2] CRAN (R 4.0.3)
#>  caret                * 6.0-86     2020-03-20 [2] CRAN (R 4.0.3)
#>  checkmate              2.0.0      2020-02-06 [2] CRAN (R 4.0.3)
#>  class                  7.3-17     2020-04-26 [2] CRAN (R 4.0.3)
#>  cli                    2.1.0      2020-10-12 [2] CRAN (R 4.0.3)
#>  clusterProfiler      * 3.18.0     2020-10-27 [2] Bioconductor  
#>  codetools              0.2-16     2018-12-24 [2] CRAN (R 4.0.3)
#>  colorspace             1.4-1      2019-03-18 [2] CRAN (R 4.0.3)
#>  cowplot                1.1.0      2020-09-08 [2] CRAN (R 4.0.3)
#>  crayon                 1.3.4      2017-09-16 [2] CRAN (R 4.0.3)
#>  data.table             1.13.2     2020-10-19 [2] CRAN (R 4.0.3)
#>  DBI                    1.1.0      2019-12-15 [2] CRAN (R 4.0.3)
#>  DelayedArray           0.16.0     2020-10-27 [2] Bioconductor  
#>  dendextend             1.14.0     2020-08-26 [2] CRAN (R 4.0.3)
#>  desc                   1.2.0      2018-05-01 [2] CRAN (R 4.0.3)
#>  DESeq                * 1.42.0     2020-10-27 [2] Bioconductor  
#>  devtools               2.3.2      2020-09-18 [2] CRAN (R 4.0.3)
#>  digest                 0.6.27     2020-10-24 [2] CRAN (R 4.0.3)
#>  DO.db                  2.9        2020-10-12 [2] Bioconductor  
#>  DOSE                 * 3.16.0     2020-10-27 [2] Bioconductor  
#>  downloader             0.4        2015-07-09 [2] CRAN (R 4.0.3)
#>  dplyr                * 1.0.2      2020-08-18 [2] CRAN (R 4.0.3)
#>  dynamicTreeCut       * 1.63-1     2016-03-11 [2] CRAN (R 4.0.3)
#>  e1071                  1.7-4      2020-10-14 [2] CRAN (R 4.0.3)
#>  ellipsis               0.3.1      2020-05-15 [2] CRAN (R 4.0.3)
#>  enrichplot             1.10.0     2020-10-27 [2] Bioconductor  
#>  evaluate               0.14       2019-05-28 [2] CRAN (R 4.0.3)
#>  fansi                  0.4.1      2020-01-08 [2] CRAN (R 4.0.3)
#>  farver                 2.0.3      2020-01-16 [2] CRAN (R 4.0.3)
#>  fastcluster            1.1.25     2018-06-07 [2] CRAN (R 4.0.3)
#>  fastmatch              1.1-0      2017-01-28 [2] CRAN (R 4.0.3)
#>  fgsea                  1.16.0     2020-10-27 [2] Bioconductor  
#>  foreach                1.5.1      2020-10-15 [2] CRAN (R 4.0.3)
#>  fs                     1.5.0      2020-07-31 [2] CRAN (R 4.0.3)
#>  genefilter             1.72.0     2020-10-27 [2] Bioconductor  
#>  geneplotter            1.68.0     2020-10-27 [2] Bioconductor  
#>  generics               0.0.2      2018-11-29 [2] CRAN (R 4.0.3)
#>  GenomeInfoDb         * 1.26.0     2020-10-27 [2] Bioconductor  
#>  GenomeInfoDbData       1.2.4      2020-10-12 [2] Bioconductor  
#>  GenomicRanges        * 1.42.0     2020-10-27 [2] Bioconductor  
#>  ggforce                0.3.2      2020-06-23 [2] CRAN (R 4.0.3)
#>  ggplot2              * 3.3.2      2020-06-19 [2] CRAN (R 4.0.3)
#>  ggraph                 2.0.3      2020-05-20 [2] CRAN (R 4.0.3)
#>  ggrepel                0.8.2      2020-03-08 [2] CRAN (R 4.0.3)
#>  glmnet                 4.0-2      2020-06-16 [2] CRAN (R 4.0.3)
#>  glue                   1.4.2      2020-08-27 [2] CRAN (R 4.0.3)
#>  GO.db                  3.12.0     2020-10-15 [2] Bioconductor  
#>  GOSemSim               2.16.0     2020-10-27 [2] Bioconductor  
#>  gower                  0.2.2      2020-06-23 [2] CRAN (R 4.0.3)
#>  graph                  1.68.0     2020-10-27 [2] Bioconductor  
#>  graphite               1.36.0     2020-10-27 [2] Bioconductor  
#>  graphlayouts           0.7.1      2020-10-26 [2] CRAN (R 4.0.3)
#>  gridExtra              2.3        2017-09-09 [2] CRAN (R 4.0.3)
#>  gtable                 0.3.0      2019-03-25 [2] CRAN (R 4.0.3)
#>  htmltools              0.5.0      2020-06-16 [2] CRAN (R 4.0.3)
#>  htmlwidgets            1.5.2      2020-10-03 [2] CRAN (R 4.0.3)
#>  httr                   1.4.2      2020-07-20 [2] CRAN (R 4.0.3)
#>  igraph                 1.2.6      2020-10-06 [2] CRAN (R 4.0.3)
#>  ipred                  0.9-9      2019-04-28 [2] CRAN (R 4.0.3)
#>  IRanges              * 2.24.0     2020-10-27 [2] Bioconductor  
#>  iterators              1.0.13     2020-10-15 [2] CRAN (R 4.0.3)
#>  jsonlite               1.7.1      2020-09-07 [2] CRAN (R 4.0.3)
#>  knitr                  1.30       2020-09-22 [2] CRAN (R 4.0.3)
#>  labeling               0.4.2      2020-10-20 [2] CRAN (R 4.0.3)
#>  lattice              * 0.20-41    2020-04-02 [2] CRAN (R 4.0.3)
#>  lava                   1.6.8      2020-09-26 [2] CRAN (R 4.0.3)
#>  lifecycle              0.2.0      2020-03-06 [2] CRAN (R 4.0.3)
#>  locfit               * 1.5-9.4    2020-03-25 [2] CRAN (R 4.0.3)
#>  lubridate              1.7.9      2020-06-08 [2] CRAN (R 4.0.3)
#>  magrittr               1.5        2014-11-22 [2] CRAN (R 4.0.3)
#>  MASS                   7.3-53     2020-09-09 [2] CRAN (R 4.0.3)
#>  Matrix                 1.2-18     2019-11-27 [2] CRAN (R 4.0.3)
#>  MatrixGenerics       * 1.2.0      2020-10-27 [2] Bioconductor  
#>  matrixStats          * 0.57.0     2020-09-25 [2] CRAN (R 4.0.3)
#>  memoise                1.1.0      2017-04-21 [2] CRAN (R 4.0.3)
#>  ModelMetrics           1.2.2.2    2020-03-17 [2] CRAN (R 4.0.3)
#>  munsell                0.5.0      2018-06-12 [2] CRAN (R 4.0.3)
#>  networkD3            * 0.4        2017-03-18 [2] CRAN (R 4.0.3)
#>  nlme                   3.1-150    2020-10-24 [2] CRAN (R 4.0.3)
#>  nnet                   7.3-14     2020-04-26 [2] CRAN (R 4.0.3)
#>  org.Hs.eg.db         * 3.12.0     2020-10-15 [2] Bioconductor  
#>  pillar                 1.4.6      2020-07-10 [2] CRAN (R 4.0.3)
#>  pkgbuild               1.1.0      2020-07-13 [2] CRAN (R 4.0.3)
#>  pkgconfig              2.0.3      2019-09-22 [2] CRAN (R 4.0.3)
#>  pkgload                1.1.0      2020-05-29 [2] CRAN (R 4.0.3)
#>  plyr                   1.8.6      2020-03-03 [2] CRAN (R 4.0.3)
#>  polyclip               1.10-0     2019-03-14 [2] CRAN (R 4.0.3)
#>  prettyunits            1.1.1      2020-01-24 [2] CRAN (R 4.0.3)
#>  pROC                   1.16.2     2020-03-19 [2] CRAN (R 4.0.3)
#>  processx               3.4.4      2020-09-03 [2] CRAN (R 4.0.3)
#>  prodlim                2019.11.13 2019-11-17 [2] CRAN (R 4.0.3)
#>  ps                     1.4.0      2020-10-07 [2] CRAN (R 4.0.3)
#>  purrr                  0.3.4      2020-04-17 [2] CRAN (R 4.0.3)
#>  qvalue                 2.22.0     2020-10-27 [2] Bioconductor  
#>  R6                     2.4.1      2019-11-12 [2] CRAN (R 4.0.3)
#>  rappdirs               0.3.1      2016-03-28 [2] CRAN (R 4.0.3)
#>  RColorBrewer           1.1-2      2014-12-07 [2] CRAN (R 4.0.3)
#>  Rcpp                   1.0.5      2020-07-06 [2] CRAN (R 4.0.3)
#>  RcppArmadillo          0.10.1.0.0 2020-10-20 [2] CRAN (R 4.0.3)
#>  RcppParallel           5.0.2      2020-06-24 [2] CRAN (R 4.0.3)
#>  RCurl                  1.98-1.2   2020-04-18 [2] CRAN (R 4.0.3)
#>  reactome.db            1.74.0     2020-10-20 [2] Bioconductor  
#>  ReactomePA           * 1.34.0     2020-10-27 [2] Bioconductor  
#>  recipes                0.1.14     2020-10-17 [2] CRAN (R 4.0.3)
#>  remotes                2.2.0      2020-07-21 [2] CRAN (R 4.0.3)
#>  reshape2               1.4.4      2020-04-09 [2] CRAN (R 4.0.3)
#>  rlang                  0.4.8      2020-10-08 [2] CRAN (R 4.0.3)
#>  rmarkdown              2.5        2020-10-21 [2] CRAN (R 4.0.3)
#>  rpart                  4.1-15     2019-04-12 [2] CRAN (R 4.0.3)
#>  rprojroot              1.3-2      2018-01-03 [2] CRAN (R 4.0.3)
#>  RSQLite                2.2.1      2020-09-30 [2] CRAN (R 4.0.3)
#>  Rtsne                  0.15       2018-11-10 [2] CRAN (R 4.0.3)
#>  rvcheck                0.1.8      2020-03-01 [2] CRAN (R 4.0.3)
#>  S4Vectors            * 0.28.0     2020-10-27 [2] Bioconductor  
#>  scales                 1.1.1      2020-05-11 [2] CRAN (R 4.0.3)
#>  scatterpie             0.1.5      2020-09-09 [2] CRAN (R 4.0.3)
#>  scGPS                * 1.4.0      2020-10-27 [1] Bioconductor  
#>  sessioninfo            1.1.1      2018-11-05 [2] CRAN (R 4.0.3)
#>  shadowtext             0.0.7      2019-11-06 [2] CRAN (R 4.0.3)
#>  shape                  1.4.5      2020-09-13 [2] CRAN (R 4.0.3)
#>  SingleCellExperiment * 1.12.0     2020-10-27 [2] Bioconductor  
#>  stringi                1.5.3      2020-09-09 [2] CRAN (R 4.0.3)
#>  stringr                1.4.0      2019-02-10 [2] CRAN (R 4.0.3)
#>  SummarizedExperiment * 1.20.0     2020-10-27 [2] Bioconductor  
#>  survival               3.2-7      2020-09-28 [2] CRAN (R 4.0.3)
#>  testthat               2.3.2      2020-03-02 [2] CRAN (R 4.0.3)
#>  tibble                 3.0.4      2020-10-12 [2] CRAN (R 4.0.3)
#>  tidygraph              1.2.0      2020-05-12 [2] CRAN (R 4.0.3)
#>  tidyr                  1.1.2      2020-08-27 [2] CRAN (R 4.0.3)
#>  tidyselect             1.1.0      2020-05-11 [2] CRAN (R 4.0.3)
#>  timeDate               3043.102   2018-02-21 [2] CRAN (R 4.0.3)
#>  tweenr                 1.0.1      2018-12-14 [2] CRAN (R 4.0.3)
#>  usethis                1.6.3      2020-09-17 [2] CRAN (R 4.0.3)
#>  vctrs                  0.3.4      2020-08-29 [2] CRAN (R 4.0.3)
#>  viridis                0.5.1      2018-03-29 [2] CRAN (R 4.0.3)
#>  viridisLite            0.3.0      2018-02-01 [2] CRAN (R 4.0.3)
#>  withr                  2.3.0      2020-09-22 [2] CRAN (R 4.0.3)
#>  xfun                   0.18       2020-09-29 [2] CRAN (R 4.0.3)
#>  XML                    3.99-0.5   2020-07-23 [2] CRAN (R 4.0.3)
#>  xtable                 1.8-4      2019-04-21 [2] CRAN (R 4.0.3)
#>  XVector                0.30.0     2020-10-27 [2] Bioconductor  
#>  yaml                   2.2.1      2020-02-01 [2] CRAN (R 4.0.3)
#>  zlibbioc               1.36.0     2020-10-27 [2] Bioconductor  
#> 
#> [1] /tmp/RtmpR6nNRV/Rinst3b47c48d809
#> [2] /home/biocbuild/bbs-3.12-bioc/R/library