[BioC] Mfuzz package

dorothyc dorothyc at bcgsc.ca
Thu Apr 2 23:06:18 CEST 2009


Hi,

I'm trying out the Mfuzz package for soft clustering.  I've  
encountered an error (please see below).  I'd be really grateful if  
you could shed some light.

 > filename <- my_table
 > my.data <- read.table(filename, sep="\t")
 > set <-new("ExpressionSet", exprs = as.matrix(my.data));
 > eset.r <- filter.NA(eset, thres = 0.25)
4 genes excluded.
 > eset.f <- fill.NA(eset.r, mode = "mean")
 > tmp <- filter.std(eset.f, min.std = 0)
1 genes excluded.
 > eset.s <- standardise(eset.f)
 > cl <- mfuzz(eset.s, c=6, m=2)
Error in sample(1:xrows, ncenters) :
   cannot take a sample larger than the population when 'replace =  
FALSE'

Thanks,
Dorothy



More information about the Bioconductor mailing list