Hi:

I am trying to run globaltest with an affy data set containing 30 
mouse430a2 chips. These are the codes I am using:
 > library(globaltest)
 > library(mouse430a2)
 > kegg <- as.list(mouse430a2PATH2PROBE)
# treatment is the factor indicating different treatments (containing 3 
time points which didn't seperate here).
 > treatment
 [1] Control Control Control Control CD137   CD137   UC10    UC10    
Control Control Control Control CD137   CD137   CD137 
[16] CD137   UC10    UC10    UC10    UC10    Control Control Control 
Control CD137   CD137   UC10    UC10    UC10    UC10  
Levels: Control CD137 UC10
 > gt.kegg.all <- globaltest(as.matrix(rma.data), treatment, kegg)

Everythig was fine when I ran all the dataset together as above. However 
when I try to use subset data and rerun the globaltest, I got following 
error:
 > gt.kegg.day8 <- globaltest(as.matrix(rma.data[ , 9:20]), 
treatment[9:20], kegg)
Error in cbind(deparse.level, ...) : number of rows of matrices must 
match (see arg 2)

I've wasted a lot of time trying to fix this error with no luck. Please 
help!

Thanks,

Xia



	[[alternative HTML version deleted]]

