Test gene expression changes along a TSCAN trajectory path

Source: R/runTSCAN.R

Wrapper for identifying genes with significant changes with respect to one of the TSCAN pseudotime paths

runTSCANDEG(inSCE, pathIndex, useAssay = "logcounts", discardCluster = NULL)

Arguments

inSCE

Input SingleCellExperiment object.

pathIndex

Path index for which the pseudotime values should be used. This corresponds to the terminal node of specific path from the root node to the terminal node. Run listTSCANTerminalNodes(inSCE) for available options.

useAssay

Character. The name of the assay to use for testing the expression change. Should be log-normalized. Default "logcounts"

discardCluster

Cluster(s) which are not of use or masks other interesting effects can be discarded. Default NULL.

Value

The input inSCE with results updated in metadata.

Author

Nida Pervaiz

Examples

data("mouseBrainSubsetSCE", package = "singleCellTK")
mouseBrainSubsetSCE <- runTSCAN(inSCE = mouseBrainSubsetSCE,
                                useReducedDim = "PCA_logcounts")
#> Fri Feb  2 14:23:09 2024 ... Running 'scran SNN clustering' with 'louvain' algorithm
#> Fri Feb  2 14:23:10 2024 ...   Identified 2 clusters
#> Fri Feb  2 14:23:10 2024 ... Running TSCAN to estimate pseudotime
#> Fri Feb  2 14:23:11 2024 ...   Clusters involved in path index 2 are: 1, 2
#> Fri Feb  2 14:23:11 2024 ...   Number of estimated paths is 1
terminalNodes <- listTSCANTerminalNodes(mouseBrainSubsetSCE)
mouseBrainSubsetSCE <- runTSCANDEG(inSCE = mouseBrainSubsetSCE,
                                   pathIndex = terminalNodes[1])