気象学研究室 - フォーラム

独立成分分析に使う基本的な数

投稿するにはまず登録を

投稿者スレッド

投稿日時: 2005-4-26 1:36

登録日: 2004-7-29

居住地: 地球

投稿: 303

独立成分分析に使う基本的な数

独立成分分析の結果の正当性を調べる際に必要な数値

非ガウス性(non-Gaussianity)は、これらの数値からのずれとして表現される。

1. 検証

> options(digits=22)
> x <- seq(-10,10,0.0001)
> sum( exp(-x^2/2) )*0.0001/sqrt(2*pi)
[1] 1.000000000000000222045

2. E{log(cosh(y))}

> x <- seq(-10,10,0.00001)
> sum( log(cosh(x)) * exp(-x^2/2) )*0.00001/sqrt(2*pi)
[1] 0.3745672074914379034993

3. E{-exp(y^2/2)}

> x <- seq(-10,10,0.00001)
> sum( -exp(-x^2/2) * exp(-x^2/2) )*0.00001/sqrt(2*pi)
[1] -0.7071067811865474617150

4. kurtosis

> sum( x^4 * exp(-x^2/2) )*0.0001/sqrt(2*pi)
[1] 3

題名	投稿者	日時
» 独立成分分析に使う基本的な数	森	2005-4-26 1:36

投稿するにはまず登録を