Learn R Programming

BSSasymp (version 1.2-3)

ASCOV_JADE: Asymptotic covariance matrix of JADE and FOBI estimates

Description

JADE solves the blind source separation problem in the case of independent components with at most one component having kurtosis values zero, while FOBI requires distinct kurtosis values. The functions compute the asymptotic covariance matrices of JADE and FOBI estimates for the mixing or the unmixing matrices.

Usage

ASCOV_JADE(sdf, supp=NULL, A=NULL, ...)

ASCOV_FOBI(sdf, supp=NULL, A=NULL, ...)

Arguments

sdf

a vector of density functions of the sources scaled so that the mean is 0 and variance is 1.

supp

a two column matrix, where each row gives the lower and the upper limit used in numerical integration for the corresponding source component which is done using integrate.

A

the mixing matrix, identity matrix as default.

arguments to be passed to integrate.

Value

A list with the following components:

W

mean of the unmixing matrix estimate.

COV_W

asymptotic covariance matrix of the unmixing matrix estimate.

A

mean of the mixing matrix estimate.

COV_A

asymptotic covariance matrix of the mixing matrix estimate.

EMD

The limiting expected value of \(n(p-1)MD^2\), see details.

Details

The order of the estimated components is fixed so that their kurtosis values are in a decreasing order. The signs of the components is fixed so that the sum of the elements of each row of the unmixing matrix is positive.

Since the unmixing matrix has asymptotic normal distribution, we have a connection between the asymptotic variances and the minimum distance index, which is defined as $$MD(\hat{W},A)=\frac{1}{\sqrt{p-1}} \inf_{P D}{||PD \hat{W} A-I||,}$$ where \(\hat{W}\) is the unmixing matrix estimate, \(A\) is the mixing matrix, \(P\) is a permutation matrix and \(D\) a diagonal matrix with nonzero diagonal entries. If \(\hat{W}A\) converges to the identity matrix, the limiting expected value of \(n(p-1)MD^2\) is the sum of the asymptotic variances of the off-diagonal elements of \(\hat{W}A\). Here \(n\) is the sample size and \(p\) is the number of components.

References

Ilmonen, P., Nevalainen, J. and Oja, H. (2010), Characteristics of multivariate distributions and the invariant coordinate system, Statistics and Probability Letters, 80, 1844--1853.

Ilmonen, P., Nordhausen, K., Oja, H. and Ollila, E. (2010), A New Performance Index for ICA: Properties, Computation and Asymptotic Analysis. In Vigneron, V., Zarzoso, V., Moreau, E., Gribonval, R. and Vincent, E. (editors) Latent Variable Analysis and Signal Separation, 229--236, Springer.

Miettinen, J., Taskinen S., Nordhausen, K. and Oja, H. (2015), Fourth Moments and Independent Component Analysis, Statistical Science, 30, 372--390.

See Also

ASCOV_JADE_est, ASCOV_FOBI_est, JADE, FOBI

Examples

Run this code
# NOT RUN {
# source components have uniform- and exponential(1)- distribution  
fu <- function(x){1/(sqrt(3)*2)}
fe <- function(x){exp(-x-1)}

supp <- matrix(c(-sqrt(3),sqrt(3),-1,Inf), nrow=2, ncol=2, byrow=TRUE)
A<-matrix(rnorm(4),2,2)

jade <- ASCOV_JADE(sdf=c(fu,fe), supp=supp, A=A)
fobi <- ASCOV_FOBI(sdf=c(fu,fe), supp=supp, A=A)

round(jade$COV_W,2)
round(fobi$COV_W,2)

jade$EMD
fobi$EMD
# }

Run the code above in your browser using DataLab