In compbiomed/singleCellTK: Comprehensive and Interactive Analysis of Single Cell RNA-Seq Data

knitr::opts_chunk$set(warning = FALSE)

````{=html}

## Introduction

Non-linear 2D embedding algorithms can be run through Single-Cell Toolkit (SCTK) using both interactive shiny application and R console. For the interactive analysis, the toolkit offers a streamlined workflow to both compute the embedding and then visualize the results. For the console analysis, SCTK offers a generic function `runDimReduce()` to compute the 2D embedding.

Methods available with the toolkit include tSNE and UMAP from *Seurat* [@Butler2018][@Stuart2019][@Satija2017][@Hao2021] package, tSNE from *Rtsne* [@Maaten2014][@Maaten2008] package and UMAP from *scater* [@McCarthy2017] package. 

To view detailed instructions on how to use these methods, please select "Interactive Analysis", from the tabs below, for the usage in Shiny application or "Console Analysis" for using these methods in R console: 

## Workflow Guide

````{=html}
<div class="tab">
  <button class="tablinks" onclick="openTab(event, 'interactive')" id="ia-button">Interactive Analysis</button>
  <button class="tablinks" onclick="openTab(event, 'console')" id="console-button">Console Analysis</button>
</div>

<div id="interactive" class="tabcontent">

In general, the UI offers options for selection of data items and choice of parameters on the left side, and a visualization panel on the right side of the interface.

1. To begin with, click on the "Feature Selection & Dimensionality Reduction" tab from the top menu. This workflow assumes that before proceeding towards computation of 2D embedding, data has been uploaded, filtered, normalized and dimension reduction has been performed through the preceding tabs.

\

1. Select "Embedding" sub-tab:

\

1. Select an input matrix to use with the tSNE or UMAP algorithm. A dimension reduction is recommended and can be obtained from Dimension Reduction tab.
2. Select a specific tSNE or UMAP method for computation. Available methods include tSNE and UMAP from Seurat package, tSNE from Rtsne package and UMAP from scater package.
3. Specify a name for the output embedding.
4. (Optional) If selected input in step 3. is an feature expression matrix, users can specify preprocessing here before the embedding is computed. HVG list refers to a feature subset which could be obtained in Feature Selection tab. Preprocessing is only available for scater UMAP and Rtsne tSNE methods.
5. Number of top dimensions (components) to use from the input dimension reduction (or the PCA generated with preprocessing specified by step 6.). Methods for determining this number have been provided in Dimension Reduction tab. And we recommend using the same number for the dimension reduction used for clustering, in order to get the best visualization of the clustering result.
6. Embedding method specific parameters. These vary when different methods are selected in step 3.
7. Set the seed of randomized computation process. A fixed seed generates reproducible result each run.
8. Start computation of selected method with the specified parameters.
9. Once computation is complete, a scatter plot of the cells on the selected 2D embedding is shown in the panel with additional options available in step 12.

\

1. Users can select another computed embedding to visualize through dropdown box.

\

````{=html}