powered by
Computes the mutation-rate vector Q_P according to three scalable schemes proposed in the reference paper.
mutRate(type = c("ScaI", "ScaII", "ScaIII"), N, L, beta.star, beta0)
Numeric vector of length 2, r, or N depending on type, summing to Q0 and proportional to the chosen strategy.
Character, one of "ScaI", "ScaII", "ScaIII".
Integer > 1, number of parallel chains.
Integer > 0, length of the parameter vector.
Target parameter vector (binary or factor).
Initial parameter vector (same length as beta.star).
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