sigmanoise_sampling function

Fitting Residule Variance Sampling Function

Function 'sigmanoise_sampling' estimates the variance of overal gene expression fitting residuals, according to an inverse-gamma distribution.

sigmanoise_sampling(Y, C, A, X, base_line, C_prior, sigma_A, sigma_baseline,
  sigma_X, alpha, beta)

Arguments

• Y: gene expression data matrix
• C: sampled TF-gene binding network
• A: sampled regulatory strength matrix
• X: sampled transcription factor activity matrix
• base_line: sampled gene expression baseline activity
• C_prior: prior TF-gene binding network
• sigma_A: variance of regulatory strength
• sigma_baseline: variance of gene expression baseline activity
• sigma_X: variance of transcription factor activity
• alpha: hyper-parameter for inverse-gamma distribution
• beta: hyper-parameter for inverse-gamma distribution