Code from:

Wood, J.L.A, Yates, M.C., and D.J. Fraser. 2016. Are heritability and selection related to population size
in nature? Meta-analysis and conservation implications. Evolutionary Applications.

################################################################################################################

library(MCMCglmm)

############HERITABILITY ANALYSIS
############h2 DATA WITH STANDARD ERRORS (META-ANALYSIS)

h2.mev=h2.data$SE^2

prior = list(R = list(V=diag(12),n=0.002),G = list(G1=list(V = 1,n=0.002),G2 = list(V=diag(3),n=0.002), G3 = list(V=diag(12),n=0.002)))

h2.model = MCMCglmm(h2~ln.har.N.centered + trait.class + analysis.type, 
random=~trt + idh(trait.class):study.ID + idh(trait.class:analysis.type):pop.ID, rcov=~idh(trait.class:analysis.type):units,
mev=h2.mev, data=h2.data, prior=prior, nitt=500000, burnin=100000, thin=50, family="gaussian")

summary(h2.model)
posterior.mode(h2.model$Sol)


############ALL h2 data (UNWEIGHTED ANALYSIS)

prior = list(R = list(V=diag(12),n=0.002),G = list(G1=list(V = 1,n=0.002),G2 = list(V=diag(3),n=0.002), G3 = list(V=diag(12),n=0.002)))

h2.model = MCMCglmm(h2~ln.har.N.centered + trait.class + analysis.type, 
random=~trt + idh(trait.class):study.ID + idh(trait.class:analysis.type):pop.ID, rcov=~idh(trait.class:analysis.type):units,
data=h2.data, prior=prior, nitt=500000, burnin=100000, thin=50, family="gaussian")

summary(h2.model)
posterior.mode(h2.model$Sol)