読者です 読者をやめる 読者になる 読者になる

自分用メモ

問題発見と解決 R

自己相関係数とか偏自己相関係数とかを曜日ごとにplotしようかと思ってて、

par(mfrow=c(4,2))

とかやって、一画面に表示していたんけど、画像として吐き出すときにそれだと小さくなりすぎる。なので、一枚一枚プロットしようかと思ったけど、めんどいので、コード書いた。

week <- c("sunday","monday","tuesday","wednesday","thursday","friday","saturday")

for(i in 1:7){
  png(paste("acf_",week[i],".png",sep=""))
  acf(data[data$week==(i-1),]$count,main=paste(week[i],"の自己相関係数",sep=""))
  dev.off()
  png(paste("pacf_",week[i],".png",sep=""))
  pacf(data[data$week==(i-1),]$count,main=paste(week[i],"の偏自己相関係数",sep=""))
  dev.off()
}

後はflickrをコマンドラインで処理できるやつを使って

flickr_upload *.png;rm *.png

とかやればよい。Rは簡単なお仕事を簡単にできるので、いいですね。