[R] Fitted Value Pareto Distribution

Markus Jäntti markus.jantti at iki.fi
Wed Jun 13 16:36:39 CEST 2007


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



More information about the R-help mailing list