SDGLM (version 0.4.0)

compute_mutation_rate: Compute Scalable Mutation-Rate Vector

Description

Implements the three mutation-rate strategies (ScaI, ScaII, ScaIII) described in the reference paper.

Usage

compute_mutation_rate(
  method = c("ScaI", "ScaII", "ScaIII"),
  beta_star,
  beta_init,
  L,
  N_chain
)

Value

List with components:

Q

Mutation-rate vector (length 2, r, or N_chain).

Q0

Base mutation rate (scalar).

r

Scalar, present only for ScaII - number of active components.

w

Vector, present only for ScaIII - normalized random weights.

Arguments

method

Character, one of "ScaI", "ScaII", "ScaIII".

beta_star

Target parameter vector (binary or factor).

beta_init

Current/initial parameter vector.

L

Integer > 0, length of the parameter vector.

N_chain

Integer > 1, number of parallel chains.

Examples

beta_star <- c(1, 0, 1, 1, 0)
beta_init <- c(1, 1, 1, 0, 0)
compute_mutation_rate("ScaII", beta_star, beta_init, L = 5, N_chain = 8)