Package: algebraic.mle 2.0.3

Alexander Towell

algebraic.mle: Algebraic Maximum Likelihood Estimators

The maximum likelihood estimator (MLE) is a technology: under regularity conditions, any MLE is asymptotically normal with variance given by the inverse Fisher information. This package exploits that structure by defining an algebra over MLEs. Compose independent estimators into joint MLEs via block-diagonal covariance ('joint'), optimally combine repeated estimates via inverse-variance weighting ('combine'), propagate transformations via the delta method ('rmap'), and bridge to distribution algebra via conversion to normal or multivariate normal objects ('as_dist'). Supports asymptotic ('mle', 'mle_numerical') and bootstrap ('mle_boot') estimators with a unified interface for inference: confidence intervals, standard errors, AIC, Fisher information, and predictive intervals. For background on maximum likelihood estimation, see Casella and Berger (2002, ISBN:978-0534243128). For the delta method and variance estimation, see Lehmann and Casella (1998, ISBN:978-0387985022).

Authors:Alexander Towell [aut, cre]

algebraic.mle_2.0.3.tar.gz
algebraic.mle_2.0.3.zip(r-4.7-any)algebraic.mle_2.0.3.zip(r-4.6-any)
algebraic.mle_2.0.3.tar.gz(r-4.7-any)algebraic.mle_2.0.3.tar.gz(r-4.6-any)
algebraic.mle_2.0.3.tgz(r-4.6-emscripten)
manual.pdf |manual.html✨
DESCRIPTION |NEWS
card.svg |card.png
algebraic.mle/json (API)

# Install 'algebraic.mle' in R:
install.packages('algebraic.mle', repos = c('https://cran.r-universe.dev', 'https://cloud.r-project.org'))

Bug tracker:https://github.com/queelius/algebraic.mle/issues

Pkgdown/docs site:https://queelius.github.io

On CRAN:

Conda:

5.46 score 8 packages 99 scripts 448 downloads 27 exports 6 dependencies

Last updated from:30b5b64b3b. Checks:6 OK. Indexed: no.

TargetResultTimeFilesSyslog
linux-develOK135
source / vignettesOK283
linux-releaseOK166
windows-develOK94
windows-releaseOK80
wasm-releaseOK98

Exports:as_distbiascdfcombineconditionalconfint_from_sigmaexpectationinv_cdfis_mleis_mle_bootjointmarginalmlemle_bootmle_numericalmsenparamsobsobserved_fimorthogonalparamspredrmapsamplerscore_valsesup

Dependencies:algebraic.distbootMASSmvtnormnumDerivR6

Dynamic failure rate model
Installation | Purpose | API Overview | Fitting exponential models | Hypothesis test and model selection

Last update: 2026-03-19
Started: 2026-01-09

Fitting Common Distributions to a DGP
Introduction | Data Simulation | Visualizing Data | Parametrically Modeling the Data | Maximum Likelihood Estimation | Performance Measures of the MLE | Invariance Property of the MLE | Delta Method | Monte-Carlo Method | Example 1 | Combining Independent MLEs | Example 2 | Bootstrapping the MLEs | Goodness-of-Fit | Prediction Intervals | Conclusion

Last update: 2026-03-19
Started: 2026-01-09

Statistics and characteristics of the MLE
Installation | Normal distribution | Monte-carlo (MC) simulation of the sampling distribution of the MLE | Bias | Variance-covariance matrix | Confidence intervals | Mean squared error matrix | Bootstrap of the sampling distribution of the MLE | Prediction intervals | Combining Independent MLEs

Last update: 2026-03-19
Started: 2026-01-09

The Algebra of MLEs
Introduction | Creating MLEs | Direct construction with mle() | From numerical optimization with mle_numerical() | From bootstrap with mle_boot() | Composing Independent MLEs | Combining Repeated Estimates | Transformations via Invariance | Univariate example: rate to mean lifetime | Multivariate example: component rates to system reliability | Bridging to Distribution Algebra | Full Pipeline

Last update: 2026-03-19
Started: 2026-03-19

Readme and manuals

Help Manual

Help pageTopics
Convert an MLE to a distribution object.as_dist.mle_fit
Convert a bootstrap MLE to an empirical distribution.as_dist.mle_fit_boot
Generic method for computing the bias of an estimator object.bias
Computes the bias of an `mle_fit` object assuming the large sample approximation is valid and the MLE regularity conditions are satisfied. In this case, the bias is zero (or zero vector).bias.mle_fit
Computes the estimate of the bias of a `mle_fit_boot` object.bias.mle_fit_boot
CDF of the asymptotic distribution of an MLE.cdf.mle_fit
Extract coefficients from an 'mle_fit' object.coef.mle_fit
Combine independent MLEs for the same parameter.combine combine.list combine.mle_fit
Conditional distribution from an MLE.conditional.mle_fit
Function to compute the confidence intervals from a variance-covariance matrixconfint_from_sigma
Function to compute the confidence intervals of `mle_fit` objects.confint.mle_fit
Method for obtained the confidence interval of an `mle_fit_boot` object. Note: This impelements the `vcov` method defined in `stats`.confint.mle_fit_boot
PDF of the asymptotic distribution of an MLE.density.mle_fit
PDF of the empirical distribution of bootstrap replicates.density.mle_fit_boot
Dimension (number of parameters) of an MLE.dim.mle_fit
Dimension (number of parameters) of a bootstrap MLE.dim.mle_fit_boot
Expectation operator applied to `x` of type `mle_fit` with respect to a function `g`. That is, `E(g(x))`.expectation.mle_fit
Quantile function of the asymptotic distribution of an MLE.inv_cdf.mle_fit
Determine if an object `x` is an `mle_fit` object.is_mle
Determine if an object is an `mle_fit_boot` object.is_mle_boot
Compose independent MLEs into a joint MLE.joint joint.mle_fit
Log-likelihood of an 'mle_fit' object.logLik.mle_fit
Method for obtaining the marginal distribution of an MLE that is based on asymptotic assumptions:marginal.mle_fit
Mean of the asymptotic distribution of an MLE.mean.mle_fit
Mean of bootstrap replicates.mean.mle_fit_boot
Constructor for making `mle_fit` objects, which provides a common interface for maximum likelihood estimators.mle
Bootstrapped MLEmle_boot
This function takes the output of `optim`, `newton_raphson`, or `sim_anneal` and turns it into an `mle_fit_numerical` (subclass of `mle_fit`) object.mle_numerical
Generic method for computing the mean squared error (MSE) of an estimator, `mse(x) = E[(x-mu)^2]` where `mu` is the true parameter value.mse
Computes the MSE of an `mle_fit` object.mse.mle_fit
Computes the estimate of the MSE of a `boot` object.mse.mle_fit_boot
Method for obtaining the number of observations in the sample used by an `mle_fit`.nobs.mle_fit
Method for obtaining the number of observations in the sample used by an `mle_fit_boot`.nobs.mle_fit_boot
Method for obtaining the number of parameters of an `mle_fit` object.nparams.mle_fit
Method for obtaining the number of parameters of an `boot` object.nparams.mle_fit_boot
Method for obtaining the observations used by the `mle_fit` object `x`.obs.mle_fit
Method for obtaining the observations used by the `mle_fit_boot`.obs.mle_fit_boot
Generic method for computing the observed FIM of an `mle_fit` object.observed_fim
Function for obtaining the observed FIM of an `mle_fit` object.observed_fim.mle_fit
Generic method for determining the orthogonal parameters of an estimator.orthogonal
Method for determining the orthogonal components of an `mle_fit` object `x`.orthogonal.mle_fit
Method for obtaining the parameters of an `mle_fit` object.params.mle_fit
Method for obtaining the parameters of an `boot` object.params.mle_fit_boot
Generic method for computing the predictive confidence interval given an estimator object `x`.pred
Estimate of predictive interval of `T|data` using Monte Carlo integration.pred.mle_fit
Print method for `mle_fit` objects.print.mle_fit
Function for printing a `summary` object for an `mle_fit` object.print.summary_mle_fit
Computes the distribution of `g(x)` where `x` is an `mle_fit` object.rmap.mle_fit
Method for sampling from an `mle_fit` object.sampler.mle_fit
Method for sampling from an `mle_fit_boot` object.sampler.mle_fit_boot
Generic method for computing the score of an estimator object (gradient of its log-likelihood function evaluated at the MLE).score_val
Computes the score of an `mle_fit` object (score evaluated at the MLE).score_val.mle_fit
Generic method for obtaining the standard errors of an estimator.se
Function for obtaining an estimate of the standard error of the MLE object `x`.se.mle_fit
Function for obtaining a summary of `object`, which is a fitted `mle_fit` object.summary.mle_fit
Support of the asymptotic distribution of an MLE.sup.mle_fit
Computes the variance-covariance matrix of `mle_fit` object.vcov.mle_fit
Computes the variance-covariance matrix of `boot` object. Note: This impelements the `vcov` method defined in `stats`.vcov.mle_fit_boot