run.CSIDE.nonparam: Runs CSIDE on a 'RCTD' object to detect nonparametric smooth...

In dmcable/RCTD: SpatialeXpressionR: Cell type identification and cell type-specific differential expression in spatial transcriptomics

Runs CSIDE on a RCTD object to detect nonparametric smooth gene expression patterns

Description

Identifies cell type specific smooth gene expression patterns. The design matrix contains thin plate spline basis functions spanning the space of smooth functions. Uses maximum likelihood estimation to estimate DE and standard errors for each gene and each cell type. Selects genes with significant nonzero DE.

Usage

run.CSIDE.nonparam(
  myRCTD,
  df = 15,
  barcodes = NULL,
  cell_types = NULL,
  cell_type_threshold = 125,
  gene_threshold = 5e-05,
  doublet_mode = T,
  weight_threshold = NULL,
  sigma_gene = T,
  PRECISION.THRESHOLD = 0.05,
  cell_types_present = NULL,
  fdr = 0.01,
  test_genes_sig = T,
  logs = F,
  test_error = F
)

Arguments

myRCTD	an RCTD object with annotated cell types e.g. from the run.RCTD function.
df	(default 15) the degrees of freedom, or number of basis functions to be used in the model.
barcodes	the barcodes, or pixel names, of the SpatialRNA object to be used when fitting the model.
cell_types	the cell types used for CSIDE. If null, cell types will be chosen with aggregate occurences of at least 'cell_type_threshold', as aggregated by aggregate_cell_types
cell_type_threshold	(default 125) min occurence of number of cells for each cell type to be used, as aggregated by aggregate_cell_types
gene_threshold	(default 5e-5) minimum average normalized expression required for selecting genes
doublet_mode	(default TRUE) if TRUE, uses RCTD doublet mode weights. Otherwise, uses RCTD full mode weights
weight_threshold	(default NULL) the threshold of total normalized weights across all cell types in cell_types per pixel to be included in the model. Default 0.99 for doublet_mode or 0.8 for full_mode.
sigma_gene	(default TRUE) if TRUE, fits gene specific overdispersion parameter. If FALSE, overdispersion parameter is same across all genes.
PRECISION.THRESHOLD	(default 0.05) for checking for convergence, the maximum parameter change per algorithm step
cell_types_present	cell types (a superset of 'cell_types') to be considered as occuring often enough to consider for gene expression contamination during the step filtering out marker genes of other cell types.
fdr	(default 0.01) false discovery rate for hypothesis testing
test_genes_sig	(default TRUE) logical controlling whether genes will be tested for significance
logs	(default FALSE) if TRUE, writes progress to logs/de_logs.txt
test_error	(default FALSE) if TRUE, exits after testing for error messages without running CSIDE. If set to TRUE, this can be used to quickly evaluate if CSIDE will run without error.

Value

an RCTD object containing the results of the CSIDE algorithm. Contains objects de_results, which contain the results of the CSIDE algorithm including 'gene_fits', which contains the results of fits on individual genes, in addition 'sig_gene_list', a list, for each cell type, of significant genes detected by CSIDE. Additionally, the object contains 'internal_vars_de' a list of variables that are used internally by CSIDE

dmcable/RCTD documentation built on Feb. 24, 2024, 11:03 p.m.