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]]

