Skip to contents

Generalized likelihood ratio test.

GLRT.Rd

Using beta-binomial to model the exon count and do generalized likelihood ratio test for two cell populations to test if they have the same psi distribution.

Usage

GLRT(
  exon_count1,
  gene_count1,
  exon_count2,
  gene_count2,
  iters = 100,
  psi_num = 500,
  p_thresh = 0.05
)

Arguments

exon_count1

a vector of exon count for cell population 1

gene_count1

a vector of gene count for cell population 1

exon_count2

a vector of exon count for cell population 2

gene_count2

a vector of gene count for cell population 2

iters

the number of iterations for bootstrap

psi_num

the sample size from two beta distributions to calculate the earth-moving distance

p_thresh

the threshold for the p-value

Value

a vector of p-value and parameters for the psi distributions in two cell populations if they have significant difference