Nothing
FOBI <-
function(X, na.action = na.fail)
{
X <- na.action(X)
X <- as.matrix(X)
n <- nrow(X)
p <- ncol(X)
Col.center <- colMeans(X)
X <- sweep(X, 2, Col.center, "-")
COV <- crossprod(X)/n
EVD.COV <- eigen(COV, symmetric=TRUE)
COV.inv.sqrt <- EVD.COV$vectors %*% tcrossprod(diag((1/EVD.COV$values)^0.5),EVD.COV$vectors)
Y <- tcrossprod(X, COV.inv.sqrt)
r <- sqrt(rowSums(Y^2))
Y <- r * Y
COV4 <- crossprod(Y)/(n*(p+2))
EVD.COV4 <- eigen(COV4, symmetric=TRUE)
W <- crossprod(EVD.COV4$vectors, COV.inv.sqrt)
#W <- diag(sign(rowMeans(W)))%*%W
W <- sweep(W, 1, sign(rowMeans(W)), "*")
S <- tcrossprod(X, W)
colnames(S) <- paste("IC.", 1:p, sep="")
res <- list(W=W, EV=EVD.COV4$values, Xmu=Col.center, S=S)
class(res) <- "bss"
return(res)
}
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.