来シーズンのアニメの本数をヒストグラムにしたい

R

http://www.moonphase.cc/Html/anime.html
からデータを取ってきて、Excelに貼り付けます。そのあと、日付の列だけコピーしてRで加工する。本当はWeb scrapingをRでやりたいんだけど、めんどくさいから入力データは事前に作っておく。確かRのドキュメントでも、入力データはPerlとかで加工しておいてね！と言ってたと思ったし。

tb <- table(as.Date(grep(".",readLines("clipboard"),value=T)))

で終わり。

> tb

2008-09-06 2008-09-07 2008-09-13 2008-09-24 2008-09-26 2008-09-28 2008-10-01 2008-10-02 2008-10-03 2008-10-04 2008-10-05 
         1          1          1          1          1          1          6          5          5          7          6 
2008-10-06 2008-10-07 2008-10-08 2008-10-09 2008-10-10 2008-10-11 2008-10-12 2008-10-15 
         2          3          1          3          1          1          1          1

情報としてはこれで十分。あとは

> barplot(tb)

でヒストグラムが出せるし、

for(i in seq(as.Date("2008-09-29"), as.Date("2008-10-10"), 1)){
    class(i) <- "Date"
    i <- as.character(i)
    s <- ifelse(is.na(tb[i]), 0, tb[i])
    writeLines(paste(i, paste(rep("*", s), collapse="")))
}

こんな感じの使い捨てスクリプトで

2008-09-29 
2008-09-30 
2008-10-01 ******
2008-10-02 *****
2008-10-03 *****
2008-10-04 *******
2008-10-05 ******
2008-10-06 **
2008-10-07 ***
2008-10-08 *
2008-10-09 ***
2008-10-10 *

こんな出力が得られます。こうしてみると、10/1〜10/5に集中してる感じかな。あとは実際に自分が見る作品と、見る時間帯を調整して計画を立てますか。

追記

自分の視聴範囲を調べてから再計算すると、こんな感じ。

2008-09-29 
2008-09-30 
2008-10-01 **
2008-10-02 ***
2008-10-03 ****
2008-10-04 *****
2008-10-05 **
2008-10-06 **
2008-10-07 ***
2008-10-08 *
2008-10-09 *
2008-10-10

１〜４辺りか。なるほど。