Почему при использовании anova_lm в Python фигуры не выровнены? - PullRequest
Новая
       72

Почему при использовании anova_lm в Python фигуры не выровнены?

0 голосов
/ 03 июня 2018

Исходя из этого вопроса введите здесь описание ссылки , я использовал statsmodels для реализации ANOVA в python.Мои данные находятся в Pandas DataFrame, а country является категориальной переменной.

def anova(data):

    mod = ols('C(country) ~ playerRank+playerGames', data=data).fit()
    aov_table = sm.stats.anova_lm(mod, typ=2)
    print aov_table

Когда я использую вышеупомянутую функцию, она показывает: 

 File "<ipython-input-32-e77ae8a55692>", line 1, in <module>
    aov_table = sm.stats.anova_lm(mod, typ=2)
  File "C:\ProgramData\Anaconda2\lib\site-packages\statsmodels\stats\anova.py", line 326, in anova_lm
    return anova_single(model, **kwargs)
  File "C:\ProgramData\Anaconda2\lib\site-packages\statsmodels\stats\anova.py", line 83, in anova_single
    robust)
  File "C:\ProgramData\Anaconda2\lib\site-packages\statsmodels\stats\anova.py", line 178, in anova2_lm_single
    cov = _get_covariance(model, None)
  File "C:\ProgramData\Anaconda2\lib\site-packages\statsmodels\stats\anova.py", line 15, in _get_covariance
    return model.cov_params()
  File "C:\ProgramData\Anaconda2\lib\site-packages\statsmodels\base\wrapper.py", line 95, in wrapper
    obj = data.wrap_output(func(results, *args, **kwargs), how)
  File "C:\ProgramData\Anaconda2\lib\site-packages\statsmodels\base\model.py", line 1180, in cov_params
    scale = self.scale
  File "C:\ProgramData\Anaconda2\lib\site-packages\statsmodels\tools\decorators.py", line 97, in __get__
    _cachedval = self.fget(obj)
  File "C:\ProgramData\Anaconda2\lib\site-packages\statsmodels\regression\linear_model.py", line 1510, in scale
    return np.dot(wresid, wresid) / self.df_resid
ValueError: shapes (10481,41) and (10481,41) not aligned: 41 (dim 1) != 10481 (dim 0)

Редактировать:

Я пытался это appraoch как сильфон: 

y, X = dmatrices('CountryCode ~ playerRank+playerGames', data=data[[['playerRank', 'playerGames', 'CountryCode']]], return_type='dataframe')
mod = sm.OLS(y, X)
res = mod.fit()
res.summary()

aov_table = sm.stats.anova_lm(res, typ=2)
print aov_table

, и это показывает следующую ошибку: 

    Traceback (most recent call last):
  File "C:\ProgramData\Anaconda2\lib\site-packages\IPython\core\interactiveshell.py", line 2881, in run_code
    exec(code_obj, self.user_global_ns, self.user_ns)
  File "<ipython-input-13-e9f27287c5da>", line 1, in <module>
    aov_table = sm.stats.anova_lm(res, typ=2)
  File "C:\ProgramData\Anaconda2\lib\site-packages\statsmodels\stats\anova.py", line 326, in anova_lm
    return anova_single(model, **kwargs)
  File "C:\ProgramData\Anaconda2\lib\site-packages\statsmodels\stats\anova.py", line 68, in anova_single
    design_info = model.model.data.design_info
AttributeError: 'PandasData' object has no attribute 'design_info'
...