Gibbs.post: Calculate the Gibbs posterior

In JustinaXie/NDATA: Network Assisted Analysis for De novo mutations

Calculate the Gibbs posterior

Description

Calculate the Gibbs posterior

Usage

Gibbs.post(
  theta_0,
  theta_1,
  nodesNeighbor,
  w,
  W = NULL,
  Y,
  N,
  mu,
  mcmc_samples = 1010,
  burnin = 10,
  thin = 1,
  initLabel = NULL
)

Arguments

theta_0	The estimated value of theta_0.
theta_1	The estimated value of theta_1.
nodesNeighbor	Neighbors of corresponding nodes in the network.
w	The parameter in the prior for network.
W	Another way to input the weight. For internal use.
Y	The vector of DNM counts for all genes.
N	The parameter in the prior for network.
mu	The vector of mutability for all genes.
mcmc_samples	The number of iterations for mcmc.
burnin	The number of burn-ins for mcmc.
thin	The number of thinning for mcmc.
initLabel	The vector of initial labels for all genes.

Value

The list containing S_mcmc the mcmc samples after burning-in and thinning, q_0 the posterior probability for calculating FDR, llk the likelihood for checking convergence

JustinaXie/NDATA documentation built on May 22, 2022, 11:44 a.m.