setwd("C:/sepp/work/fabia/experiments") library(fabia) spa <- function(ax,n) { mo <- list() for (j in 1:8) { mo[[j]] <- sum(ax^j)/n } return(mo) } data(Breast_A) X <- as.matrix(XBreast) X <- X- rowMeans(X) XX <- (1/ncol(X))*tcrossprod(X) dXX <- 1/sqrt(diag(XX)+0.001*as.vector(rep(1,nrow(X)))) X <- dXX*X ax <- as.numeric(X) n <- length(ax) spx <- spa(ax,n) exps <- "breast: " write.table(t(as.vector(c(exps,spx))), file = "sparseness.txt", quote = FALSE, sep = "\t",row.names = FALSE,col.names = FALSE,append =FALSE) data(Multi_A) X <- as.matrix(XMulti) X <- X- rowMeans(X) XX <- (1/ncol(X))*tcrossprod(X) dXX <- 1/sqrt(diag(XX)+0.001*as.vector(rep(1,nrow(X)))) X <- dXX*X ax <- as.numeric(X) n <- length(ax) spx <- spa(ax,n) exps <- "multi: " write.table(t(as.vector(c(exps,spx))), file = "sparseness.txt", quote = FALSE, sep = "\t",row.names = FALSE,col.names = FALSE,append =TRUE) data(DLBCL_B) X <- as.matrix(XDLBCL) X <- X- rowMeans(X) XX <- (1/ncol(X))*tcrossprod(X) dXX <- 1/sqrt(diag(XX)+0.001*as.vector(rep(1,nrow(X)))) X <- dXX*X ax <- as.numeric(X) n <- length(ax) spx <- spa(ax,n) exps <- "dlbcl: " write.table(t(as.vector(c(exps,spx))), file = "sparseness.txt", quote = FALSE, sep = "\t",row.names = FALSE,col.names = FALSE,append =TRUE) m <- 100 for (i in 1:m) { load(file=paste("exp_",as.character(i),".RData",sep = "")) X <- dat[[1]] X <- X- rowMeans(X) XX <- (1/ncol(X))*tcrossprod(X) dXX <- 1/sqrt(diag(XX)+0.001*as.vector(rep(1,nrow(X)))) X <- dXX*X ax <- as.numeric(X) n <- length(ax) spx <- spa(ax,n) exps <- paste("exp_",as.character(i),": ",sep = "") write.table(t(as.vector(c(exps,spx))), file = "sparseness.txt", quote = FALSE, sep = "\t",row.names = FALSE,col.names = FALSE,append =TRUE) }