Statsmodels is a Python package that provides a complement to scipy for statistical computations including descriptive statistics and estimation of statistical models.

regression: Generalized least squares (including weighted least squares and least squares with autoregressive errors), ordinary least squares.

glm: Generalized linear models with support for all of the one-parameter exponential family distributions.

discrete choice models: Poisson, probit, logit, multinomial logit

rlm: Robust linear models with support for several M-estimators.

tsa: Time series analysis models, including ARMA, AR, VAR

nonparametric : (Univariate) kernel density estimators

datasets: Datasets to be distributed and used for examples and in testing.

PyDTA: Tools for reading Stata*.dta files into numpy arrays.

System Message: WARNING/2 (string, line 20);backlink

Inline emphasis start-string without end-string.

stats: a wide range of statistical tests

sandbox : There is also a sandbox which contains code for generalized additive models (untested), mixed effects models, cox proportional hazards model (both are untested and still dependent on the nipy formula framework), generating descriptive statistics, and printing table output to ascii, latex, and html. There is also experimental code for systems of equations regression, time series models, panel data estimators and information theoretic measures. None of this code is considered production ready.

Development branches will be on Github. This is where to go to get the most up to date code in the trunk branch. Experimental code is hosted here in branches and in developer forks. This code is merged to master often. We try to make sure that the master branch is always stable.

Source download of stable tags will be on SourceForge.

System Message: WARNING/2 (string, line 55)

The official documentation is hosted on SourceForge.

The sphinx docs are currently undergoing a lot of work. They are not yet comprehensive, but should get you started.

Our blog will continue to be updated as we make progress on the code.

The source distribution for Windows includes a htmlhelp file (statsmodels.chm). This can be opened from the python interpreter

All chatter will take place on the or scipy-user mailing list. We are very interested in receiving feedback about usability, suggestions for improvements, and bug reports via the mailing list or the bug tracker at

to discuss development and design issues that are deemed to be too specialized for the scipy-dev/user list.

scikits.statsmodels has been ported and tested for Python 3.2. Python 3 version of the code can be obtained by running over the entire statsmodels source. The numerical core of statsmodels worked almost without changes, however there can be problems with data input and plotting. The STATA file reader and writer in iolib.foreign has not been ported yet. And there are still some problems with the matplotlib version for Python 3 that was used in testing. Running the test suite with Python 3.2 shows some errors related to foreign and matplotlib.

Changes that break backwards compatibility

Added for importing. So the new convention for importing is

System Message: WARNING/2 (string, line 140)

import scikits.statsmodels.api as sm

Importing from modules directly now avoids unnecessary imports and increases the import speed if a library or user only needs specific functions.

sandbox/ – iolib/

lib/ – iolib/ (Now contains Stata .dta format reader)

Datasets Load class is now load function. – regression/ – discrete/ – robust/ – genmod/

t() method – tvalues attribute (t() still exists but raises a warning)

Time Series Analysis model (tsa) – Vector Autoregression Models VAR (tsa.VAR) – Autogressive Models AR (tsa.AR) – Autoregressive Moving Average Models ARMA (tsa.ARMA) :

System Message: ERROR/3 (string, line 164)

optionally uses Cython for Kalman Filtering use install with option –with-cython

System Message: WARNING/2 (string, line 166)

Block quote ends without a blank line; unexpected unindent.

Baxter-King band-pass filter (tsa.filters.bkfilter)

Hodrick-Prescott filter (tsa.filters.hpfilter)

Christiano-Fitzgerald filter (tsa.filters.cffilter)

Improved maximum likelihood framework uses all available scipy.optimize solvers

Refactor of the datasets sub-package.

Removed RPy dependency for running the test suite.

Refactored codebase/directory structure.

Support for offset and exposure in GLM.

Removed data_weights argument to for Binomial models.

New statistical tests, especially diagnostic and specification tests

General Method of Moment framework in sandbox

renames for more consistency RLM.fitted_values – RLM.fittedvalues GLMResults.resid_dev – GLMResults.resid_deviance

GLMResults, RegressionResults: lazy calculations, convert attributes to properties with _cache

expanded examples in examples directory

add PyDTA to — functions for reading Stata*.dta binary files and converting them to numpy arrays

System Message: WARNING/2 (string, line 197);backlink

Inline emphasis start-string without end-string.

made tools.categorical much more robust

add_constant now takes a prepend argument

fix GLS to work with only a one column design

add four new Maximum Likelihood Estimators for models with a discrete dependent variables with examples

System Message: ERROR/3 (string, line 213)

add sandbox.tsa (time series analysis) and sandbox.regression (anova)

add principal component analysis in

add Seemingly Unrelated Regression (SUR) and Two-Stage Least Squares for systems of equations in sandbox.sysreg.Sem2SLS

Builds marked with a lock icon are only available via PyPM to users with a currentActivePython Business Editionsubscription.

ActivePython Enterprise Editionguarantees priority access to technical support, indemnification, expert consulting and quality-assured language builds.

Get re-distribution rights and eliminate legal risks withActivePython OEM Edition.