Re: [R] Fitted Value Pareto Distribution

From: Markus Jäntti <markus.jantti_at_iki.fi>
Date: Wed, 13 Jun 2007 17:36:39 +0300

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
markus.jantti_at_iki.fi
http://www.iki.fi/~mjantti

______________________________________________
R-help_at_stat.math.ethz.ch 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.
Received on Wed 13 Jun 2007 - 15:25:33 GMT

Archive maintained by Robert King, hosted by the discipline of statistics at the University of Newcastle, Australia.
Archive generated by hypermail 2.2.0, at Wed 13 Jun 2007 - 15:31:42 GMT.

Mailing list information is available at https://stat.ethz.ch/mailman/listinfo/r-help. Please read the posting guide before posting to the list.