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

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

投稿するにはまず登録を

投稿者	スレッド
森	投稿日時: 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(2pi) [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(2pi) [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(2pi) [1] -0.7071067811865474617150` 4. kurtosis `> sum( x^4 * exp(-x^2/2) )0.0001/sqrt(2pi) [1] 3`

投稿するにはまず登録を