Correlate two genes across samples

Computes Pearson or Spearman correlation between two genes in an expression matrix.

Usage

correlate_genes(expr_matrix, gene_x, gene_y, method = c("pearson", "spearman"))

Arguments

expr_matrix

Numeric matrix (genes x samples) of transformed expression values (e.g., VST)

gene_x

Gene name for x-axis

gene_y

Gene name for y-axis

method

Correlation method: "pearson" (default) or "spearman"

Value

List with components: gene_x, gene_y, method, estimate, p_value, n_samples, expr_x (numeric vector), expr_y (numeric vector)

See also

Other gene-correlation: correlate_genes_batch(), plot_gene_correlation()