pseudotimeDE | R Documentation |
Test if one gene is differentially expressed along pseudotime.
pseudotimeDE(
gene,
ori.tbl,
sub.tbl,
mat,
assay.use = "counts",
model = c("nb", "zinb", "gaussian", "auto", "qgam"),
k = 6,
knots = c(0:5/5),
fix.weight = TRUE,
aicdiff = 10,
seed = 123,
quant = 0.5,
usebam = FALSE,
formula = NULL,
seurat.assay = "RNA"
)
gene |
A string of gene name. It should be one of the row names in sce. |
ori.tbl |
A tibble or dataframe which contains the original cells and pseudotime as two columns. |
sub.tbl |
A list of tibbles or dataframes where each is the fit of a subsample. Each element is the same format as |
mat |
The input expression data. It can be: (1) A SingleCellExperment object which contain the expression data; (2) An matrix; (3) A Seurat object which contain the expression data. Its row names should be genes and col names should be cells. |
assay.use |
The |
model |
A string of the model name. One of |
k |
A integer of the basis dimension. Default is 6. The reults are usually robust to different k; we recommend to use k from 5 to 10. |
knots |
A numeric vector of the location of knots. Default is evenly distributed between 0 to 1. For instance, if your k = 6, and your range is [0, 10], then the position of knots should be |
fix.weight |
A logic variable indicating if the ZINB-GAM will use the zero weights from the original model. Used for saving time since ZINB-GAM is computationally intense. |
aicdiff |
A numeric variable of the threshold of model selection. Only works when |
seed |
A numeric variable of the random seed. It mainly affects the parametricfitting of null distribution. |
quant |
The quantile of interest for quantile regression (qgam), range from 0 to 1, default as 0.5. |
usebam |
A logical variable. If use |
formula |
An (optional) custom formula to be passed to |
seurat.assay |
The |
A list with the components:
fix.pv
The p-value assuming the pseudotime is fixed
emp.pv
The permutation p-value
para.pv
The permutation p-value by fitting a parametric null distribution
ad.pv
P-value of Anderson-Darling test on comparing null distribution and its parametric fit
rank
The estimated effect degree of freedom of the original model
gam.fit
The fitted gam model on original data
zinf
Whether the model is zero inflated
aic
The AIC of the orignial model
expv.quantile
Quantiles of the log counts plus 1
expv.mean
Mean of the log counts plus 1
expv.zero
Zero proportions of counts
Dongyuan Song, Shiyu Ma
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.