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

Question

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

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

如果女人有孩子，则母亲= 1。如果年份是1993年或之后，则= 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 < / strong>来自residuals^2 = B0 + B1 mother ...