
library(lmtest)
getwd()
setwd("C:/UseE2-1rs/花寻觅/Documents/xuexuexue")
a=read.csv("C:/Users/花寻觅/Documents/xuexuexue/上机作业2.csv")
h<-lm(Y~X+T,data=a)
summary(h)
h1<-lm(Y~T,data=a)
E1<-residuals(h1)
h2<-lm(X~T,data=a)
E2<-residuals(h2)
(h3<-lm(E1~0+E2,data=a))#E1~E2-1
summary(h3)#h=0.08645与h3=0.08645
e2<- residuals(h)^2
white<- lm(e2~X*T+I(T^2)+I(X^2)+X+T,data=a)#省略也可以
summary(white)
18*0.6698>11.0705#拒绝H0
bptest(h,~X*T+I(T^2)+I(X^2),data=a)#White检验若P<0.05,拒绝H0
#x1<- a$X
#x2<-a$T
#par(mfrow=c(1,2))
#plot(x1,e2,main="残差图")
#plot(x2,e2,main="残差图")
h4<-lm(log(Y)~log(X)+log(T),data=a)
summary(h4)#0.26832
bptest(h4)


雷达卡



京公网安备 11010802022788号







