2018-09-01から1ヶ月間の記事一覧

多変量回帰と相関

多変量解析では、 (0)単独変数の調整 (0-1)カテゴリ化すると楽 cut関数 (0-2)少なすぎる変数の除去 (1)説明変数間の相関を調べる →とりあえず、主要っぽい変数を10個くらいとりだすとわかりやすい http://kusanagi.hatenablog.jp/entry/2014/07/23/2…

Rでgroup lasso の実験 一般化線形モデルの正則化

R

x1<-sample(1:100) x2<-sample(1:100) y<-2*x1+0.01*x2+rnorm(100) p<-data.frame(y,x1,x2) lm(y~x1+x2,data=p) f<-function(r){ sum((y-r[1]*x1-r[2]*x2)^2) } optim(c(1,3),f) ##Lasso f2<-function(r){ sum((y-r[1]*x1-r[2]*x2)^2)+20000*sum(abs(r)) } o…

EXELEファイルのRの読み込み。着工件数の年次データをえる

R

特定の都道府県の賃貸物件の着工数をえる。 国交省の【住宅】利用関係別 構造別 建て方別 都道府県別 戸数 (Excel ファイル 3,225KB) をダウンロード http://www.mlit.go.jp/sogoseisaku/jouhouka/sosei_jouhouka_tk4_000002.htmlreadxlパッケージの利用 htt…