import wooldridge as wooimport pandas as pdimport statsmodels.api as smimport statsmodels.formula.api as smfphillips = woo.dataWoo('phillips')T = len(phillips)# define yearly time series beginning in 1948:date_range = pd.date_range(start='1948', periods=T, freq='Y')phillips.index = date_range.year# estimation of both Phillips curve models:yt96 = (phillips['year'] <= 1996)phillips['inf_diff1'] = phillips['inf'].diff()reg_s = smf.ols(formula='Q("inf") ~ unem', data=phillips, subset=yt96)reg_ea = smf.ols(formula='inf_diff1 ~ unem', data=phillips, subset=yt96)results_s = reg_s.fit()results_ea = reg_ea.fit()# DW tests:DW_s = sm.stats.stattools.durbin_watson(results_s.resid)DW_ea = sm.stats.stattools.durbin_watson(results_ea.resid)print(f'DW_s: {DW_s}\n')print(f'DW_ea: {DW_ea}\n')
此处可能存在不合适展示的内容,页面不予展示。您可通过相关编辑功能自查并修改。
如您确认内容无涉及 不当用语 / 纯广告导流 / 暴力 / 低俗色情 / 侵权 / 盗版 / 虚假 / 无价值内容或违法国家有关法律法规的内容,可点击提交进行申诉,我们将尽快为您处理。