キーワード調査の前に（２）

前のエントリーで、グラフを重ねたほうがいいと思ったので少し改造。こういうことができるからRは楽しい。

plot.datestep <- function(d, v){ #d:date(x) v:value(y)
    mar.old <- mar.new <- par("mar")
    mar.new[1] = mar.new[1] + 1
    mar.new[4] = mar.new[4] + 1
    color1 <- "#1e5692"
    color2 <- "#da2025"
    par(mar=mar.new)
    plot(d, v, type="s", col=color1, ylab="", lwd=2, axes=F)
    axis(2, col=color1)
    par(new=T)
    plot(d, cumsum(v)/sum(v), type="l", col=color2, ylab="", lwd=2, axes=F)
    axis(4, col=color2)
    axis.Date(1, at=seq(min(d), max(d), "2 days"), format="%m/%d", las=2)
    par(mar=mar.old)
}

l <- readLines("clipboard")
date <- strptime(l, "%m月%d日")
date <- as.Date(date[!is.na(date)])
plot.datestep(unique(date), table(as.Date(date)))

こんな感じ。左がカウント数、右が累積分布。分かりやすくなったような、分かりにくくなったような・・・。