I would like to fit a gamma/weibull/pareto distribution to a data set of values to produce a CDF. I know I can do this using the fitdist function but I would also like to force the distribution through a given value at a specified percentile to shape the tail of the data. Is this possible? More information: My data:28.744,385.714,20.595,99.350,31.864,77.713,264.408,21.204,31.937,0.900,18.762,173.276,23.707. I also have an additional value of 500.612 which is at the 95th percentile. I know this is a 'bad' way to fit a graph but we need this constraint to shape the tale of the fit. So I would like to fit a gamma distribution to this dataset but the gamma curve must reach 500.612 at the 95th percentile Thanks!