用R语言写邹检验，但是它不输出结果，我崩溃了！

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

setwd("F:/Movies/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初级的讲义资料/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吗？

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：输出结果 R语言 Section RStudio length Movies 资料

天王腹黑 发表于 2014-6-18 00:30
setwd("F:/Movies/R/方匡南Ｒ/R初级的讲义资料/Radv/data/section 1")
dd=read.csv(file="dd.csv",header= ...

是邹至庄的“邹氏检验”吗

JiaBinXu 发表于 2014-6-18 00:50
是邹至庄的“邹氏检验”吗

是的！

最后输出结构的字体颜色很奇怪，麻烦调整一下，方便他人查阅