I have to create charts (similar to growth charts) for children of ages 5 to 15 years (only 5,6,7 etc; there are no fractional values like 2.6 years) for a health variable which is non-negative, continuous and in the range of 50-150 (with only a few values outside this range). I have to create 90th, 95th and 99th percentile curves and also create tables for these percentiles. The sample size is about 8000.

I checked and found following possible ways:

1. Find quantiles and then use loess method to get a smooth curve from these quantiles. The degree of smoothness can be adjusted by 'span' parameter.

2. Use LMS (Lambda-Mu-Sigma) method (e.g. using gamlss or VGAM packages in R).

3. Use quantile regression.

What is the best way to do it? By 'best' I mean either the ideal method which is the standard method for creation of such growth curves and would be acceptable to all. Or an easier and simpler to implement method, which may have some limitations, but is an acceptable, quicker method. (For example using loess on percentile values is much faster than using LMS of gamlss package).

Thanks for your help.