The Breusch-Pagan (BP) test is one of the most
common tests for heteroskedasticity. It begins by allowing the
heteroskedasticity process to be a function of one or more of
your independent variables, and it’s usually applied by
assuming that heteroskedasticity may be a linear function of
all the independent variables in the model. This assumption
can be expressed as