##Read and check the data data <- read.table(file.choose(), header=T) summary(data) head(data) # scale variables data$zdamAge<-scale(data$damAge) data$zsireAge<-scale(data$sireAge) data$zMaxAge<-scale(data$MaxAge) # Main analysis mod<-MCMCglmm(LRS~zdamAge/sex+sex/zsireAge+zMaxAge, random=~dam+sire+cohort, family="poisson", data=data)