全部版块 我的主页
论坛 提问 悬赏 求职 新闻 读书 功能一区 经管百科 爱问频道
4061 3
2014-06-18
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吗?
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

全部回复
2014-6-18 00:50:42
天王腹黑 发表于 2014-6-18 00:30
setwd("F:/Movies/R/方匡南R/R初级的讲义资料/Radv/data/section 1")
dd=read.csv(file="dd.csv",header= ...
是邹至庄的“邹氏检验”吗
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

2014-6-18 01:18:11
JiaBinXu 发表于 2014-6-18 00:50
是邹至庄的“邹氏检验”吗
是的!
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

2014-6-18 20:39:07
最后输出结构的字体颜色很奇怪,麻烦调整一下,方便他人查阅
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

相关推荐
栏目导航
热门文章
推荐文章

分享

扫码加好友,拉您进群