I get the statistic value easily:
statsmodels.stats.stattools.durbin_watson(residuals, axis=0)
Any idea how I can get the p-value? Thank you very much.