Hi, I'm running a tutorial ("Meta-analyses of data from two (or more) microarray data sets"), which use wgcna package. I have an error in the function modulePreservation (it is below). I'm using R2.13 Can you help me? Do you know, what is happens?
Thanks Raquel multiExpr = list(A = list(data=t(badea)),B = list(data=t(mayo))) # two independent datasets (dim = 13447 x 36) mp = modulePreservation(multiExpr,multiColor,referenceNetworks=1,verbose=3,networkType="signed",nPermutations=2,maxGoldModuleSize=100,maxModuleSize=400) ..checking data for excessive amounts of missing data.. Flagging genes and samples with too many missing values... ..step 1 Flagging genes and samples with too many missing values... ..step 1 ..unassigned 'module' name: grey ..all network sample 'module' name: gold ..calculating observed preservation values ..calculating permutation Z scores ..Working with set 1 as reference set ....working with set 2 as test set ......working on permutation 1 ......working on permutation 2 Error en zAll[, stat] = z : número de items para para sustituir no es un múltiplo de la longitud del reemplazo (number of items to replace is not a multiple of&nbs p;replacement length) -- ******************************************** Raquel Martinez Garcia, Graduate Student Gastrointestinal Cancer Clinical Research Unit & Structural Biology and BioComputing Programme Spanish National Cancer Research Center (CNIO) Melchor Fernandez Almagro, 3. 28029 Madrid, Spain. Phone: +34 91 732 80 00 #3015 rmartin...@cnio.es http://www.cnio.es ******************************************** **NOTA DE CONFIDENCIALIDAD** Este correo electrónico, y ...{{dropped:7}}
______________________________________________ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.