squared_epsilon: Paired-sine model squared error term

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

Description

This function returns the squared error term as described in Leng et al. 2015, section Oscope: paired-sine model.

Usage

squared_epsilon(psi, x1, x2)

Arguments

psi	Phase shift between the gene expression of the two genes.
x1	Vector of gene expression values for gene 1.
x2	Vector of gene expression values for gene 2.

Details

The definition of the error term is:

ε_{1,2}^2 = ∑_s ( x_{1,s}^2 + x_{2,s}^2 - 2 x_{1,s} x_{2,s} cos(ψ) - sin(ψ)^2)^2

Value

Squared error term, i.e., ε_{1,2}^2.

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