時間のベクトルを data.frame の列に入れます。以下のサンプル データには、現在の時刻から 2 年後のランダムな時刻の 1000 のタイムスタンプがあります。
dat = data.frame(timestamp = Sys.time() + sort(round(runif(1000, (24*3600), (2*365*24*3600)))))
次のステップは、タイムスタンプがどの月と年にあるかを識別する新しい列を作成することです。
dat$month = strftime(dat$time, "%b")
dat$year = strftime(dat$time, "%Y")
count
これで、plyr
パッケージを使用して、各年の月ごとのタイムスタンプを数えることができます。
library(plyr)
timestamps_month = count(dat, vars = c("month","year"))
ggplot2 でヒストグラムを作成します。
library(ggplot2)
ggplot(data = timestamps_month) + geom_bar(aes(x = month, y = freq, fill = year), stat="identity", position = "dodge")
結果のプロットがどのように見えるかの例については、この SO 投稿を参照してください。
CSV時間データを使用してRでヒストグラムを作成するには?