commu=read.table(file="commu.txt",header=T)
y=commu[,2]
x1=commu[,3]
x2=commu[,4]
lm.com=lm(log(y)~x1+x2)
library(bstats)
bptest(lm.com)#这个包自动做BP TEST，p=0.2199
library(car)
ncvTest(lm.com)
u1=lm.com$residual
u2=u1^2
u2
lm.u2=lm(u2~x1+x2+I(x1^2)+I(x2^2))
summary(lm.u2)
#lm.u2的可决系数是0.4089
tr=0.4089*9
tr
dchisq(tr,4)
pchisq(tr,4)
ttt=1-pchisq(tr,4)
ttt#我自己算出来p=0.4510


想请问一下为什么结果差了一倍？我编的代码哪里错了？