Displaying 1 result from an estimated 1 matches for "mcmcmd".
Did you mean:
mc_cmd
2001 Aug 30
1
MCMC coding problem
...ng error message:
> error in rep(data,t1): invalid number of copies in "rep".
This happens no matter what values that I give the parameters when I call
the function.
Any suggestions on clearing this up would be appreciated! The function is
printed below.
Many Thanks,
Monnie McGee
mcmcmd <- function(d, w, mean, phi2, varwt, ntimes)
{
# A function to run a multi dimensional MCMC for a bivariate
# normal distribution
# d = dimension of problem
# w = weight assigned to larger of two modes in target distribution
# mu = mean of second mode
# phi2 = multiplier of variance factor 2.3...