Rことはじめ（自分用メモ）

ツマグロヒョウモン交尾

proxy

proxyが必要な場合は，起動時に --internet2 という引数をつける
- http://cse.naro.affrc.go.jp/takezawa/r-tips/r/08.htmlより
- Sys.putenv("http_proxy"="http://プロキシーの IP アドレス:ポート番号") を実行と書いてあるページもあったが*1，うまくいかなかった．

エクセルのファイルを直接読み込むには

http://cse.naro.affrc.go.jp/takezawa/r-tips/r/41.html

chooseCRANmirror()
install.packages("gregmisc")
install.packages("gmodels")

メニューからディレクトリ変更してから，

library(gregmisc)
x <- read.xls("Eh_climate.xls", sheet=1)

ここで，エクセルファイルの1行目に日本語が入っていたら

5 は不正なマルチバイト文字です

というエラーになったので，1 byte文字に書き直して再度実行．

重回帰

result1 <- lm(LogSW ~ Alt+Lat+Lon+Temp+Rain+haplotype, data=x)
summary(result1)
result2 <- step(result1)
summary(result2)

とかそんな感じ．カテゴリー変数もそのまま扱えるっぽい．ただ，JMPでいうところの「by」による分割はどうするんだろう？

HN <- split (x, x$Geni.type)
resultH <- lm(LogSW ~ Alt+Lat+Lon+Temp+Rain+haplotype, data=HN$H)
resultN <- lm(LogSW ~ Alt+Lat+Lon+Temp+Rain+haplotype, data=HN$N)

かな？
RでもJMPと同じで標準偏回帰係数は出してくれないのか．標準偏回帰係数 $\hat{\beta}_j^\ast = \hat{\beta}_j \sqrt{s_{jj}/s_{yy}}$ （ $s_{jj}$ はxの分散， $s_{yy}$ はyの分散）で求められるらしい*2．