setwd("F:/Movies/R/方匡南R/R初级的讲义资料/Radv/data/section 1")
dd=read.csv(file="dd.csv",header=T)
betd=read.csv(file="beta_daily.csv",header=T)
y=betd[,2]
dd=as.data.frame(dd)
n=seq(1:30)
for (i in 2:31){
n[i-1]=length(which(dd[,i]==0))
}
r1=seq(1:30)
for (i in 3:32){
y1=y[1:n[i-2]]
x=betd[,i]
x1=x[1:n[i-2]]
lm1=lm(y1~x1)
r=as.matrix(as.vector(resid(lm1)))
r1[i-2]=t(r)%*%r
r2=seq(1:30)
for (i in 3:32){
y2=y[n[i-2]+1,length(y)]
x=betd[,i]
x2=x[n[i-2]+1,length(x)]
lm1=lm(y2~x2)
r=as.matrix(as.vector(resid(lm2)))
r2[i-2]=t(r)%*%r
r3=seq(1:30)
for (i in 3:32){
y3=y[1,length(y)]
x=betd[,i]
x3=x[1,length(x)]
lm1=lm(y3~x3)
r=as.matrix(as.vector(resid(lm3)))
r3[i-2]=t(r)%*%r
ff=seq(1:30)
for (i in 1:30){
ff(i)=(r3-(r1+r2))/2/((r1+r2)/300)
}
ff
qf(0.95,2,300)
我在R-STUDIO里面运行后它不输出结果,显示的是:
setwd("F:/Movies/R/方匡南R/R初级的讲义资料/Radv/data/section 1")+ dd=read.csv(file="dd.csv",header=T)+ betd=read.csv(file="beta_daily.csv",header=T)+ y=betd[,2]+ dd=as.data.frame(dd)+ n=seq(1:30)+ for (i in 2:31){+ n[i-1]=length(which(dd[,i]==0))+ }+ + r1=seq(1:30)+ for (i in 3:32){+ y1=y[1:n[i-2]]+ x=betd[,i]+ x1=x[1:n[i-2]]+ lm1=lm(y1~x1)+ r=as.matrix(as.vector(resid(lm1)))+ r1[i-2]=t(r)%*%r+ + r2=seq(1:30)+ for (i in 3:32){+ y2=y[n[i-2]+1,length(y)]+ x=betd[,i]+ x2=x[n[i-2]+1,length(x)]+ lm1=lm(y2~x2)+ r=as.matrix(as.vector(resid(lm2)))+ r2[i-2]=t(r)%*%r+ + r3=seq(1:30)+ for (i in 3:32){+ y3=y[1,length(y)]+ x=betd[,i]+ x3=x[1,length(x)]+ lm1=lm(y3~x3)+ r=as.matrix(as.vector(resid(lm3)))+ r3[i-2]=t(r)%*%r+ + ff=seq(1:30)+ for (i in 1:30){+ ff(i)=(r3-(r1+r2))/2/((r1+r2)/300)+ }+ ff
为什么我的R不输出结果了? Rstudio的BUG吗?