statsample-timeseries 0.0.3 → 0.3.0

Sign up to get free protection for your applications and to get access to all the features.
Files changed (36) hide show
  1. checksums.yaml +7 -0
  2. data/.gitignore +19 -0
  3. data/.travis.yml +13 -10
  4. data/Gemfile +2 -21
  5. data/History.md +4 -0
  6. data/LICENSE.txt +1 -1
  7. data/README.md +62 -0
  8. data/Rakefile +12 -17
  9. data/lib/statsample-timeseries.rb +3 -13
  10. data/lib/statsample-timeseries/arima.rb +72 -74
  11. data/lib/statsample-timeseries/arima/kalman.rb +20 -40
  12. data/lib/statsample-timeseries/arima/likelihood.rb +3 -4
  13. data/lib/statsample-timeseries/daru_monkeys.rb +78 -0
  14. data/lib/statsample-timeseries/timeseries/pacf.rb +47 -38
  15. data/lib/statsample-timeseries/utility.rb +105 -133
  16. data/lib/statsample-timeseries/version.rb +5 -0
  17. data/statsample-timeseries.gemspec +31 -0
  18. data/test/helper.rb +6 -29
  19. data/test/test_acf.rb +41 -0
  20. data/test/test_arima_ks.rb +28 -12
  21. data/test/test_arima_simulators.rb +59 -42
  22. data/test/test_matrix.rb +1 -1
  23. data/test/test_pacf.rb +7 -2
  24. data/test/test_wald.rb +7 -3
  25. metadata +81 -132
  26. data/README.rdoc +0 -72
  27. data/VERSION +0 -1
  28. data/bin/bio-statsample-timeseries +0 -74
  29. data/features/acf.feature +0 -31
  30. data/features/pacf.feature +0 -42
  31. data/features/step_definitions/bio-statsample-timeseries_steps.rb +0 -0
  32. data/features/step_definitions/step_definitions.rb +0 -37
  33. data/features/step_definitions/step_definitions_acf.rb +0 -8
  34. data/features/support/env.rb +0 -15
  35. data/lib/statsample-timeseries/timeseries.rb +0 -291
  36. data/test/test_tseries.rb +0 -103
data/test/test_tseries.rb DELETED
@@ -1,103 +0,0 @@
1
- require(File.expand_path(File.dirname(__FILE__)+'/helper.rb'))
2
-
3
- class StatsampleTestTimeSeries < MiniTest::Unit::TestCase
4
- include Statsample::Shorthand
5
-
6
- # All calculations are compared to the output of the equivalent function in R
7
-
8
- def setup
9
- # daily closes of iShares XIU on the TSX
10
- @xiu = Statsample::TimeSeries::Series.new [17.28, 17.45, 17.84, 17.74, 17.82, 17.85, 17.36, 17.3, 17.56, 17.49, 17.46, 17.4, 17.03, 17.01,
11
- 16.86, 16.86, 16.56, 16.36, 16.66, 16.77], :scale
12
- end
13
-
14
- def test_acf
15
- acf = @xiu.acf
16
-
17
- assert_equal 14, acf.length
18
-
19
- # test the first few autocorrelations
20
- assert_in_delta 1.0, acf[0], 0.0001
21
- assert_in_delta 0.852, acf[1], 0.001
22
- assert_in_delta 0.669, acf[2], 0.001
23
- assert_in_delta 0.486, acf[3], 0.001
24
- end
25
-
26
- def test_lag
27
- #test of default lag
28
- lag1 = @xiu.lag
29
-
30
- assert_in_delta 16.66, lag1[lag1.size - 1], 0.001
31
- assert_in_delta 16.36, lag1[lag1.size - 2], 0.001
32
-
33
- #test with different lagging unit
34
- lag2 = @xiu.lag(2)
35
-
36
- assert_in_delta 16.36, lag2[lag2.size - 1], 0.001
37
- assert_in_delta 16.56, lag2[lag2.size - 2], 0.001
38
- end
39
-
40
- def test_delta
41
- diff = @xiu.diff
42
-
43
- assert_in_delta 0.11, diff[@xiu.size - 1], 0.001
44
- assert_in_delta 0.30, diff[@xiu.size - 2], 0.001
45
- assert_in_delta -0.20, diff[@xiu.size - 3], 0.001
46
- end
47
-
48
- def test_ma
49
- # test default
50
- ma10 = @xiu.ma
51
-
52
- assert_in_delta ma10[-1], 16.897, 0.001
53
- assert_in_delta ma10[-5], 17.233, 0.001
54
- assert_in_delta ma10[-10], 17.587, 0.001
55
-
56
- # test with a different lookback period
57
- ma5 = @xiu.ma 5
58
-
59
- assert_in_delta ma5[-1], 16.642, 0.001
60
- assert_in_delta ma5[-10], 17.434, 0.001
61
- assert_in_delta ma5[-15], 17.74, 0.001
62
- end
63
-
64
- def test_ema
65
- # test default
66
- ema10 = @xiu.ema
67
-
68
- assert_in_delta ema10[-1], 16.87187, 0.00001
69
- assert_in_delta ema10[-5], 17.19187, 0.00001
70
- assert_in_delta ema10[-10], 17.54918, 0.00001
71
-
72
- # test with a different lookback period
73
- ema5 = @xiu.ema 5
74
-
75
- assert_in_delta ema5[-1], 16.71299, 0.0001
76
- assert_in_delta ema5[-10], 17.49079, 0.0001
77
- assert_in_delta ema5[-15], 17.70067, 0.0001
78
-
79
- # test with a different smoother
80
- ema_w = @xiu.ema 10, true
81
-
82
- assert_in_delta ema_w[-1], 17.08044, 0.00001
83
- assert_in_delta ema_w[-5], 17.33219, 0.00001
84
- assert_in_delta ema_w[-10], 17.55810, 0.00001
85
- end
86
-
87
- def test_macd
88
- # MACD uses a lot more data than the other ones, so we need a bigger vector
89
- data = File.readlines(File.dirname(__FILE__) + "/fixtures/stock_data.csv").map(&:to_f).to_time_series
90
-
91
- macd, signal = data.macd
92
-
93
- # check the MACD
94
- assert_in_delta 3.12e-4, macd[-1], 1e-6
95
- assert_in_delta -1.07e-2, macd[-10], 1e-4
96
- assert_in_delta -5.65e-3, macd[-20], 1e-5
97
-
98
- # check the signal
99
- assert_in_delta -0.00628, signal[-1], 1e-5
100
- assert_in_delta -0.00971, signal[-10], 1e-5
101
- assert_in_delta -0.00338, signal[-20], 1e-5
102
- end
103
- end