在RcppRoll中使用負權

是否可以在R包RcppRoll中使用負權？在RcppRoll中使用負權

例如，假設我們有一系列dd，我們希望找到一個系列的第一個區別：

library(RcppRoll) 
set.seed(1) 
dd <- 1:100 + runif(100)

所以我們知道dd[100] - dd[99]是0.79 ...寫的另一種方式，這是的當然，(-1L) * dd[99] + (1L) * dd[100]。

第二種方法的優點是使權重顯而易見。因此，權重爲c(-1L, 1L)的加權總和應返回一系列差異;至少這似乎是正確的數學/邏輯。

diff_dd <- roll_sum(dd, n=2, weights = c(-1,1), align='right')

......我們發現，這是所有NaN ...

tail(diff_dd) 
[1] NaN NaN NaN NaN NaN NaN

這是怎麼回事？

來源

2017-07-15 ricardo

問題是normalize參數默認設置爲TRUE。權重被歸一化，並且因爲使用的權重是c(-1, 1)/0，所以1 + -1 = 0事情變得不合時宜。試試這個，你應該得到你想要的。

> diff_dd <- roll_sum(dd, n=2, weights = c(-1,1), align='right', normalize=FALSE) 
> diff_dd 
[1] 1.1066152 1.2007295 1.3353544 0.2934741 1.6967078 1.0462856 0.7161225 0.9683163 0.4326722 1.1441883 
[11] 0.9705822 1.5104661 0.6970809 1.3857377 0.7278578 1.2199193 1.2742876 0.3881291 1.3974100 1.1572600 
[21] 0.2774373 1.4395312 0.4738813 1.1416656 1.1188934 0.6272762 1.3689976 1.4873029 0.4706582 1.1417311 
[31] 1.1174857 0.8939755 0.6926763 1.6411557 0.8410934 1.1257731 0.3137038 1.6157673 0.6875635 1.4096719 
[41] 0.8261139 1.1358726 0.7701035 0.9766833 1.2596367 0.2339750 1.4538989 1.2550837 0.9604178 0.7848881 
[51] 1.3835899 0.5768876 0.8067002 0.8258818 1.0287871 1.2168055 1.2023626 1.1433708 0.7448251 1.5060457 
[61] 0.3807274 1.1654624 0.8733289 1.3184758 0.6071463 1.2205285 1.2877654 0.3179362 1.7910744 0.4637516 
[71] 1.5003674 0.5072431 0.9870914 1.1425763 1.4158471 0.9721411 0.5256501 1.3873312 1.1832973 0.4740415 
[81] 1.2778552 0.6874797 0.9253578 1.4317350 0.4456051 1.5084290 0.4105707 1.1237966 0.8978159 1.0963250 
[91] 0.8193050 1.5833539 1.2339810 0.9026455 1.0183941 0.6579656 0.9548096 1.4007862 0.7940630

來源

2017-07-15 11:18:13 ekstroem

在RcppRoll中使用負權

回答

相關問題