統計学への誘い - かずのたしなみ（統計環境Rによる） @Wiki

################################################
##國友直人『現代統計学』（1992）第1章
################################################
##Step1　データの入力
################################################

##(a)作業フォルダ（ディレクトリ）を指定

getwd("/Users/osamu/Documents/R/kunitomo/chap1")

##(b)データ「表１．１　成績データ」を直接入力し、seisekiという名前を付ける。

seiseki<-c(50,82,70,54,75,72,70,33,80,52,98,73,62,73,43,86,93,85,61,89,98,28,87,100,50,87,66,82,85,70,75)

##(c)データ「表１．２　収益率データ」をcsvから読み込み，shuueki_rituという名前をつける

shuueki_ritu<-read.csv("data1.2.csv")

##(c)読み込んだデータを確認

seiseki
[1]  50  82  70  54  75  72  70  33  80  52  98  73  62  73   43  86  93  85  61  89  98  28  87 100  50  87
[27]  66  82  85  70  75

shuueki_ritu
   X X1980 X1981 X1982 X1983 X1984
1   1   9.2   2.8  -0.6  -2.8   0.0
2   2   2.3  -1.4 -11.8   9.3  -5.7
3   3  -6.5  17.6   3.5  11.4  10.6
4   4   9.0  17.8   1.9   3.0  -0.6
5   5   5.3   5.5  -5.5  -7.5 -11.2
6   6  -4.3  -1.9  -9.1   2.5  -3.8
7   7  -3.7   1.9  -5.7  -0.6  -5.2
8   8   7.0   9.0   2.3   1.8   6.2
9   9   7.6 -10.3  -4.9   5.1  -4.2
10 10   1.4 -10.3  -0.8  -2.3   2.1
11 11  -3.4  -7.7   8.0  -6.0   0.6
12 12   0.7   6.5   6.7  10.6   4.7

##(d)収益率データは時系列データなので、時系列データに変換する。
##まず各年の月毎のデータ（各列に入っている）をつなぎ合わせる

temp<-c(shuueki_ritu[,2],shuueki_ritu[,3],shuueki_ritu[,4],shuueki_ritu[,5],shuueki_ritu[,6])
temp
[1]   9.2   2.3  -6.5   9.0   5.3  -4.3  -3.7   7.0   7.6   1.4  -3.4   0.7   2.8  -1.4  17.6  17.8   5.5
[18]  -1.9   1.9   9.0 -10.3 -10.3  -7.7   6.5  -0.6 -11.8   3.5   1.9  -5.5  -9.1  -5.7   2.3  -4.9  -0.8
[35]   8.0   6.7  -2.8   9.3  11.4   3.0  -7.5   2.5  -0.6   1.8   5.1  -2.3  -6.0  10.6   0.0  -5.7  10.6
[52]  -0.6 -11.2  -3.8  -5.2   6.2  -4.2   2.1   0.6   4.7

##各年の月毎のデータをつなぎ合わせたものを、ts()を使って時系列データに変換する
##１年に12回（毎月）という頻度（freq=12）、最初は1980年（start=1980）という情報を与えてやる。

t_shuueki_ritu<-ts(temp,freq=12,start=1980)
t_shuueki_ritu
       Jan   Feb   Mar   Apr   May   Jun   Jul   Aug   Sep   Oct   Nov   Dec
1980   9.2   2.3  -6.5   9.0   5.3  -4.3  -3.7   7.0   7.6   1.4  -3.4   0.7
1981   2.8  -1.4  17.6  17.8   5.5  -1.9   1.9   9.0 -10.3 -10.3  -7.7   6.5
1982  -0.6 -11.8   3.5   1.9  -5.5  -9.1  -5.7   2.3  -4.9  -0.8   8.0   6.7
1983  -2.8   9.3  11.4   3.0  -7.5   2.5  -0.6   1.8   5.1  -2.3  -6.0  10.6
1984   0.0  -5.7  10.6  -0.6 -11.2  -3.8  -5.2   6.2  -4.2   2.1   0.6   4.7

→國友直人（1992）『経済学入門シリーズ　現代統計学（上・下）』（日経文庫）へもどる