I would code these as lognormally distributed, with the parameters being the natural logs of the mean and standard deviation. As I said above: Something like: lnweight_m=rand('NORMAL', 5.24597, 4.0792); weight_lbs_m=round(exp(lnweight_m),1); where the parameters passed to the rand function are the natural logs of the values you currently use. This might be more useful. But be warned, the resulting data will be skewed. No getting around that. And I mean, there is no getting around the fact that the real data are not normally distributed. If it were, there would be a nonzero probability of people with negative heights or weights. Steve Denham
... View more