View source: R/sin_loglik_cellcycle.R
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)
1 | sin_loglik_cellcycle(cycle_data, cell_times, amp, phi, sigma, freq = 1)
|
cycle_data |
a N x G matrix, where N is number of cells, G number of genes |
cell_times |
a N x 1 vector of cell times |
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. |
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.