傀儡師の館.Python

Ruby

(0)

甘い物

Profile

kugutsushi

フォローする

Free Space

設定されていません。

< 新しい記事

新着記事一覧(全1864件)

過去の記事 >

2007.10.02

R を使い始める (Swivel への言及もあり)

カテゴリ：カブロボ

しばらくオープンソースの統計解析システム R を地道に使ってみることにした。

R とは (RjpWiki)

TPCI - TIOBE Programming Community Index を見ると 50位で 0.072% で、プログラミング言語としてはメジャーとはとらえられないものの、統計処理とかの世界では知らない人はいないだろう。ニッチのメジャーシステムというところか。実際にどんなグラフが描けるかについては R graphical manuals -- Image Browser のあたりを見ると分かる。

Swivel は期待が持てるサービスや統計情報のグラフ化とデータコミュニティなどで書いた Swivel では Inside Swivel: Graphing を見ると、R もグラフ作成のパッケージの一つとして当初検討対象にされたようだ。もっとも、プレゼンテーション面で融通が利かないところやたくさんの人がアクセスするようなウェブサイトで使うのには向いていないということで使われなかったようだが。他には、Ruby のグラフライブラリ Gruff も検討されたが大規模サイトには向かないということで選択肢からはずれている。 GNUplot も速度的にはよいがプレゼンテーション面で足りず、 Ploticus rploticus を使うようになったようだ。Swivel って Ruby 使っているのか。

その他、視覚化の観点からすると IBM's Many Eyes とか prefuse とかもいいらしい。確かにぱっと見た目良い感じ。

すでに値がある株価の動き等を時間軸に沿ってウェブサイト上で表示させるなら Timeplot もよさそうな感じ。 Energy Pricing みたいなものが簡単にできるようだから株価の推移などにはもってこいか。DHTMLベースで AJAX を使っていて、マウスカーソルでポイントしたところの付加情報を表示させるなどといったことができる。圧巻！AJAX+DHTMLのグラフライブラリ「Timeplot」に日本語の紹介がある。

グラフを描くということだけなら GnuPlot の方が使うのは楽だろうが、やりたいのは単にプロットするだけでなく計算もさせたいので R に再挑戦。前にもいっったん手をつけたことがあるのだが、途中でめげた。今度こそある程度のレベルまで使いこなす。日本語資料も豊富にあるので時間さえかければなんとかなるかなと。でも、最初は R Commander のような GUI ベースで学んでから ( Rコマンダー(1) 、 Rコマンダー(2) 等参照) を使ってから、生の R を使うとかいう方向の方が楽かもしれないとかも思う。とりあえずマニュアルの類をざっと見ながら動かしてみる。最近、頭が一段と固くなった。。。。。とりあえず、先はまだ長い。


MA <- function (x, day) {
results = c()
# 移動平均が出せないところは NULL を入れておく
for(i in 1:day-1) {
results[i] = NULL
}
# 移動平均を求めて結果に入れていく
for(i in day:length(x)){
start = i - day + 1
results[i] = mean(x[start:i])
}
results
}

試してみる。


> d = c(1:50)　# 1,2,3,4,5....50 のデータを代入する
> d
> MA(d, 5)
[1] NA NA NA NA  3  4  5  6  7  8  9 10 11 12 13 14 15 16 17 18 19 20 21 22 23
[26] 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48

> MA(d, 10)
[1]   NA   NA   NA   NA   NA   NA   NA   NA   NA  5.5  6.5  7.5  8.5  9.5 10.5
[16] 11.5 12.5 13.5 14.5 15.5 16.5 17.5 18.5 19.5 20.5 21.5 22.5 23.5 24.5 25.5
[31] 26.5 27.5 28.5 29.5 30.5 31.5 32.5 33.5 34.5 35.5 36.5 37.5 38.5 39.5 40.5
[46] 41.5 42.5 43.5 44.5 45.5

> summary(MA(d, 10))
Min. 1st Qu.  Median    Mean 3rd Qu.    Max.    NA's
5.5    15.5    25.5    25.5    35.5    45.5     9.0

# 日経平均の終値を使って 5日移動平均値をプロットしてみる
# nikkei225.csv の中に株価情報が入っている
> nikkei = read.table("nikkei225.csv", header=T) # データを読み込む
> summary(nikkei)
日付         始値            高値            安値
2007年10月1日: 1   Min.   :15477   Min.   :15907   Min.   :15262
2007年7月20日: 1   1st Qu.:16075   1st Qu.:16260   1st Qu.:15915
2007年7月23日: 1   Median :16460   Median :16532   Median :16366
2007年7月24日: 1   Mean   :16595   Mean   :16696   Mean   :16460
2007年7月25日: 1   3rd Qu.:16950   3rd Qu.:17037   3rd Qu.:16823
2007年7月26日: 1   Max.   :18149   Max.   :18223   Max.   :18125
(Other)      :44
終値
Min.   :15274
1st Qu.:16150
Median :16428
Mean   :16584
3rd Qu.:16920
Max.   :18158

> plot(MA(nikkei$終値, 5), ylab="5日移動平均", xlab="最近の50日",
+ main="日経平均の5日移動")

こうやってみると、16500円近辺をうろうろしていて、とりあえず 16920円超えができるかどうかで次のステージに入れるかどうかの展開が変わってきそうということで、今日は 17000円突破なのでとりあえずがんばれの展開か。とりあえず為替相場も今日は良い感じ。今週末まである程度のところで維持されればよいのだが。。。。なんとなく、明日以降は一服感でうろうろ相場になりそうな気もする。

いちおう、ダブルボトムは形成ってことで下値不安は低くなったのかな。再び異変がない限りは。しかし、こうやってみるとよく下げたもんだ。