Learn R Programming

SDGLM (version 0.4.0)

mutRate: Scalable Mutation-Rate Strategies for Sca-MCMC

Description

Computes the mutation-rate vector Q_P according to three scalable schemes proposed in the reference paper.

Usage

mutRate(type = c("ScaI", "ScaII", "ScaIII"), N, L, beta.star, beta0)

Value

Numeric vector of length 2, r, or N depending on type, summing to Q0 and proportional to the chosen strategy.

Arguments

type

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

N

Integer > 1, number of parallel chains.

L

Integer > 0, length of the parameter vector.

beta.star

Target parameter vector (binary or factor).

beta0

Initial parameter vector (same length as beta.star).

Examples

Run this code
beta.star <- c(1, 0, 1, 1, 0)
beta0 <- c(1, 1, 1, 0, 0)
mutRate("ScaII", N = 8, L = 5, beta.star, beta0)

Run the code above in your browser using DataLab