Rの関数いくつか - 日常メモ

今日は規則性の無いデータに無闇に関数を作用させるだけの話。
なので、内容が無いよーうです。

> setwd("D:\\100_work\\data")
> data<-read.csv("data.csv")
> head(data)
  No       Date One Two Three Four Five Six Option
1  1  2000/10/5   2   8    10   13   27  30     39
2  2 2000/10/12   1   9    16   20   21  43      5
3  3 2000/10/19   1   5    15   31   36  38     13
4  4 2000/10/26  16  18    26   27   34  40     13
5  5  2000/11/2   9  15    21   23   27  28     43
6  6  2000/11/9   6  12    23   25   28  38     22
・・・・・・・・・・・・・・・・・・・・・・・・・

それぞれを別々の変数に格納。

> data1<-data[3]
> head(data1)
  One
1   2
2   1
3   1
4  16
5   9
6   6
> data2<-data[4]
> data3<-data[5]
> data4<-data[6]
> data5<-data[7]
> data6<-data[8]
> dataOption<-data[9]

時系列データに変換。

> data1.ts<-ts(data1)
> data1.ts<-ts(as.numeric(data1.ts))
> data2.ts<-ts(data2)
> data2.ts<-ts(as.numeric(data2.ts))
> data3.ts<-ts(data3)
> data3.ts<-ts(as.numeric(data3.ts))
> data4.ts<-ts(data4)
> data4.ts<-ts(as.numeric(data4.ts))
> data5.ts<-ts(data5)
> data5.ts<-ts(as.numeric(data5.ts))
> data6.ts<-ts(data6)
> data6.ts<-ts(as.numeric(data6.ts))
> dataOption.ts<-ts(dataOption)
> dataOption.ts<-ts(as.numeric(dataOption.ts))

描画領域を7行1列に分割して、行の左から順番にグラフを描く。
mfrowパラメータで分割する行数・列数を指定し、pltパラメータで余白部分を調整。

> par(mfrow=c(7,1),plt = c(0.1, 0.98, 0, 0.95)) 
> plot(data1.ts, ylim=c(0,50))
> plot(data2.ts, ylim=c(0,50))
> plot(data3.ts, ylim=c(0,50))
> plot(data4.ts, ylim=c(0,50))
> plot(data5.ts, ylim=c(0,50))
> plot(data6.ts, ylim=c(0,50))

ヒストグラムを描く。

> par(mfrow=c(7,1),plt = c(0.1, 0.98, 0, 0.95)) 
> hist(data1.ts, xlim=c(0,45),ylim=c(0,600),xaxp=c(0,40,40),breaks=45,col="orange")
> hist(data2.ts, xlim=c(0,45),ylim=c(0,600),xaxp=c(0,40,40),breaks=45,col="orange")
> hist(data3.ts, xlim=c(0,45),ylim=c(0,600),xaxp=c(0,40,40),breaks=45,col="orange")
> hist(data4.ts, xlim=c(0,45),ylim=c(0,600),xaxp=c(0,40,40),breaks=45,col="orange")
> hist(data5.ts, xlim=c(0,45),ylim=c(0,600),xaxp=c(0,40,40),breaks=45,col="orange")
> hist(data6.ts, xlim=c(0,45),ylim=c(0,600),xaxp=c(0,40,40),breaks=45,col="orange")