« PNAS | 最新のページに戻る | Science 4/16 »
■ 来訪者プロット
来訪者数のプロットを作って、RのGAM(general additive model)でfittingしてみた。横軸が日を、縦軸が来訪者数を示す。クロスが各時点での値。Fittingした値とプラスマイナスseを曲線で表示。GAMでの計算は以下のとおり。
# diary.txtとして一列目に日記開始からの日数を、二列目に各日の来訪者数のファイルを作成しとく。
library(mgcv)
DIARY<-read.table("diary.txt",header=FALSE)
x<-DIARY[,1];
y<-DIARY[,2];
b3=gam(y~s(x),family=poisson)
pred3<-predict.gam(b3,se.fit=TRUE)
se3a<-exp(pred3$fit+pred3$se.fit)
se3b<-exp(pred3$fit-pred3$se.fit)
write.table(cbind(x,se3a,b3$fitted.value,se3b),file="diary3.txt", quote = FALSE, sep = " ")
めちゃシンプル。ほとんどデフォルトの設定のままだが、縦軸は0以上しかないので、Poissonモデルを使用した。GAMはあまり細かい変動には影響されずにfittingしてくれるので都合がよい。データからは、いくつか山があって、いまは多少停滞しているということがわかる。