求用DFA方法求Hurst指数的R程序，谢谢

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

求用DFA方法或DMA法求Hurst指数的R程序，谢谢

扫码加我 拉你入群

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

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

分享0 收藏1 回帖

关键词：Hurst指数 hurst R程序 DFA Urs

1. ###DFA方法估计hurst
   
2. ##获得分组
   
3. getmatrix<-function(x,s){#x是时间序列,s是子列的长度
   
4. x<-cumsum(x-mean(x))
   
5. y<-c(x[1:(round(length(x)/s)*s)],x[(length(x)-round(length(x)/s)*s+1):length(x)])
   
6. N<-length(y)/(2*s)
   
7. ymatrix<-t(matrix(y,s))#分组2*N行s列的矩阵，每一行为一组
   
8. return(ymatrix)
   
9. }
   
10. 
    
11. ###计算消除趋势序列均值
    
12. del<-function(yi,n,s){#yi是各组子列,n是多项式阶数
    
13. xn<-c(1:s)
    
14. model<-lm(yi ~ poly(xn,n))
    
15. return(mean(yi-predict(model)))
    
16. }
    
17. 
    
18. ###计算均值平方根
    
19. fs<-function(x,ymatrix,s,n){#ymatrix是样本矩阵，s是子列的长度，n是多项式阶数
    
20. jisu<-0
    
21. y<-c(x[1:(round(length(x)/s)*s)],x[(length(x)-round(length(x)/s)*s+1):length(x)])
    
22. for(i in 1:(length(y)/s)){
    
23. jisu<-jisu+(del(ymatrix[i,],n,s))^2
    
24. }
    
25. return(sqrt(jisu/(length(y)/(2*s))))
    
26. }
    
27. 
    
28. ###计算hurst指数
    
29. s<-10
    
30. n<-3
    
31. H<-log(fs(x,getmatrix(x,s),s,n),2)/log(s,2)

复制代码

我用数列x<-c(1/(1:100)),n=3,s=10测试，H=-15.22023

1. ###DFA方法估计hurst
   
2. ##获得分组
   
3. getmatrix<-function(x,s){#x是时间序列,s是子列的长度
   
4. x<-cumsum(x-mean(x))
   
5. y<-c(x[1:(round(length(x)/s)*s)],x[(length(x)-round(length(x)/s)*s+1):length(x)])
   
6. N<-length(y)/(2*s)
   
7. ymatrix<-t(matrix(y,s))#分组2*N行s列的矩阵，每一行为一组
   
8. return(ymatrix)
   
9. }
   
10. 
    
11. ###计算消除趋势序列均值
    
12. del<-function(yi,n,s){#yi是各组子列,n是多项式阶数
    
13. xn<-c(1:s)
    
14. model<-lm(yi ~ poly(xn,n))
    
15. return(mean(yi-predict(model)))
    
16. }
    
17. 
    
18. ###计算均值平方根
    
19. fs<-function(x,ymatrix,s,n){#ymatrix是样本矩阵，s是子列的长度，n是多项式阶数
    
20. jisu<-0
    
21. y<-c(x[1:(round(length(x)/s)*s)],x[(length(x)-round(length(x)/s)*s+1):length(x)])
    
22. for(i in 1:(length(y)/s)){
    
23. jisu<-jisu+(del(ymatrix[i,],n,s))^2
    
24. }
    
25. return(sqrt(jisu/(length(y)/(2*s))))
    
26. }
    
27. 
    
28. ###计算hurst指数
    
29. s<-10
    
30. n<-3
    
31. H<-log(fs(x,getmatrix(x,s),s,n),2)/log(s,2)

复制代码

我用x<-c(1/(1:100)),n=3,s=10测试，最后H=-15.22023

###DFA方法估计hurst
##获得分组
getmatrix<-function(x,s){#x是时间序列,s是子列的长度
x<-cumsum(x-mean(x))
y<-c(x[1:(round(length(x)/s)*s)],x[(length(x)-round(length(x)/s)*s+1):length(x)])
N<-length(y)/(2*s)
ymatrix<-t(matrix(y,s))#分组2*N行s列的矩阵，每一行为一组
return(ymatrix)
}

###计算消除趋势序列均值
del<-function(yi,n,s){#yi是各组子列,n是多项式阶数
xn<-c(1:s)
model<-lm(yi ~ poly(xn,n))
return(mean(yi-predict(model)))
}

###计算均值平方根
fs<-function(x,ymatrix,s,n){#ymatrix是样本矩阵，s是子列的长度，n是多项式阶数
jisu<-0
y<-c(x[1:(round(length(x)/s)*s)],x[(length(x)-round(length(x)/s)*s+1):length(x)])
for(i in 1:(length(y)/s)){
jisu<-jisu+(del(ymatrix[i,],n,s))^2
}
return(sqrt(jisu/(length(y)/(2*s))))
}

###计算hurst指数
s<-10
n<-3
H<-log(fs(x,getmatrix(x,s),s,n),2)/log(s,2)

x<-c(1/(1:100)),n=3,s=10,最后H=-15.22023

乐文义 发表于 2018-2-8 17:25
###DFA方法估计hurst
##获得分组
getmatrix

请问第一行的命令输入后出现，此上下文不允许函数定义，是咋回事