请问R语言哪些包可以实现分数阶差分呢？

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

最近在看ARFIMA的东西，需要做分数阶差分。
自己写了一个，也不知道对不对……R语言自带的diff()函数不支持分数差分。。
请教下有哪些包带这个功能呢？找来对比下看看写对了没。
顺便附上写的代码，这里x是已经读取好的序列：

1. get_sita <- function(k,d){# 使用函数返回分数差分的系数值。
   
2. return(gamma(k-d) / (gamma(-d)*gamma(k+1)))
   
3. }
   
4. 
   
5. 
   
6. d = 0.5
   
7. n = 170 #最多只有171个系数可以计算，再多就超出边界了。
   
8. sita = array(0,n+1) #这就是分数差分项的系数
   
9. for( i in 0:n){
   
10. sita[i+1] = get_sita(i,d)
    
11. }
    
12. 
    
13. 
    
14. N = length(x)
    
15. y = array(0, N) #存放自定义的分数差分结果
    
16. temp = 0
    
17. 
    
18. for( i in 1:N){
    
19. temp = i
    
20. if( i >= n)
    
21. temp = n#
    
22. for( j in 1:temp )
    
23. y[i] = y[i] + sita[j]*x[i-j+1]
    
24. }

复制代码

扫码加我 拉你入群

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

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

分享0 收藏1 回帖

关键词：R语言 ARFIMA Temp Diff sit

同问？

我想可以用package(fracdiff) 里的函数diffseries(x, d)