fitdist_fr_wei: Fit Frechet-Weibull distribution to non-censored data
In jmw86069/salsa: Single-cell Amalgamation by Latent Semantic Analysis (SALSA)

Description Usage Arguments Details Value See Also Examples

Fit Frechet-Weibull distribution to non-censored data

1
2
3

fitdist_fr_wei(x, param_fr_wei = params_fr_wei(), method = c("mle",
  "mme", "qme", "mge"), optim.method = "Nelder-Mead", gof = "KS",
  log = FALSE, verbose = FALSE, ...)

`x`	numeric vector
`param_fr_wei`	`data.frame` of parameter bounds, output from `params_fr_wei()`, which defines the start, upper, and lower bounds for each of the five parameters.
`method`	fitting method, argument passed to `fitdistrplus::fitdist()`.
`optim.method`	character `method` argument passed to `fitdistrplus::mledist` or corresponding function based upon the `method` argument, but ultimately passed to `stats::optim()`.
`gof`	character string compatible with the relevant downstream function goodness of fit method,, currently the default `"KS"` uses the Kolmogorov-Smirnov distance. See `fitdistrplus::mgedist()` for more details.
`verbose`	logical indicating whether to print verbose output.
`...`	additional arguments are passed to `fitdistrplus::fitdist()`.

This function is a simple wrapper around fitdistrplus::fitdist() for the five-parameter Frechet-Weibull distribution. Most customizations are performed by the corresponding support function, for example when method="mle" the distribution fit function is fitdistrplus::mledist(), which then calls stats::optim(). Arguments are passed to these functions using ... where applicable.

TODO: This function currently expects to use method="mle" which utilizes the upper,lower,start parameters for each of the five parameters. Other methods may have different requirements which are not formally handled. In future we will evaluate and support additional methods, modifying the call to fitdistrplus::fitdist() as needed.

object of class "fitdist" described in fitdistrplus::fitdist().

Other SALSA distribution fit functions: dfr_wei, fitdist_fr, params_fr_wei

library(salsa);
data(oz2_numi_per_cell);
x <- oz2_numi_per_cell$count[oz2_numi_per_cell$count >= 16];
fit_fr_wei <- fitdist_fr_wei(x);
# print coefficients
print(coef(fit_fr_wei));