bump_loglik_cellcycle: Loglikelihood of of sinusoidal genes in cell cycle

Description Usage Arguments

View source: R/bump_loglik_cellcycle.R

Description

Computes the loglikelihood of all the cells in the cycle under sinusoidal gene expression patterns.

@author Kushal K Dey

@export

@examples G <- 500; num_cells <- 400; amp_genes <- rep(10, G); phi_genes <- runif(G, 0, 2*pi) sigma_genes <- rchisq(G, 4); cell_times_sim <- sample(seq(0,2*pi, 2*pi/(num_cells-1)), num_cells, replace=FALSE); cycle_data <- sim_sinusoidal_cycle(G, amp_genes, phi_genes, sigma_genes, cell_times_sim);

loglik_cell_cycle(cycle_data, cell_times, amp_genes, phi_genes, sigma_genes)

Usage

1
bump_loglik_cellcycle(data, cell_times, amp, phi, sigma)

Arguments

cell_times

: a N x 1 vector of cell times

cycle_data:

a N x G matrix, where N is number of cells, G number of genes

amp:

the amplitude vector (G x 1) over the genes

phi:

the G x 1 vector of phase values over genes

sigma:

the G x 1 vector of gene variation

freq:

frequency of the signal of the gene. Defaults to 1.


kkdey/cellcycleR documentation built on May 20, 2019, 10:33 a.m.