fICA: Symmetric and Deflation-based FastICA Methods for Independent...

In fICA: Classical, Reloaded and Adaptive FastICA Algorithms

Description

The symmetric and deflation-based FastICA methods for the independent component problem. The function estimates the unmixing matrix using a single nonlinearity g.

Usage

fICA(X, g="tanh", dg=NULL, G=NULL, init=NULL, n.init = 1, method="sym2",
     inR=TRUE, eps=1e-06, maxiter=1000)

Arguments

X	a numeric data matrix. Missing values are not allowed.
g	the nonlinearity, tanh by default, see details.
dg	the first derivative of the nonlinearity, see details.
G	the integral function of the nonlinearity, see details.
init	a numeric matrix for the initial value of the algorithm
n.init	a positive integer for the number of initial values in symmetric algorithms, see details.
method	squared symmetric ("sym2"), symmetric ("sym") or deflation-based ("def")

.

inR	a logical which indicates whether R or C is used for computations, see details.
eps	convergence tolerance.
maxiter	maximum number of iterations.

Details

The deflation-based FastICA estimate depends on the initial value of the unmixing matrix. This is due to the fact that the algorithm does not always find the global maximum of the objective function, but it may stop to local maxima as well. Hence, the deflation-based FastICA estimate is affine equivariant only if the initial value is affine equivariant. Therefore, we recommend the use of reloaded_fICA or adapt_fICA.

The standard nonlinearities can be chosen by g="pow3", g="tanh" or g="gaus". These estimates can be computed either in R or in C (except the squared symmetric). In order to use some other nonlinearity, one has to give the nonlinearity and its derivative as functions (and the integral function when squared symmetric FastICA is computed), and the computations have to be performed in R.

If the symmetric or squared symmetric algorithm does not converge, one can choose n.init=k for some k>1. Then up to k random initial values are tried, and if none of them gives a convergent run of the algorithm, the function returns the matrix which gave the largest value of the objective function over all k\times maxiter steps.

Value

A list with class 'bss' containing the following components:

W	estimated unmixing matrix.
g	nonlinearity used.
method	symmetric or deflation-based.
S	estimated source components standardized to have mean 0 and unit variances.

Author(s)

Jari Miettinen

References

Hyvarinen, A. and Oja, E. (1997), A fast fixed-point algorithm for independent component analysis, Neural Computation, vol. 9, 1483–1492.

Miettinen, J., Nordhausen, K., Oja, H., Taskinen, S. and Virta, J. (2015), The squared symmetric FastICA estimator, Signal Processing, vol. 131, 402–411.

See Also

adapt_fICA, reloaded_fICA, nonlinearities

Examples

# creating some toy data
A<- matrix(rnorm(9),3,3)
s1 <- rt(1000,6)
s2 <- rexp(1000,1)
s3 <- runif(1000)

S <- cbind(s1,s2,s3)
X <- S %*% t(A)

# tanh is the default nonlinearity
res1<-fICA(X,method="sym")
coef(res1)
plot(res1) 
require(JADE)
MD(coef(res1),A)

# nonlinearity pow3 is chosen as follows
res2<-fICA(X,g="pow3",method="sym")
coef(res2)
require(JADE)
MD(coef(res2),A)

# nonlinearity from gf is chosen as follows
res3<-fICA(X,g=gf[[6]],dg=dgf[[6]],method="sym")
coef(res3)
require(JADE)
MD(coef(res3),A)

fICA documentation built on Dec. 11, 2021, 9:58 a.m.