R

グラムシュミッドの正規直交化を書いてみてRcppを比較する

Eigenあっ, どーも僕です. そういえば, 理工系学者, 学生にとっては垂涎ものであるFreeFEMの日本語参考書がでましたね. FreeFEMの参考書は言語を問わず少ないので, 工学部系のほとんどの研究室で購入しているんじゃないでしょうか. Rcpp sourceCpp 訂正:cpp…

アルマジロを可愛がることにした

Armadillo: C++ linear algebra library あっ, どーも僕です. 数あるRcpp系パッケージのなかでRcppArmadilloを少し勉強したのでまとめます. RcppArmadilloによりRcpp経由で, C++ライブラリのArmadilloを使用できます. 以下にRcppArmadilloの機能をまとめます…

Rで実装する多層パーセプトロン

R

あっ, どーも僕です. 修士論文書きました. 公聴会まで少し時間があるのでブログ書きます. 本題に入る前に Rのdplyrパッケージが爆速なのですね. デフォルトパッケージ推しのわたしもさすがに乗り換えました. ところで, dplyrで外生的に分解軸を設定するには…

Rでk-means法とその拡張3 改良k-means編

R

あっ, どーも僕です。 Rでk-means法とその拡張の最後は, 改良k-meansです。概要は前々回の記事をみてください。また, アルゴリズムは論文をみてみてください。フリーで入手できます。さっそく, 結果ですが, クラスター数を3にしたときがこちら。相変わらずir…

Rでk-means法とその拡張2 x-means編

R

あっ, どーも僕です。注意!!本エントリーの記述はRソースを含めて全面的に書き直しています!!(なぜかこのエントリーはアクセス数が多いので(汗))@2014/9/11Python版はこちらから.github.com クラスター数が自動で決まるkmeans法 前回はk-means法の実…

Rでk-means法とその拡張1 k-means編

R

あっ, どーも僕です。 先週, 「明日から頑張る」と宣言したはずの後輩が, その日から一度も研究室に来てくれません。 まだ, 卒論に書くことないがほとんどはずなのに.....先輩としてどうしたらいいものか。 k-means法によるクラスリングの拡張 k-means法とそ…

Rで普通の数値解析その4

R

加藤あいの結婚に衝撃を受けたのですが, みんな関心なくてひとり浮いてました。あっ, どーも僕です。会社員と会社役員の差って思ったよりもデカそうです。 追記 スクリプトの最新版はこちらから。 Rで普通の数値解析 Rで普通の数値計算その4 やっぱりeigen…

Rで普通の数値解析その3

R

あっ, どーも僕です. そういえば, 2019年にハンドボール女子世界選手権が日本開催に決定したのでここで宣伝しときますね!! 日本代表女子の試合はyoutubeに動画が上がっています!! ハンドボールをもっと楽しみたい人は本場のハンドボールを見てください!…

Rで普通の数値解析その2

R

あっ, どーも僕です。計算による線形代数 (工系数学講座 2)に次のように書かれていました。 初心者が作成したプログラムは, 計算精度や計算コストの点で十分な精度を得ることができないことが多い. したがって, 実際に利用する場合には, 自信で作成するより…

ガッキーより竹内結子

R

なのでダンダリン見てる。 あっ, どーも僕です。リアップは劇薬じゃなくてよかったです。 rstanでちょこちょこ こちらのブログにrstanの導入&入門編が紹介されていました。MCMCの計算にStanを使ってみた(超基礎・導入編) - 銀座で働くデータサイエンティ…

Windows7にRStanをインストールするためにしたよくわからんこと

R

あっ, どーも僕です。 Stan2を使って自分もベイジアンの端くれになりたいのですが、RStanがなかなかインストールできずに苦戦していました。JAGSでもいいのですが、いまいち自分がやりたいことが実装できなかったのでなんとかStanに移動したかったです。イン…

Rで普通の数値解析その1

R

塩には必要だが。あっ, どーも僕です。 来年から おかげさまで、来年からの食い扶持が決まりました。 特に@beroberoさんは見ず知らずの私に声をかけてくださりありがとうございました。励みになりました。内定先ではデータ解析とかとはあまり関係ない仕事か…

最近、日本の出版社からR本出過ぎ

R

海外の出版社もですが。 あ、どーも僕です。 上級ハンドブックの目次みれるようになりましたね。 角度データなる初耳の単語があったり、苦手なネットワーク関連があったりしてとても楽しみです。 これ出る前にワンちゃん本読み終わりたいです。R言語上級ハン…

Rで確率的ニューラルネットワーク

R

滑舌悪いので集団討論でアウト. あっ, どーも僕です. 昨日Tokyo.Rがあったらしいですね. 一度は参加したり発表してみたいと思いつつも, 片道10時間の7000円なのでちょっと厳しいですね. 発表のネタとしてはこないだ自分の研究のために適合的棄却サンプリング…

満島さん薄幸似合いすぎ. フォルダ5聴き直すわ.

R

あっ, どーも僕です. 下痢と肩こりが最近の悩みです. flexmixとかいうパッケージすげぇ 昨日モデルベースのEMアルゴリズムを使った論文を読みました. EMアルゴリズムは一年くらい前に勉強して, 結局実装できなくて苦手意識バリバリだったのでうげええええっ…

RでタグチメソッドのMTA法を組んでirisで試してみる

R

水着よりもいまは浴衣姿を拝みたい. あっ, どーも僕です. みなさんこちらの記事見ましたか? Webデータ分析&データサイエンスで役立つ統計学・機械学習系の分析手法10選 - 道玄坂で働くデータサイエンティストのブログRをはじめて一年あまりだそうです. すご…

Rで95%信頼区間の95%具合を感じてみる

R

大学院までいかせてもらってるのに就職できないオレwwwwwwwwwwwwwwwwあっ, どーも僕です. また, 落ちましたとさ. 誰か助けて. 95%信頼区間のシュミレーション t検定における信頼係数αとした信頼区間のお話. わたしが通っている学科は土木系…

平均値0, 分散1にした対称な連続分布いろいろ

R

ハゲてきたので, 禿げ上がる前にもう一度告白したい. あっ, どーも僕です. 異常値検知は正規分布でイカンのか? 対照な分布でも正規分布とは限らないって話. こないだ, はてなのホットエントリーに2σ管理による異常値検知の話が上がって来ました. 「n日間移…

Rで指数乱数を使って中心極限定理を感じてみる

R

エレベーターに酸っぱい匂いが充満してました. あっ, どーも僕です. 中心極限定理 中心極限定理とは, 同一の確率分布から得られた独立な確率変数の期待値の分布は, 漸近的に標準正規分布に従うように変換できるというあれです. で, 統計学の本を学びなおして…

3つの標準正規乱数の発生方法

R

嫌いな女性のタイプは靴の踵を踏んじゃってる人. あっ, どーも僕です. 3つ手法による標準正規分布の発生 こないだの本に従って, 一様乱数が得られるという仮定のもと, 3つ方法で標準正規乱数を発生させてみました. めんどくさいので, 理論も数式も書きません…

相関がある二変量正規乱数を発生させてみる

R

あっ, どーも僕です. やたら同級生(女)に相談されるのでチャンスかと勘違いしちゃいました(・∀・)友達一人減ったかしら(・∀・) 相関がある二変量正規乱数 今日家に帰ってきたら中古で2500円だったのでついポチってしまいまった本が届いていたので, さっそく…

Rを勉強するときに買ってよかった成書5冊

R

好きなタイプは日南響子と答えてます. あっ, どーも僕です. 春だ!Rを勉強しよう! Rを真面目に使いはじめて一年半くらいです. 4月になり研究室の後輩にRを教えるにあたりオススメの本を何かと考えたので載せておきます. これらを読むことでR入門からR中級く…

Rのplyrパッケージをちょっとだけ復習

R

第一志望落ちて正直死にたくなりましたがあまちゃんが可愛くて面白くてもっと見たくて生きています. あまちゃんがこの時期に始まらなければ学部2年の時のように黒歴史に突入でした. あっ, どーも僕です. plyrパッケージ 復習 Tokyo.Rの演題に「@teramonagi p…

R 3.0.0のインストールとライブラリーのアップデート(泥臭く。。。。)

R

あっ, どーも僕です. R 3.0.0出ましたね!!!アップデートの方法がいままでとちょっと違ったので誰かの参考になればと思い載せておきますね. といっても, すでにRPubsに紹介されていたので, 私がやった泥臭い版を. Rのアップデート installrというパッケージで…

そろそろ同級生の結婚が増えてきたので母親からの圧力やばい

R

オレが高校生の恋心をいつまでも引きずっている ストーカー野郎だということを 一生独身が3割を超える現代社会の流れに乗っていることを 母はまだ知らない....... あっ, どーも僕です. カレーづくり 実は結城浩先生のカレーの問題に挑戦していましたのです…

恋愛禁止されていないけど

R

23年間自粛中です. あっ, どーも僕です. 就活諦めモードです. お腹いたいです. R RbloggersでtwitteRの使い方が載っていので なぞったのを載せたいと思います. twitteRの存在は前から知っていたのですが, cygwinというものを知らず windowsのRstudio上ででき…

ハサミを使うタイミングをやってみた

R

あっ, どーも僕です. こちらの問題がおもしろそうなのでやってみました. 第6回 ハサミを使うタイミング~どうやって端を見つけるの?─結城浩からの問題|gihyo.jp … 技術評論社 ただ, 解答期間が過ぎていたので提出できないです. 残念. まぁ方針はゴリ押しな…

Rで「フィールドデータによる統計モデリングとAIC」のP.95のシュミレーションをしてみた

R

鳴らない携帯に月7000円を払うなど. あっ, どーも僕です. 現実が地獄でとにかく逃避したくてタイトルにあるシュミレーションをRでやってみました. フィールドデータによる統計モデリングとAIC (ISMシリーズ:進化する統計数理)作者: 島谷健一郎出版社/メーカ…

RでLagrange polynomial interpolation

R

私はまだ知らない. 尊敬している人の鼻毛がでているときどう伝えればいいのか. あっ, どーも僕です. Rでラグランジュ補間多項式 研究の息抜きにラグランジュ補間多項式のRでの仕方を調べてみました. まぁ息抜きというか現実逃避ですかね. ・・・ラグランジュ…

RでLepage Test

R

血尿でるまで研究頑張ろうとおもったけど先にア●ルから血がでてきました. あっ, どーも僕です. みなさんは元気ですか. Lepage Test レポートには課されていないのですが, マンケンドール検定とともにせっかく教わったので(・・・というか研究がうまくいかな…