find_optimal_phase_shift: Optimal phase shift

In ramiromagno/oscillation: A Statistical Pipeline for Identifying Oscillatory Genes in Unsynchronized Single Cell RNA-Seq

Description

This function determines the optimal phase shift between two gene expression profiles x1 and x2 according to the paired-sine model, see Leng et al. 2015, section Oscope: paired-sine model.

Usage

find_optimal_phase_shift(x1, x2, psi_start = 0, method = "L-BFGS-B")

Arguments

x1	Vector of gene expression values for gene 1.
x2	Vector of gene expression values for gene 2.
psi_start	The phase shift (psi) starting value for the optimisation search.
method	The name of the optimisation method to be used by optim.

Value

A tibble where rows are for each pair of genes, and five columns: gene 1, gene 2, psi, se_epsilon and score.

ramiromagno/oscillation documentation built on April 20, 2020, 10:37 a.m.