多次元データをクラスタリングする際に、それらのデータを2次元データに落とし込んで可視化させたいことがあります。そんな時に便利なのが「多次元尺度法」という手法です。 個々のデータ間の距離／類似度が分かっている場合に、それらのデータの座標を求めて、データ構造を復元するようなものです。 詳しい説明は割愛します。知りたい人はwikipediaと金先生の連載を読んで下さい。 体で覚えるタイプなので、とにかく何かデータを処理してみます。 「山手線」の地図を再現 さっそく試してみます。 山手線の各駅同士の直線距離を測っておいて、そのデータから実際の位置関係を復元できるか実験してみます。 山手線全駅の距離を測るのはめんどいので、適当に抜粋してしらべました。 以下のような表になりました。単位はメートルです。 さてさて、この距離表からどのようなデータ構造が再現されるでしょうか？ このデータを統計解析ソフトRで

モデル式 2.01×がく片長-12.57≧0のときバージニアアヤメと判別 2.01×がく片長-12.57＜0のときヘンショクアヤメと判別 （このモデル式では、バージニアアヤメは標本50個中37個、ヘンショクアヤメは50個中36個が正しく判別されている。） 最尤推定値は漸近的には正規分布することが知られている。今回の標本50個ずつのデータで出した最尤推定値（切片: −12.57、がく片長の係数: 2.01）が、どの程度正規分布に近いか、ブートストラップ法で以下のように調べることができる。 元データから n 個の標本を復元抽出する。このとき n は元データの標本数である。 最尤法でロジスティック回帰モデルに当てはめる。 このブートストラップ抽出を何度も（B 回）繰り返す。 こうして計算された「推定量の標本分布」は、本来の標本分布の近似になっている。 下図は10000回のブートストラップ抽出によ

Tsukuba.R#9での発表に向けてggplot2の勉強を始めた．これまで調べたことを軽くまとめてみる． ggplot2とは 高度な2Dグラフィックの描画について統一的な操作を提供してくれるライブラリ．「The Grammar of Graphics」という思想に基づいていて，ユーザがデータのどの部分に注目したいかを記述していくだけで簡単にきれいなグラフを得られるようになっている． http://had.co.nz/ggplot2/ ggplot2のインストール いつも通り． > install.packages("ggplot2") ggplot2で図を描画する ggplot2で図を描画する方法は主に以下の2通り． qplot関数を使う 自分が得たいグラフに合わせてレイヤーオブジェクトを組み合わせる ここでは以下のggplot2パッケージ付属のdiamondsサンプルデータを使い，それ

ツイッターでggplot2で日本語使えるの？的な話があったので、昔書きかけた記事を中途半端に描き上げてみました。 結論から言うと使えるんですが、OSとかデバイスが何か、とかでいろいろ違うので、完全に理解するのは結構大変です。 あとRのグラフィックスはbaseベースとgridベースがあるんですが、入門書では絶対にbaseベースなので、gridベースになじめない、このあたりも誤解を招く原因になってるような気がする。描画の自由度はgridベースのほうが大きいです。お手軽なプロットならbaseだったんだけど、gridベースのグラフィックスパッケージが出てきてgridでもお手軽にできるようになりつつある、という現状なのかな。 おまじない 面倒な解説は不要、という人は以下の例で。 library(ggplot2) d <- data.frame(x=1:2, y=1:2, t=c("あ","い"))

毎日ブログを書くようにしたら体調が良くなった ある日ふと、「最近の自分はアウトプットに臆病になっているなー」と感じ、こういうエントリを書いた。 daiksy.hatenablog.jp 内容は、雑なアウトプットで「だいくしーさん考えが浅いな」みたいな「化けの皮が剥がれる」ことを恐れているが、自分が過剰に意識するほど…

ここあママさんの「おかん秘伝の☆スペアリブ」の圧力鍋編です。 このレシピの生い立ち ここあママさんの「おかん秘伝の☆スペアリブ」がおいしすぎて、漬け込む時間を短縮できないかな～と思って、圧力鍋で煮込む方法を考えました。水以外は同じ分量です。 ここあママさんの「おかん秘伝の☆スペアリブ」の圧力鍋編です。 このレシピの生い立ち ここあママさんの「おかん秘伝の☆スペアリブ」がおいしすぎて、漬け込む時間を短縮できないかな～と思って、圧力鍋で煮込む方法を考えました。水以外は同じ分量です。

長いあいだ、プラモデルは表面だけを作るイメージがあったのです。中身は中空。 でも最近のプラモデルはそうじゃないからなあ。層を重ねて成形して、肌色の赤みを差したりするくらいだからなあ。

In 2010, each day for 30 consecutive days I attempted to make something creative. I hoped to force myself into making quick decisions, creating things instinctively. You can find below an index of all 30 Days in the Speed Creating project linked to all the day’s pages. For more information about the project visit this page. I exhibited some of this project in Milan, here are some pictures. You can