livia wrote: > I would like to fit a Pareto Distribution and I am using the following codes. > > 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 a single value for > all. > > fit=vglm(ycf1 ~ 1, pareto1(location=alpha), trace=TRUE, crit="c") > fitted(fit) > > The result is > fitted(fit) > [,1] > [1,] 0.07752694 > [2,] 0.07752694 > [3,] 0.07752694 > [4,] 0.07752694 > [5,] 0.07752694 > [6,] 0.07752694 > [7,] 0.07752694 > [8,] 0.07752694 > [9,] 0.07752694 > [10,] 0.07752694 > [11,] 0.07752694 > [12,] 0.07752694 > [13,] 0.07752694 > > Could anybody give me some advice? >
Your model only includes an intercept, so the fitted value is supposed to be the same for all units (there is nothing in your model that allows the fitted value to vary across units). markus -- Markus Jantti Abo Akademi University [EMAIL PROTECTED] http://www.iki.fi/~mjantti ______________________________________________ [email protected] mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.
