bptest结果与人工计算结果不同

Question

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

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

如果妇女有孩子，母亲=1。如果年份为1993年或以后，则为1。Y是一个虚拟变量，如果一个女人工作，它等于1。我运行了模型，我得到了结果。在那之后我进行了一次测试-

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 

我得到的F值和BP检验不一样。F值大到足以表示异方差，但与BP值不同。为什么会这样呢？

Answer 1

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