Skip to main content

Generalized M-Estimation

Project description

delicatessen

Delicatessen

tests version arXiv docs Downloads

The one-stop sandwich (variance) shop in Python. delicatessen is a Python 3.8+ library for the generalized calculus of M-estimation.

Citation: Zivich PN, Klose M, Cole SR, Edwards JK, & Shook-Sa BE. (2022). Delicatessen: M-Estimation in Python. arXiv:2203.11300 [stat.ME]

M-Estimation and Estimating Equations

Here, we provide a brief overview of M-estimation theory. For more detailed introductions to M-estimation, see Ross et al. (2024), Stefanski & Boos (2002), or Chapter 7 of Boos & Stefanski (2013). M-estimation is a generalization of likelihood-based methods. M-estimators are solutions to estimating equations. To apply the M-estimator, we solve the estimating equations using observed data. This is similar to other approaches, but the key advantage of M-Estimators is variance estimation via the empirical sandwich variance estimator.

While M-Estimation is a powerful tool, the derivatives and matrix algebra can quickly become unwieldy. This is where delicatessen comes in. delicatessen takes estimating functions and data, and solves for the parameter estimates, computes the derivatives, and performs the matrix algebra calculations. Therefore, M-estimators can be more easily adopted without having to perform by-hand calculations. In other words, we can let the computer do the math for us.

To further ease use, delicatessen also comes with a variety of built-in estimating equations. See the delicatessen website for details on the available estimating equations, how to use them, and practical examples.

Installation

Installing:

You can install via python -m pip install delicatessen

Dependencies:

The dependencies are: numpy, scipy

To replicate the tests located in tests/, you will additionally need to install: panda, statsmodels, and pytest

While versions of delicatessen prior to v1.0 were compatible with older versions of Python 3 and NumPy and SciPy, the v1.0+ releases are only available for Python 3.8+ with NumPy v1.18.5+ and SciPy v1.9.0. This change was made to use a better numerical approximation procedure for the derivative. If you want to use with older versions of those packages or older versions of Python, install v0.6 instead.

Getting started

Below is a simple demonstration of calculating the mean with delicatessen

import numpy as np
y = np.array([1, 2, 3, 1, 4, 1, 3, -2, 0, 2])

Loading the M-estimator functionality, building the corresponding estimating equation for the mean, and printing the results to the console

from delicatessen import MEstimator

def psi(theta):
    return y - theta[0]

estr = MEstimator(psi, init=[0, ])
estr.estimate()

print(estr.theta)     # Estimate of the mean
print(estr.variance)  # Variance estimate

For further details on using delicatessen, see the full documentation and worked examples available at delicatessen website.

References

Boos DD, & Stefanski LA. (2013). M-estimation (estimating equations). In Essential Statistical Inference (pp. 297-337). Springer, New York, NY.

Stefanski LA, & Boos DD. (2002). The calculus of M-estimation. The American Statistician, 56(1), 29-38.

Ross RK, Zivich PN, Stringer JS, & Cole SR. (2024). M-estimation for common epidemiological measures: introduction and applied examples. International Journal of Epidemiology, 53(2).

Zivich PN, Klose M, Cole SR, Edwards JK, & Shook-Sa BE. (2022). Delicatessen: M-Estimation in Python. arXiv preprint arXiv:2203.11300.

Project details


Download files

Download the file for your platform. If you're not sure which to choose, learn more about installing packages.

Source Distribution

delicatessen-3.1.tar.gz (92.3 kB view details)

Uploaded Source

Built Distribution

delicatessen-3.1-py2.py3-none-any.whl (132.9 kB view details)

Uploaded Python 2Python 3

File details

Details for the file delicatessen-3.1.tar.gz.

File metadata

  • Download URL: delicatessen-3.1.tar.gz
  • Upload date:
  • Size: 92.3 kB
  • Tags: Source
  • Uploaded using Trusted Publishing? No
  • Uploaded via: twine/4.0.2 CPython/3.9.4

File hashes

Hashes for delicatessen-3.1.tar.gz
Algorithm Hash digest
SHA256 5b5ee76e8813e4af81a251a3daa4edabf39070ad9946ff85bc6d761031961040
MD5 b640bb5d8b78a6d95ebc652fbbbb5167
BLAKE2b-256 8b21533eada096f493290b358052db43f199d282cfc6c41c44e9c3e218add54c

See more details on using hashes here.

File details

Details for the file delicatessen-3.1-py2.py3-none-any.whl.

File metadata

  • Download URL: delicatessen-3.1-py2.py3-none-any.whl
  • Upload date:
  • Size: 132.9 kB
  • Tags: Python 2, Python 3
  • Uploaded using Trusted Publishing? No
  • Uploaded via: twine/4.0.2 CPython/3.9.4

File hashes

Hashes for delicatessen-3.1-py2.py3-none-any.whl
Algorithm Hash digest
SHA256 6ff20a975d6ebec5b4b0e3f2f96b28ed2efd624b672af9155cde4f3fb9f0688a
MD5 79a6cc8f7d01ca4b11ef2d75e2444dca
BLAKE2b-256 09883eee3841d07212797d0b0157ef9038ffc5b91365cfa45436638e6050c0a5

See more details on using hashes here.

Supported by

AWS Cloud computing and Security Sponsor Datadog Monitoring Fastly CDN Google Download Analytics Pingdom Monitoring Sentry Error logging StatusPage Status page