simu1c work in progress,  Purpuse:
# 1-Test the stability of hertability estimate by decreasing the number of subjects until reaching
# critical point where the script fail to retrive heritability from the simulated data.
# 2-define the critcal mz/dz ration for different sampe size