icm_randsearch: Estimate the parameters iteratively

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

Estimate the parameters iteratively

Description

Estimate the parameters iteratively

Usage

icm_randsearch(
  nodesNeighbor,
  w = NULL,
  W = NULL,
  Y,
  mu,
  N,
  theta_0_init,
  theta_1_init,
  maxiter = 100,
  threshold = 0.001,
  lambda = 0.5,
  initLabel = NULL
)

Arguments

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.
mu	The vector of mutability for all genes.
N	The parameter in the prior for network.
theta_0_init	The initial value of theta_0.
theta_1_init	The initial value of theta_1.
maxiter	The maximum number of iteration.
threshold	The threshold for convergence.
lambda	The tuning parameter for estimation of theta_0.
initLabel	The vector of initial labels for all genes.

Value

The list containing the estimated values for theta_0_est,theta_1_est and the pseudo-likelihood pseudollk.

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