Simulate mean per age
simMean(age)
age
the age variable
return predicted means
## Not run: x <- simMean(a) ## End(Not run)