I would like to fit a Pareto Distribution and I am using the following codes. First, I thought the fitted (fit1) should be the fitted value for the data, is it correct? As the result of the "fitted" turns out to be the same value. fit=vglm(ycf1 ~ 1, pareto1(location=alpha), trace=TRUE, crit="c") coef(fit, matrix=TRUE) summary(fit) fitted(fit) Secondly, how can I plot the density for the fitted distribution? Many thanks. -- View this message in context: http://www.nabble.com/Pareto-Distribution-tf3908751.html#a11082669 Sent from the R help mailing list archive at Nabble.com.