web-dev-qa-db-ja.com

時系列、月次データを四半期ごとに変更

今、私は次のようないくつかの月次データを持っています：

1/1/90 620
2/1/90,591
3/1/90,574
4/1/90,542
5/1/90,534
6/1/90,545
#...etc

Ts（）関数を使用すると、データを次のような時系列構造にするのは簡単です。

         Jan Feb Mar ... Nov Dec
   1990  620 591 574 ... 493 464
   1991  100 200 300 ...........

次のように四半期ごとに繰り返すように変更する可能性はありますか？

            1st 2nd 3rd 4th
   1990-Q1  620 591 574 464
   1990-Q2  100 200 300 400
   1990-Q3  ...
   1990-Q4  ...
   1991-Q1  ...

変えてみた

   ts(mydata,start=c(1990,1),frequency=12)

に

   ts(mydata,start=c(as.yearqrt("1990-1",1)),frequency=4)

しかし、それは機能していないようです。

誰か助けてもらえますか？どうもありがとうございました。

rdatetime-series
9
user2152814
 monthly <- ts(mydata, start = c(1990, 1), frequency = 12)
 quarterly <- aggregate(monthly, nfrequency = 4)
9
Rob Hyndman

私はこれについてHyndmanに同意しません。 Hyndmanは通常間違いを犯すことができないので、これはまれです。しかし、彼の解決策がOPに彼が望むものを与えていないことをあなたに示すことができます。

test<-c(1:100)
test_ts <- ts(test, start=c(2000,1), frequency=12)
test_ts

     Jan Feb Mar Apr May Jun Jul Aug Sep Oct Nov Dec
2000   1   2   3   4   5   6   7   8   9  10  11  12
2001  13  14  15  16  17  18  19  20  21  22  23  24
2002  25  26  27  28  29  30  31  32  33  34  35  36
2003  37  38  39  40  41  42  43  44  45  46  47  48
2004  49  50  51  52  53  54  55  56  57  58  59  60
2005  61  62  63  64  65  66  67  68  69  70  71  72
2006  73  74  75  76  77  78  79  80  81  82  83  84
2007  85  86  87  88  89  90  91  92  93  94  95  96
2008  97  98  99 100     

test_agg <-  aggregate(test_ts, nfrequency=4)
test_agg

2000    6   15   24   33
2001   42   51   60   69
2002   78   87   96  105
2003  114  123  132  141
2004  150  159  168  177
2005  186  195  204  213
2006  222  231  240  249
2007  258  267  276  285
2008  294

ええと、待ってください、その第1四半期は、3か月の平均ではなく、その合計です。 （1 + 2 + 3 = 6ですが、平均= 2を表示する必要があります）。したがって、少し変更する必要があります。

test_agg <-  aggregate(test_ts, nfrequency=4)/3
# divisor is (old freq)/(new freq) = 12/4 = 3
     Qtr1 Qtr2 Qtr3 Qtr4
2000    2    5    8   11
2001   14   17   20   23
2002   26   29   32   35
2003   38   41   44   47
2004   50   53   56   59
2005   62   65   68   71
2006   74   77   80   83
2007   86   89   92   95
2008   98

これで、四半期ごとに書き込まれた月次データの平均が表示されます。ここでの秘訣は除数です。毎週（freq = 52）で、四半期（freq = 4）が必要な場合は、52/4 = 13で割ります。

7
NonReformed Bayesianist

合計ではなく平均が必要な場合は、「平均」を追加するだけです。

quarterly <- aggregate(monthly, nfrequency=4,mean)
3
Cheuk