![[PukiWiki]](image/pukiwiki.png "[PukiWiki]")
フリーの統計解析ソフト。普通にこの手のソフトを買うと10万円とかとられるので懐にやさしい。かなり高機能。ソフト名で検索しにくいのが最大の欠点。
以下のCRAN(The Comprehensive R Archive Network)ミラーのうち、お好きなところで。
以下を参照。「R-x.x.x-win32.exe」のような名前のファイル(32bit版Windowsの場合)をダウンロードして実行するだけ。
インストールしたら、必要ならR User Configurationで設定を変更しておこう。また、RコマンダーというGUI環境が用意されているので、使いたい場合は
install.packages("Rcmdr")
library(Rcmdr)
のように入力（インストールは初めてのときのみ）するとRコマンダーが起動する。最初の起動には分析に必要なパッケージを沢山ダウンロード＆インストールするので、多少時間がかかる。Rコマンダーについては
を参照。
他にもぐぐるといっぱい出てくる。CRANには他にもマニュアル等が
あたりに置いてある。英語が多いが、日本語・中国語などのもあり。本については
を参照。
「data.csv」というcsvファイルを「dataframe」という名前で読み込む。
dataframe <- read.csv("data.csv", stringsAsFactors = FALSE)
data.csvという名前で保存する場合、以下のようにする。
write.table(dataframe, "data.csv", sep=",", na="", quote = T, col.names = T, row.names = F, append = F)
column列が「hoge」となっている行を抽出。
dataframe2 <- subset(dataframe, column == "hoge")
with(dataframe, ・・・)
kaiki.kekka <- lm(y ~ x, data = dataframe)
「expand.grid」や「outer」を使うといいらしい。
標準のwilcox.testは近似計算になるらしいので、
install.packages("exactRankTests")
library(exactRankTests)
wilcox.exact(A群,B群)
のようにすると正確な値が出る。
これも古くなったらしく、代わりに
install.packages("coin")
library(coin)
wilcox_test(c(A群,y) ~ factor(c(rep("A群",length(A群)),rep("B群",length(B群)))), distribution="exact")
とするらしい。
必要なパッケージ類は
からまとめてインストールすれば楽。計量経済学のためには
を入れる。入れ方は
install.packages("ctv")
library(ctv)
install.views("Econometrics")
でおけ。
以下のように書く。
kaiki <- lm(y ~ x1 + x2 + x3, data = hogedata) coeftest(kaiki, vcov=vcovHC(kaiki, type = "HC1"))
オプションのデフォルトはHC3。HC1にするとStataと同じ結果を得る(Stata13と比較)。
plmを使う。モデルはwithin, random, ht, between, pooling, fdを指定できる。ランダム効果モデルではswar(標準), amemiya, walhus, nerlove, kinlaの5つを指定できる。
outcome <- plm(y ~ x1 + x2, data = hoge, model = "random", random.method = "amemiya")
のように書く。
ファイル - ファイルの保存 で保存できる。
適当なリンク。
