2012-05-09から1日間の記事一覧

RでUNIXTIMEに変換

R

例１ そのままの場合 as.numeric(as.POSIXlt("2007-07-03 13:32:06 JST")) 例2 日付が2007.07.03だったら？ hidu<-paste(unlist(strsplit("2007.07.03","\\.")),collapse="-") as.numeric(as.POSIXlt(paste(sep="",hidu," 13:32:06 JST")) 例３ ベクトルを一…

Rでバーコード状のグラフを書く（ラグプロット）

R

データ， data 1の部分だけ黒い縦線を引く． バーコード状のプロット． 1がクラスタしているか確かめられる． rug関数を用いる． #data<-rbinom(size=1,prob=0.5,1000) plot(data,col=0) rug(which(data==1),ticksize=1)

Rで移動平均

R

Rで移動平均 moving_average<-function(ts123,lag){ tmp<-0 for(i in 1:lag){ tmp<-tmp+ts123[(i):(length(ts123)-(lag-i))]/lag } tmp<-c(rep(NA,lag-1),tmp) tmp } j<-moving_average(ll,32) plot(j)