Displaying 1 result from an estimated 1 matches for "pigs3".
Did you mean:
pigs
2008 Jul 23
1
R2WinBUGS problem
...pow(sigma.alpha,2)
sigma.alpha ~ dunif(0,100)
sigma.AB[2,2] <- pow(sigma.beta,2)
sigma.beta ~ dunif(0,100)
sigma.AB[1,2] <- rho*sigma.alpha*sigma.beta
sigma.AB[2,1] <- sigma.AB[1,2]
rho ~ dunif(-1,1)
}",
file="ml3_pig.bug")
pigs3.inits <- function () {
list (AB=array(c(runif(16,0,10),runif(16,-30,0)),c(16,2)),
mu.alpha=runif(1,0,10), mu.beta=runif(1,-30,0),
sigma.dep=runif(1), sigma.alpha=runif(1), sigma.beta=runif(1),
rho=runif(1))
}
pigs.parameters <- c ("mu.alpha", "mu.beta", "sigm...