xy データの単純なテスト セットで statsmodels のロバスト線形モデルを使用しようとしています。ただし、model.params の戻り値として、単一の値しか取得しません。フィットの勾配と切片を取得するにはどうすればよいですか? 最小限の例（外れ値を適合から除外しようとしているため、rlm）：
import statsmodels.api as sm
a = np.array([1,2,3,4,5,6,7,8,9])
b = 2. * np.array([1,2,9,4,5,6,7,13,9])
model = sm.RLM(b, a, M=sm.robust.norms.HuberT()).fit()
model.params
最後の行は のみを返しますarray([2.])。同じパッケージの ols で同じことを試してみましたが、見返りにインターセプトとスロープが得られました。