gene_pseudotime_test: Test a single gene as a function of pseudotime

In kieranrcampbell/embeddr: Laplacian eigenmaps & principal curves for pseudotemporal ordering of single-cell RNA-seq data

Description

Test a single gene as a function of pseudotime

Usage

gene_pseudotime_test(sce, gene, full_model = NULL)

Arguments

sce	An object of class SCESet
gene	Name or index of the gene to be tested.
full_model	If a full pseudotime model has already been calculated, use that. Otherwise (if NULL, default), recalculate one on the fly.

Value

A p-value for the given gene

Examples

library(scater)
data('sc_example_counts') ; sce <- newSCESet(countData = sc_example_counts)
sce <- embeddr(sce)
sce <- fit_pseudotime(sce)
pval <- gene_pseudotime_test(sce, 1) # select first gene

kieranrcampbell/embeddr documentation built on May 20, 2019, 9:24 a.m.