info192
Messages postés9Date d'inscriptionvendredi 1 février 2013StatutMembreDernière intervention12 décembre 2013
-
Modifié par pijaku le 13/12/2013 à 09:57
Bonsoir,
S'il vous plaît qui peut m'expliquer les instructions de cet exemple bicbin.R
#library(SparseM) #mat <- matrix(sample(c(rep.int(0, 100), 1), 800*25000, replace=TRUE), 800, 25000) #mat <- as.matrix.csr(mat) #mat[10:15, 10:15] <- 1 #M <- nrow(mat) # TFs #N <- ncol(mat) # Genes #alpha <- 0.5 #beta <- 0.5 #p <- sum(attr(mat, "ra")) / (M*N) # #max <- 0; #for (i in 1:300) { # print(i) # res <- BicBin(mat, alpha, beta, p, proc_genes=TRUE) # if (res > max) { max <- res } # res <- BicBin(mat, alpha, beta, p, proc_genes=FALSE) # if (res > max) { max <- res } #}