bptest 结果与手动计算结果不同

Question

我有一个模型“model1”，它是一个线性概率模型。它衡量了一名妇女在 1993 年之前和之后工作的机会，同时考虑了她是否有孩子。 所以模型 1 是 -

Y = B0 + B1 mother + B2 after + B3 mother * after + epsilon

如果妇女有孩子，则母亲 = 1。如果年份是 1993 年或之后，则 after = 1。 Y 是一个虚拟变量，如果女性工作则等于 1。我运行了模型并得到了结果。 之后我进行了 bptest -

bptest(model1,studentize=TRUE)

studentized Breusch-Pagan test

data:  model1
BP = 20.764, df = 3, p-value = 0.0001178

但是，当我运行回归时 -

residuals(squared) = B0 + B1 mother + B2 after + B3 after mother 

我没有得到与 BP 测试相同的 F 值。 F 值大到足以表示异方差，但它与 BP 值不同。 这是为什么呢？

Answer 1

你从lmtest::bptest得到的值不是F，而是LM = n * R^2的结果，其中R^2 来自residuals^2 = B0 + B1 mother ...