2014-04-01から1ヶ月間の記事一覧

Rのlmとglmとglmnetの違い

SVMのカーネル比較

モデルによって最適なカーネルが異なる > kernel.vec <- c('linear','poly','radial','sigmoid') > scores <- data.frame() > for(word in kernel.vec){ + iris.svm <- svm(Species~.,iris,kernel=word,cross=5) + prediction <- predict(iris.svm,iris[,1:4…

2峰性のGMM(混合ガウス分布)をRのEMアルゴリズムで求める

http://stats.stackexchange.com/questions/59237/applying-bayes-estimating-a-bimodal-distributionmixtoolsを使う library(mixtools) data(faithful) plot(density(faithful$waiting), las=1, col="steelblue", lwd=2, main="") out <- normalmixEM(faith…

Rでbayesmixパッケージを使う

Rjagsのエラーが出る http://dr-k-lo.blogspot.jp/2013/03/problems-loading-jags-into-r.htmlパッケージインストーラでrjagsを入れてからbayesmixを入れる library(rjags) #success install.packages('bayesmix')

階層ベイズモデル

メモ http://www.slideshare.net/yokkuns/r-16189019 あと久保さんのPDF

Rで最大となる列名を取り出す

>head(df) setosa versicolor virginica 1 1 0 0 2 1 0 0 3 0 1 0 4 0 1 0 5 0 0 1 6 0 0 1 上の行列を下のベクトルに変形したい (最大となる列名を取り出す) 1 setosa 2 setosa 3 versicolor 4 versicolo 5 virginica 6 virginica max.colを使う colnames…

Rで使える機械学習の手法の比較

使うのはアヤメのデータ(iris) 縦軸: Sepal.Width, 横軸: Sepal.Lengthでプロットする ggplot(iris,aes(x=Sepal.Length,y=Sepal.Width,color=rev(Species)))+geom_point() k-meansクラスタリング R言語プログラミング: クラスター分析 - k-means http://d.h…

イェンスのブログの記事をMDSで可視化

library(tm) library(RMySQL) dbconnector <- dbConnect(md,dbname='jens_blog',user='root',password='******') jens.table <- dbGetQuery(dbconnector,'set names utf8') jens.table <- dbGetQuery(dbconnector,'select * from articles') all.jens <- as.…

ブログエントリの関連度をMDSで処理して図示

解析対象はshoheiaoki.comの記事 library(RMySQL) library(RMeCab) md <- dbDriver('MySQL') dbconnector <- dbConnect(md,dbname='ghana_blog',user='root',password='*****') shohei.table <- dbGetQuery(dbconnector,'set names utf8') #utf8に変換 shohe…

RMeCabで単語ベクトルからTDM(単語文書行列)を作成する

docMatrixDFを使う hoge.tdm <- docMatrixDF(hoge.vec)

RMySQLで日本語が文字化けするエラー

http://qiita.com/14021939/items/ccd5284070ce4f448a12 で、日本語が文字化けするSET NAMES utf8 を実行すればいいらしい http://stackoverflow.com/questions/12869778/fetching-utf-8-text-from-mysql-in-r-returns > library(RMySQL) > md <- dbDriver("…

rpartで決定木を綺麗に書く

オプションがわかりにくい uniform=T use.n=T pretty=1 あとprettyでラベルを表示する具合に応じてmarginを指定する (ラベルの文字列が長いとはみだすので) > z <- data.frame(Titanic) >Titanic1 <- data.frame(Class = rep(z[, 1], z[, 5]), Sex = rep(z…

ggplotで複数のグラフを重ねる

データは以下のもの >head(df) Ozone Solar.R Wind Temp Month Day Date 1 41 190 7.4 67 5 1 2014-05-01 2 36 118 8.0 72 5 2 2014-05-02 3 12 149 12.6 74 5 3 2014-05-03 4 18 313 11.5 62 5 4 2014-05-04 5 NA NA 14.3 56 5 5 2014-05-05 6 28 NA 14.9 6…

Rで列を削除する

colnamesと%in%を上手く使う もっと良いやり方がありそう(-演算子を使ったり?) (追記:この例だとdf[,1:4]のスライシングでできた)dfは以下の記事の続き http://d.hatena.ne.jp/seinzumtode/20140419/1397876421 >head(df) Ozone Solar.R Wind Temp Mon…

Rで文字列から日付をつくる

sprintfを使って文字列を作り、それをas.Dateで日付に変換するairqualityを使う >data(airquality) >df <- airquality >head(df) Ozone Solar.R Wind Temp Month Day 1 41 190 7.4 67 5 1 2 36 118 8.0 72 5 2 3 12 149 12.6 74 5 3 4 18 313 11.5 62 5 4 5 …

Shinyの基礎

ui.R library(shiny) shinyUI(fluidPage( titlePanel(h1("あいりすぷろったー ver.1.0",align="center")), br(), h4("@moyahima",align="center"), br(), h3(textOutput("hogeText"),align='center'), br(), fluidRow( column(9,offset=2,mainPanel(plotOutp…

shinyでggplot2を使う

ここを参考に http://www.stanford.edu/~cengel/cgi-bin/anthrospace/building-my-first-shiny-application-with-ggplotui.R library(shiny) shinyUI( pageWithSidebar( headerPanel("iris shiny"), sidebarPanel(), mainPanel(plotOutput("irisPlot")) ) ) …

Mac OSXにrsrubyをインストール

http://hrst.hatenablog.jp/entry/20100204/1265288817Rを普通にインストール環境変数を追加 export R_HOME=/Library/Frameworks/R.framework/Resourcesgemでrsrubyをインストール $ gem install rsruby -- --with-R-include=/Library/Frameworks/R.framewor…

Chrome extensionでクロスドメイン通信を行なう

あとで読む http://shokai.org/blog/archives/5387

jsshellが便利すぎる

jQueyrのターミナル jsshell https://chrome.google.com/webstore/detail/jsshell/kmgmkbicahmbceidoidjbkbpkfogaldh使い方 http://www.chrome-life.com/chrome/3792/ http://tera911.blog105.fc2.com/blog-entry-165.html カスタムスクリプトでSSL-VPNの自…