Hello Guys, I'm using the package RankProduct to determine the differential expression of genes in 2 classes. The importance of the thing is to know: - How can I determine which class is the class of up-regulated genes and which to the down regulated gene - How can i interpret the results : $pfp, $RPs, $RPrank, Orirank, $Orirank$"class1 > class 2" $Orirank$"class1 < class 2" $AveFC Table1: Genes called significant under class1 < class2 Table2: Genes called significant under class1 > class2 Thanks for your help. Example: Formula: RP.out<-RP(data, pt.cl, g.names, num.perm = 100, logged = TRUE, na.rm = FALSE, plot = FALSE, rand =123) $pfp class1 < class2 class1 > class 2 1 0.97115385 0.002000000 $RPs class1 < class2 class1 > class 2 1 37.54242 21.63239 $RPrank class1 < class2 class1 > class 2 1 26.0 5.0 Orirank $Orirank$"class1 < class2" [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [,11] [,12] …………………………………………….. $Orirank$"class1 > class 2" …………………………………….. $AveFC log/unlog(class1/class2) 1 0.0694052288 ……………………………… #Genes by controlling pfp <0.05 topGene(RP.out, cutoff = 0.05, logged = TRUE, logbase = 2, g.names) Table1: Genes called significant under class1 < class2 Table2: Genes called significant under class1 > class2 $Table1 gene.index RP/Rsum FC:(class1/class2) pfp Gen n 51 15.8595 0.8204 0.0000 Gen n-1 65 18.3764 0.9582 0.0250 Gen n-2 75 18.6893 0.8916 0.0167 Gen n-3 52 19.0778 0.8928 0.0150 ……………………………………………. $Table2 gene.index RP/Rsum FC:(class1/class2) pfp Gen n-3 98 10.0913 1.4331 0.0000 Gen n-2 39 15.2657 1.4062 0.0000 Gen n-1 2 19.8796 1.1414 0.0000 Gen n 97 20.3225 1.1446 0.0000 ………………………………………………... --------------------------------- [[alternative HTML version deleted]]