ロバスト線形モデル

Norms にリストされているM推定量をサポートするロバスト線形モデル。

コマンドと引数については、 モジュール リファレンス を参照してください。

# Load modules and data
In [1]: import statsmodels.api as sm

In [2]: data = sm.datasets.stackloss.load()

In [3]: data.exog = sm.add_constant(data.exog)

# Fit model and print summary
In [4]: rlm_model = sm.RLM(data.endog, data.exog, M=sm.robust.norms.HuberT())

In [5]: rlm_results = rlm_model.fit()

In [6]: print(rlm_results.params)
const       -41.026498
AIRFLOW       0.829384
WATERTEMP     0.926066
ACIDCONC     -0.127847
dtype: float64

詳細な例は次の場所にあります :

技術文書

参考文献

  • PJ Huber著、 ‘Robust Statistics’ John Wiley and Sons, Inc., New York. 1981年.

  • PJ Huber著. 1973年, ‘The 1972 Wald Memorial Lectures: Robust Regression: Asymptotics, Conjectures, and Monte Carlo.’ The Annals of Statistics, 1.5, 799-821.

  • R Venables, B Ripley著. ‘Modern Applied Statistics in S’ Springer, New York,

  • C Croux, PJ Rousseeuw著, 'Time-efficient algorithms for two highly robust estimators of scale' Computational statistics. Physica, Heidelberg, 1992年.

モジュールリファレンス

モデルクラス

RLM(endog, exog[, M, missing])

ロバスト線形モデル

モデルの結果

RLMResults(model, params, ...)

RLM 結果を含むクラス

Norms

AndrewWave([a])

M 推定のためのアンドリュー波。

Hampel([a, b, c])

M 推定のための Hampel 関数。

HuberT([t])

M 推定のためのHuber の T。

LeastSquares()

M 推定のための最小二乗 rho とその派生関数。

MQuantileNorm(q, base_norm)

基底ノルムに基づく M-分位数目的関数

RamsayE([a])

M 推定のための Ramsay の Ea。

RobustNorm()

ロバスト回帰に使用されるノルムの親クラス。

TrimmedMean([c])

M 推定のためのトリミング平均関数。

TukeyBiweight([c])

M 推定のための Tukey の二重重み関数。

estimate_location(a, scale[, norm, axis, ...])

self.norm を使用した位置の M 推定量とスケールの現在の推定量。

Scale

Huber([c, tol, maxiter, norm])

場所と規模を同時に推定するためのHuberの提案2。

HuberScale([d, tol, maxiter])

ロバスト線形モデルをフィッティングするための Huber のスケーリング。

mad(a[, c, axis, center])

配列の指定された軸に沿った絶対偏差の中央値

hubers_scale

ロバスト線形モデルをフィッティングするための Huber のスケーリング。

iqr(a[, c, axis])

配列の指定された軸に沿った正規化された四分位範囲

qn_scale(a[, c, axis])

スケールの Qn ロバスト推定量を計算する


最終更新日: 2025年01月28日