ryu9_2007_2

bloｇ開設した方は、ここに登録されているか確認してください。

学籍番号とリンク先が対応しているか確認してください。間違っていたら申し出てくさい。修正します。

～～～～～～～～～～～～～～～～～～～～～～
～～～～～～～～～～～～～～～～～～～～～～
Rを使ってみよう。

RWiki
・RWikiはRubyで開発されたWikiで、Rの日本語に関する情報のほとんどが網羅されている。
・Rの日本語文献は最近数多く出版されているが、
The　R　Book（WebCatPlus）
がおすすめ。琉大図書館にはないが、沖国、沖大、沖縄高専にはあるようだ。不思議。
他の参考文献に関してはRWikiを参考にしてください。

ところでSASにしろSPSSにしろ、Sも目の前のマシンにインストールしなければならない。Rもしかり。とことろがRには"R on Web"というサイトがあり、ある程度のことは体験することができる。
そこで、Rをインストールする前にブラウザさえあれば、どこででも実行できる"R on Web"を使ってみることにします。


R　on　Web

上記リンクを開いてください。こんなかんじになっていると思います。


この部分をマウスで範囲選択してコピーします。


今度はコードを記述するボックス部分に貼り付けて、ボックス下の「Submit」ボタンを押してください。


しばらくすると計算結果が出力されます。


コードの説明

重用なのは、言語仕様を押さえる、ということ。

x ＜- rnorm(100)
・ｒnormは正規乱数を発生させる関数で、ここでは100個の正規乱数を発生させ、ｘという変数に代入しています。
・代入は「＝」ではなく「＜－」で記述します。

y ＜- exp(x) + rnorm(100)
・ｘを指数変換し、それに100個の正規乱数を加える。こうしてできた系列をｙという変数に代入。

result ＜- lm(y ~ x)
・lmは線形回帰の関数で「ｙ＝a+ｂ・ｘ」という式を当てはめ、aとｂを推計する。推計結果をｒｅｓｕｌｔという変数の代入。
・このように推計式の場合、Rは「＝」でなく「～」とイコールでなく、チルダを用いる。

summary(result)
・summaryはその名の通り、結果を要約して出力する。回帰結果が保持されているresultの内容を表示する。

plot(x,y)
・ｘとｙをプロット（散布図）する。

abline(result)
・resultの中にある回帰係数aとｂを使った回帰線を描き、プロットグラフに追加する。

lines(lowess(x,y), col=2)
・Rには多くの平滑化関数があるが、その中の一つで、colは色を指定し、結果をlineで描画している。

hist(result$residuals)
・histはその名の通りヒストグラムを描画する関数で、回帰計算の結果が保持されているresultから残差系列を抜き出して、描画している。
・residualsはｌｍ関数の結果として出力される数値系列で、その名の通り残差系列である。
Rの書式では、ある処理を施したあと、その結果として出力される数値等に対して、「キャラクターダラー＄その変数名：ここではresult$residuals」というように指定することで結果の一部を取り出す。