plotGeneCount {tradeSeq} | R Documentation

## Plot gene expression in reduced dimension.

### Description

Plot the gene in reduced dimensional space.

### Usage

    plotGeneCount(
      curve,
      counts = NULL,
      gene = NULL,
      clusters = NULL,
      models = NULL,
      title = NULL
    )

### Arguments

`curve` | A `SlingshotDataSet` object. The output from trajectory inference using Slingshot.
`counts` | The count matrix, genes in rows and cells in columns.
`gene` | The name of gene for which you want to plot the count or the row number of that gene in the count matrix. Alternatively, one can specify the `clusters` argument.
`clusters` | The assignation of each cell to a cluster. Used to color the plot. Either `clusters` or `gene` and `counts` must be supplied.
`models` | The fitted GAMs, typically the output from `fitGAM`. Used to display the knots.
`title` | Title for the plot.

### Details

If both `gene` and `clusters` arguments are supplied, the plot will be colored according to gene count level.

### Value

A `ggplot` object

### Examples

    set.seed(97)
    library(slingshot)
    data(crv, package="tradeSeq")
    data(countMatrix, package="tradeSeq")
    rd <- slingshot::reducedDim(crv)
    cl <- kmeans(rd, centers = 7)$cluster
    lin <- slingshot::getLineages(rd, clusterLabels = cl, start.clus = 4)
    crv <- slingshot::getCurves(lin)
    counts <- as.matrix(countMatrix)
    gamList <- fitGAM(counts = counts,
     pseudotime = slingPseudotime(crv, na = FALSE),
     cellWeights = slingCurveWeights(crv))
    plotGeneCount(crv, counts, gene = "Mpo")

[Package tradeSeq version 1.3.13 Index]