## Code adapted from Ben Bolker

fruitdat1 = read.csv("../DATA/wildflower_pods_complete.csv")

library(MCMCglmm)
st3 <- system.time(nfit2 <- MCMCglmm(toF ~State+Pods,
                                     random=~idh(1):Year+idh(Pods):PlantID,
                                     family="categorical", data=fruitdat1,
                                     nitt = 101000, burnin = 1000, thin = 100))

save(nfit2, file = 'example_MCMCglmm_sample.Rdata')
summary(nfit2)
st3["elapsed"]
xyplot(as.mcmc(nfit2$Sol),layout=c(2,3))
xyplot(as.mcmc(nfit2$VCV),layout=c(2,3))