Hi all, I have expression data from 3 different laboratories i wanted perform meta analysis using rank product package i wrote following code having 21 samples library(RankProd) #loading normalized data si<-read.table("affy_all.txt",sep="\t",header=TRUE,as.is=T) #reading affy ids si.gnames<-si[,1] #Creating class variables si1<-7 si2<-14 si.cl<-rep(c(0,1),c(si1,si2)) 0 0 0 0 0 0 0 1 1 1 1 1 1 1 1 1 1 1 1 1 1 #Creating origin variable si.origin<-c(rep(1,6),rep(2,3),rep(3,12)) 1 1 1 1 1 1 2 2 2 3 3 3 3 3 3 3 3 3 3 3 3 #Identifying regulated genes RP.adv.out<-RPadvance(si,si.cl,si.origin,num.perm=100,na.rm=FALSE,gene .names =si.gnames,rand=123) *But it gives error :( as follows* *Error in RPadvanceV1(data, cl, origin, num.perm, logged, na.rm, gene.names, : * * Number of classes and/or origins should match the columns in the data* I checked my data it is correct ...couldnt figure out problem Can any any body help Thanks in Advance, Sainitin -- Sainitin D PhD student Biotechnology Center Technische Universität Dresden Tatzberg 47/49 01307 Dresden, Germany [[alternative HTML version deleted]]