滑走窓方式データ整形 Rscript

滑走窓方でなくても使えるスクリプトです。

とりあえず部分時系列を得たい方推奨です。

ans<-spliting(data,1,3500,20)

kmeans(ans,3)

data:時系列データ。複数行不可 

start:与えた行列のスタートしたい地点。1が推奨

end:与えた行列の終わりにしたい地点。行列最終行推奨

point:分割地点、僕の場合は20とか15とか

分割で余った部分は全て切り捨てられる仕組みです。

最終的にデータフレームにしています。それとラベルが好みでない方はご自分で。

spliting <- function(data,start,end,point){

ncolnum<-end%/%point
print(ncolnum)
kazu<-point*ncolnum

#program return this array
#moto <- matrix(start:kazu, nrow=point, ncol=ncolnum)


moto<-c()
#print(moto)

kaisu <- 1
name = c()

ma<-c()

for(ii in 1:ncolnum){


newcol<-paste("R",ii)
name <- rbind (name, newcol)



mb<-c() #初期化

for(iii in 1:point){
#moto[ii,iii]<-data[kaisu]
kaisu<-kaisu+1
mb<-c(mb, data[kaisu])
}

moto<-rbind(moto,mb)

}

rownames(moto) <- c(name) # 列の名前
# 列の名前
moto <- as.data.frame(moto)
return(moto)

}