statsmodels.graphics.regressionplots.plot_ccpr¶
- statsmodels.graphics.regressionplots.plot_ccpr(results, exog_idx, ax=None)[source]¶
Plot CCPR against one regressor.
Generates a component and component-plus-residual (CCPR) plot.
- Parameters:
- results
result
instance
A regression results instance.
- exog_idx{
int
,str
} Exogenous, explanatory variable. If string is given, it should be the variable name that you want to use, and you can use arbitrary translations as with a formula.
- ax
AxesSubplot
,optional
If given, it is used to plot in instead of a new figure being created.
- results
- Returns:
Figure
If ax is None, the created figure. Otherwise the figure to which ax is connected.
See also
plot_ccpr_grid
Creates CCPR plot for multiple regressors in a plot grid.
Notes
The CCPR plot provides a way to judge the effect of one regressor on the response variable by taking into account the effects of the other independent variables. The partial residuals plot is defined as Residuals + B_i*X_i versus X_i. The component adds the B_i*X_i versus X_i to show where the fitted line would lie. Care should be taken if X_i is highly correlated with any of the other independent variables. If this is the case, the variance evident in the plot will be an underestimate of the true variance.
References
http://www.itl.nist.gov/div898/software/dataplot/refman1/auxillar/ccpr.htm
Examples
Using the state crime dataset plot the effect of the rate of single households (‘single’) on the murder rate while accounting for high school graduation rate (‘hs_grad’), percentage of people in an urban area, and rate of poverty (‘poverty’).
>>> import statsmodels.api as sm >>> import matplotlib.pyplot as plt >>> import statsmodels.formula.api as smf
>>> crime_data = sm.datasets.statecrime.load_pandas() >>> results = smf.ols('murder ~ hs_grad + urban + poverty + single', ... data=crime_data.data).fit() >>> sm.graphics.plot_ccpr(results, 'single') >>> plt.show()
(Source code, png, hires.png, pdf)