Displaying 1 result from an estimated 1 matches for "stateabb".
Did you mean:
statbb
2013 Mar 28
3
problem with plots with short example.
...mcmc(FE.out$mcmc)
exchange.chains <- combine.mcmc(exchange.out$mcmc)
mu.bar <- apply(FE.chains[, grep("mu\\[", colnames(FE.chains))], 2, mean)
mu.bar2 <- apply(hier.chains[, grep("mu\\[", colnames(hier.chains))], 2, mean)
ns <- aggregate(therms$numstate, list(therms$stateabb), length)
plot(mu.bar, mu.bar2, cex=sqrt(ns[,2])/3,
xlab = "FE mu[j]",
ylab = "Hierarchical mu[j]")
abline(a=0, b=1)
###################################################
### code chunk number 7: dotchart
###################################################
fe.mu <- FE.cha...