请问如下代码怎么就下标越界了

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

1. z<-rnorm(200,0,1)
   
2. y<-3+z
   
3. x<-c(1:200)
   
4. fit<-lm(y~x)
   
5. plot(x,y,type="l")
   
6. abline(fit,col=2)
   
7. zwave<-residuals(fit)
   
8. 
   
9. 
   
10. S0t<-matrix()
    
11. S1t<-matrix()
    
12. SS0<-matrix()
    
13. SS1<-matrix()
    
14. SSS0<-c()
    
15. SSS1<-c()
    
16. theata<-c()
    
17. for (i in 40:160) {
    
18.   for (i1 in 1:(200-i))
    
19.   {
    
20.    
    
21.     if(i1 == 1)
    
22.       {S1t[1,i]<-zwave[i+1]}
    
23.     else
    
24.       {S1t[i1,i]<-S1t[(i1-1),i]+zwave[i+i1]}
    
25.    
    
26.    
    
27.     SSS1[i]<-sum(S1t[i1,i]^2)
    
28.    }
    
29.   for(i2 in 1:i)
    
30.   {
    
31.    
    
32.     if(i2 == 1){S0t[1,i]<-zwave[1]}
    
33.     else {S0t[i2,i]<-S0t[(i2-1),i]+zwave[i2]}
    
34.    
    
35.   
    
36.     SSS0[i]<-sum(S0t[i2,i]^2)
    
37.   }
    
38.   
    
39. 
    
40.   theata[i]<-((1-i/200)*200)^(-2)*sum(SSS1)/(i^(-2)*sum(SSS0))
    
41. }

复制代码

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：下标越界