Displaying 3 results from an estimated 3 matches for "logcosh".
Did you mean:
locos
2008 Jul 17
2
fastICA
...repeated runs of fastICA produce quite significantly different
mixing matrices (not only in terms of sign and row order). I'm not a
specialist, so would appreciate any advice on whether this should really be
the case:
> res3 =
> fastICA(af[,2:20],4,alg.typ="parallel",fun="logcosh",alpha=1,method="C",row.norm=TRUE)
colstandard
> res4 =
> fastICA(af[,2:20],4,alg.typ="parallel",fun="logcosh",alpha=1,method="C",row.norm=TRUE)
colstandard
> res3$A
[,1] [,2] [,3] [,4] [,5] [,6]
[,7...
2009 Nov 11
1
Loadings and scores from fastICA?
...tICA.pdf)
if(require(MASS))
{
x <- mvrnorm(n = 1000, mu = c(0, 0), Sigma = matrix(c(10, 3, 3, 1), 2, 2))
x1 <- mvrnorm(n = 1000, mu = c(-1, 2), Sigma = matrix(c(10, 3, 3, 1), 2, 2))
X <- rbind(x, x1)
a <- fastICA(X, 2, alg.typ = "deflation", fun = "logcosh", alpha = 1, method = "R", row.norm = FALSE, maxit = 200, tol = 0.0001, verbose = TRUE)
par(mfrow = c(1, 3))
plot(a$X, main = "Pre-processed data")
plot(a$X%*%a$K, main = "PCA components")
plot(a$S, main = "ICA components")
}
Best...
2009 Nov 04
2
PCA with tow response variables
Hi all,
I'm new to PCA in R, so this might be a basical thing, but I cannot find anything on the net about it.
I need to make a PCA plot with two response variables (df$resp1 and df$resp2) against eight metabolites (df$met1, df$met2, ...) and I don't have a clue how to do... and I've only used the simplest PCAs before, like this:
pcaObj=prcomp(t(df[idx, c(40:47)]))