get_log_p_D_KL_continuous function

Estimates the significance of the observed Kullback-Leibler divergence by comparing to randomizations for the continuous version of haystack.

get_log_p_D_KL_continuous(
  D_KL.observed,
  D_KL.randomized,
  all.coeffVar,
  train.coeffVar,
  output.dir = NULL,
  spline.method = "ns"
)

Arguments

• D_KL.observed: A vector of observed Kullback-Leibler divergences.
• D_KL.randomized: A matrix of Kullback-Leibler divergences of randomized datasets.
• all.coeffVar: Coefficients of variation of all genes. Used for fitting the Kullback-Leibler divergences.
• train.coeffVar: Coefficients of variation of genes that will be used for fitting the Kullback-Leibler divergences.
• output.dir: Optional parameter. Default is NULL. If not NULL, some files will be written to this directory.
• spline.method: Method to use for fitting splines "ns" (default): natural splines, "bs": B-splines.

Returns

A vector of log10 p values, not corrected for multiple testing using the Bonferroni correction.