Felix Distribution Family Function
Estimates the parameter of a Felix distribution by maximum likelihood estimation.
felix(lrate = extlogitlink(min = 0, max = 0.5), imethod = 1)
lrate |
Link function for the parameter,
called a below;
see |
imethod |
See |
The Felix distribution is an important basic Lagrangian distribution. The density function is
f(y;a) = (1 / ((y-1)/2)!) * y^((y-3)/2) * a^((y-1)/2) * exp(-ay)
where y=1,3,5,… and 0 < a < 0.5. The mean is 1/(1-2a) (returned as the fitted values). Fisher scoring is implemented.
An object of class "vglmff" (see vglmff-class).
The object is used by modelling functions such as vglm
and vgam.
T. W. Yee
Consul, P. C. and Famoye, F. (2006). Lagrangian Probability Distributions, Boston, USA: Birkhauser.
fdata <- data.frame(y = 2 * rpois(n = 200, 1) + 1) # Not real data! fit <- vglm(y ~ 1, felix, data = fdata, trace = TRUE, crit = "coef") coef(fit, matrix = TRUE) Coef(fit) summary(fit)
Please choose more modern alternatives, such as Google Chrome or Mozilla Firefox.