[求解]请大家帮忙看看这个boot函数程序。

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

第一次用boot函数，还不太懂。现在想要计算如下Diff函数的频数分布。程序总是报错，新手改不好。。。

1. b<-data.frame(a$profit)
   
2. 
   
3. library(boot)  
   
4. beta<-function(data,indices){  
   
5. d<-data[indices,]
   
6. Max<-0;
   
7. Diff<-0;
   
8. SUM<-0;
   
9. Diff=matrix(0,1,839);
   
10. Max=matrix(0,1,840);
    
11. SUM=matrix(0,1,839);
    
12. 
    
13. for(i in 1:838){
    
14.                   SUM[1,i+1]=SUM[1,i]+data[i+1,1]
    
15.                   }
    
16. 
    
17.   for(j in 1:839){
    
18.                    if(SUM[1,j]>Max[1,j])
    
19.                    {Diff[1,j]<-0
    
20.                    Max[1,j+1]<-SUM[1,j]
    
21.                    } else {
    
22.                    Diff[1,j]<-Max[1,j]-SUM[1,j]
    
23.                    Max[1,j+1]<-Max[1,j]
    
24.                    }
    
25.                   }
    
26.             
    
27. return(density(d))  
    
28. }  
    
29. 
    
30. result<-boot(data=b,statistic=beta,R=1000)  
    
31. boot.ci(result);

复制代码

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：boot Diff DIF IFF 程序 function library return matrix

不懂，看看

不知道楼主这个函数的目的是什么。
看代码的话第14行data似乎应该改成d。
最后返回值是density(d)，函数本身用不改变d，这样的话其实第6行到第25行都是没用的。

mymei 发表于 2014-8-4 15:42
不知道楼主这个函数的目的是什么。
看代码的话第14行data似乎应该改成d。
最后返回值是density(d)，函数本 ...

额。复制错了，应该是
b<-data.frame(a$profit)

library(boot)  
beta<-function(data,indices){  
d<-data[indices,]
Max<-0;
Diff<-0;
SUM<-0;
Diff=matrix(0,1,839);
Max=matrix(0,1,840);
SUM=matrix(0,1,839);

for(i in 1:838){
                  SUM[1,i+1]=SUM[1,i]+d[i+1,1]
                  }

  for(j in 1:839){
                   if(SUM[1,j]>Max[1,j])
                   {Diff[1,j]<-0
                   Max[1,j+1]<-SUM[1,j]
                   } else {
                   Diff[1,j]<-Max[1,j]-SUM[1,j]
                   Max[1,j+1]<-Max[1,j]
                   }
                  }
            
return(density(Diff))  
}  

result<-boot(data=b,statistic=beta,R=1000)  
boot.ci(result);

axiaohuohu 发表于 2014-8-5 10:36
额。复制错了，应该是
b

还是不明白你的目的。是否能说得更清楚点？
boot函数是对某个统计量进行bootstrap，但beta函数返回的不是统计量，而是个list，因此会出错。
参看帮助文档：
statistic: A function which when applied to data returns a vector containing the statistic(s) of interest.

mymei 发表于 2014-8-5 19:31
还是不明白你的目的。是否能说得更清楚点？
boot函数是对某个统计量进行bootstrap，但beta函数返回的不是 ...

想要得到定义的Diff数据的频率分布。所以这里beta函数只能返回统计量，不能是频率分布的list是吗？这里出错了是吧？

axiaohuohu 发表于 2014-8-6 10:45
想要得到定义的Diff数据的频率分布。所以这里beta函数只能返回统计量，不能是频率分布的list是吗 ...

对的。