Description Usage Arguments Details References
View source: R/ggamma_flexfit.R
ggamma_flexfit is used to fit a Generalized Gamma distribution to a strictly positive response variable. The scale parameter may be specified either as a function of covariates or as a constant estimated using the response variable alone. If the scale parameter is specified to be a function of covariates, the canonical log link function is used.
1 2 3 |
formula |
An object of class "formula" (or one that can be coerced to that class): a symbolic description of the model to be fitted. |
data |
An optional data frame, list or environment (or object coercible by as.data.frame to a data frame) containing the variables in the model. If not found in data, the variables are taken from environment(formula), typically the environment from which ggamma_flexfit is called. |
weights |
An optional vector of weights to be used in the fitting process. Should be NULL or a numeric vector. |
subset |
An optional vector specifying a subset of observations to be used in the fitting process. |
ownstart |
An optional list containing starting values for the maximum likelihood estimation procedure. If a model with an intercept has been specified, the list must be of the form ownstart = list(alpha = , phi = , beta0 = , beta1 =, …); if a mode with no intercept has been specified the list must be of the form ownstart = list(alpha = , phi = , beta1 = , …); if the rate parameter is not a function of covariates the list must be of the form ownstart = list(alpha = , mu = , phi = ). It is important that the list have as many elements as there are parameters in the model, and that these be supplied in the order set out above. |
key |
A logical parameter dictating whether a key is produced alongside the model’s output. |
warnings |
A logical parameter dictating whether warnings from the maximum likelihood estimation procedure are produced alongside the model’s output. |
... |
Additional arguments to be passed to the function optim within the maximum likelihood estimation procedure. Useful arguments include the gradient descent algorithm to be used and bounds on parameter values; see the stats package. |
This function uses the same parametrization of the Fisk distribution as is used in Stacy (1962); note that this is differs significantly to the parametrization used in many common R packages. Starting values for the maximum likelihood estimation procedure are obtained using the procedure set out in Achcar et al. (2017). The probability probability density function is used is:
f(y) = [α/Γ(φ)] • μ^αφ y^αφ-1 exp(-(μy)^α)
μ is a scale parameter, while φ and α are shape parameters.
When the argument formula specifies a full model with an intercept, μ takes the following form and is estimated via a two step (least squares and maximum likelihood) procedure:
μ = exp(β0 + β1x1 + …. + βkxk)
When the formula argument specifies a model without an intercept, μ takes the bellow form and is estimated via a two step (least squares and maximum likelihood) procedure. Unless theory suggests that an intercept should not be used, users are advised to use a model with an intercept as the maximum likelihood estimation procedure is more stable.
μ = exp(β1x1 + …. + βkxk)
When a null model is specified (formula = y ~ 0) μ is not estimated as a function of covariates. The starting value for the maximum likelihood estimation procedure is obtained by calling the function ggamma_alpha.
Alberto Achcar, Jorge, Pedro Luiz Ramos, and Edson Zangiacomi Martinez. "Some Computational Aspects To Find Accurate Estimates For The Parameters Of The Generalized Gamma Distribution." Pesquisa Operacional Vol.37(2) (2017): n. pag. Print.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.