Simulate sd per age
simSD(age)
age
the age variable
return predicted sd
## Not run: x <- simSD(a) ## End(Not run)