tests/one-component.R

## failed in 1.1-9
library(fastICA)
set.seed(49394)
S <- matrix(runif(10000), 5000, 2)
A <- matrix(c(1, 1, -1, 3), 2, 2, byrow = TRUE)
X <- S%*%A
a <- fastICA(X, 1, alg.typ = "parallel", fun = "logcosh", alpha = 1,
             method = "C", row.norm = FALSE, maxit = 200,
             tol = 0.0001, verbose = TRUE) # good
a <- fastICA(X, 1, alg.typ = "parallel", fun = "logcosh", alpha = 1,
             method = "R", row.norm = FALSE, maxit = 200,
             tol = 0.0001, verbose = TRUE) # fails

Try the fastICA package in your browser

Any scripts or data that you put into this service are public.

fastICA documentation built on Nov. 27, 2023, 9:06 a.m.