Dear users,

I use the LPE package to detect differentially expressed genes but I'm pretty surprised that I didn't find the same genes according the order of my groups

>library(LPE)
>ind_g1<-c(2,4:6)
>ind_g2<-c(1,3,7)
>group1<-data[,g1]
>group2<-data[,g2]
# Finding the baseline distribution of condition 1 and 2.
>var.1 <- baseOlig.error(group1, q=0.01)
>var.2 <- baseOlig.error(group2, q=0.01)
>lpe.result <- lpe(group1,group2, var.1, var.2,probe.set.name=rownames(data))
>finalresult<-fdr.adjust(lpe.result,adjp="resamp")
>finalresult
      target.fdr z.critical
 [1,]      0.001   5.606208
 [2,]      0.010   5.606208
 [3,]      0.020   5.606208
 [4,]      0.030   5.606208
 [5,]      0.040   4.720516
 [6,]      0.050   4.123739
 [7,]      0.060   2.823956
 [8,]      0.070   2.593902
 [9,]      0.080   2.534597
[10,]      0.090   2.482537
[11,]      0.100   2.375929
[12,]      0.150   2.015756
[13,]      0.200   1.630128
[14,]      0.500   1.500434

>lpe.inv<-lpe(group2,group1,var.2,var.1,probe.set.name=rownames(data))
>final.inv<- fdr.adjust(lpe.inv, adjp="resamp")
>final.inv
      target.fdr z.critical
 [1,]      0.001   3.096275
 [2,]      0.010   2.793027
 [3,]      0.020   2.593902
 [4,]      0.030   2.482537
 [5,]      0.040   2.482537
 [6,]      0.050   2.368713
 [7,]      0.060   2.368713
 [8,]      0.070   2.216946
 [9,]      0.080   2.091183
[10,]      0.090   1.774938
[11,]      0.100   1.639910
[12,]      0.150   1.547681
[13,]      0.200   1.537522
[14,]      0.500   1.459746

When I change the order of my groups, I didn't find the same z-critical. Do you think that's right?

Thanks

Sabrina


Sabrina Carpentier
Service Bioinformatique
Institut Curie - Bat. Trouillet Rossignol (4e étage)
26 rue d'Ulm - 75248 Paris Cedex 5 - FRANCE
Tel : +33 1 42 34 65 21 
	[[alternative HTML version deleted]]

