![[PukiWiki]](https://xn--p8ja5bwe1i.jp:443/wiki/image/banner.png "[PukiWiki]")
- バックアップ一覧
- 差分 を表示
- 現在との差分 を表示
- ソース を表示
- バイオ・データ・マイニング/Rで独立成分分析する へ行く。
この記事はまだ書きかけです．
はじめに †
ここでは，Rを使って独立成分分析を行います．
準備 †
Rのインストールについては，次のページを見てください．
独立成分分析は，fastICAパッケージに含まれています． そこで，fastICAパッケージをインストールします． install.packagesコマンドを実行すると，パッケージをダウンロードするサーバーを聞かれますので，リストからJapanのどれかを選択します．
> install.packages("fastICA") > library(fastICA)
独立成分分析 †
独立成分分析 (Independent Component Analysis, ICA) は，多変量解析の一手法で，高次元の数値データを（独立に生起している）複数の成分に分解する手法です．
複数の人が同時に話をすると，音声信号が混ざってしまいます． これを複数のマイクで拾って，元の成分に分解するのによく使われています．
[math]n[/math] 個の信号源から発生し（て混ざっ）た信号 [math]s_1, \dots, s_n[/math] を [math]m[/math] 個のセンサーで観測し [math]x_1, \dots, x_m[/math] を得たとします．
独立成分分析は，観測した [math]m[/math] 個のデータ [math] x_1, \dots, x_m[/math] から [math]n[/math] 個の（独立な）信号源から発生した [math]n[/math] 個のデータ [math]s_1, \dots, s_n[/math] に復元します．
例として，2つの信号源（音源）からそれぞれ次のような信号（音） [math]s_1[/math] と [math]s_2[/math] が発生しているとします．
二つの信号源からちょうど同じ距離の地点でこの信号を観測すると，この信号が1:1の割合で混ざって観測されます． \[ x_1 = 0.5 s_1 + 0.5 s_2 \]
また，二つの信号源からの距離が4:1の地点でこの信号を観測すると，次のようになります． \[ x_2 = 0.4 s_1 + 0.1 s_2 \]
このような条件を満たす地点は二つの信号源を結ぶ直線と直行する直線上に無数に存在し，係数の大きさは信号源からどのくらい離れているかに比例します． ここでは，上のような係数だとします．
上の図の信号をこのようにして観測すると，次のような信号が得られます．
独立成分分析を行う †
標準で使えるirisデータセットに対して，主成分分析を行ってみましょう．
> data(iris)
このデータセットは，アヤメの種類（Species）を花びらの長さ（Sepal.Length），幅（Lepal.Width），がくの長さ（Petal.Length），幅（Petal.Width）によって分類する問題です． 長さと幅は連続値，種類はsetosa, versicolor, virginicaのいずれかをとる離散値です．
このデータセットには，setosa, versicolor, virginicaという3種類のアヤメについて，それぞれ50個ずつ，合計150個のデータが含まれています． ランダムに10個のデータを選択して，見てみましょう．
> iris[sort(sample(1:150,10)),] Sepal.Length Sepal.Width Petal.Length Petal.Width Species 4 4.6 3.1 1.5 0.2 setosa 22 5.1 3.7 1.5 0.4 setosa 65 5.6 2.9 3.6 1.3 versicolor 97 5.7 2.9 4.2 1.3 versicolor 100 5.7 2.8 4.1 1.3 versicolor 108 7.3 2.9 6.3 1.8 virginica 116 6.4 3.2 5.3 2.3 virginica 122 5.6 2.8 4.9 2.0 virginica 136 7.7 3.0 6.1 2.3 virginica 146 6.7 3.0 5.2 2.3 virginica
> iris.ic <- fastICA(iris[1:4], 4)
> pairs(iris.ic$S, pch=as.integer(iris$Species))
