日ごとの時系列で季節性をチェックしたい。問題は、月の長さ (または頻度) が等しくないことです。31 日、28 日、30 日という月があります。オブジェクトを宣言するときts、固定の頻度しか指定できないため、正しくありません。
> x <- data.frame(d = as.Date("2013-01-01") + 1:365 , v = runif(365))
> tapply(as.numeric(format(x$d,"%d")) , format(x$d,"%m") , max)
01 02 03 04 05 06 07 08 09 10 11 12
31 28 31 30 31 30 31 31 30 31 30 31
後で分解して季節性を確認できる時系列オブジェクトをrで作成するにはどうすればよいですか?
ピボット テーブルを作成して ts に変換することはできますか?