From the two t-tests: > fc.Th[ order(t.pval.Th) ] [1:10] SLC7A5 KIF11 C1orf29 KIAA0101 CDC20 GARNL4 MICB -1.1060657 0.6082505 3.6931245 1.4301633 0.4454738 0.4562688 0.6838960 RRM2 OAS1 OAS1 1.4466054 1.6530430 1.9240594 > fc.CTL[ order(t.pval.CTL) ] [1:10] DUSP8 PIGF GTF2E1 IFI16 CCR5 CHI3L2 C10orf6 -0.9767053 0.6449114 0.8589276 1.1565018 2.2925576 1.7866050 0.2963409 SLC7A5 TXNDC9 HIMAP4 -1.3253478 1.0762418 1.6573616 From the two-way ANOVA: > names( HIV.pval[ order(HIV.pval) ] ) [1:20] [1] "SLC7A5" "KIAA0101" "OAS1" "HIMAP4" "OAS1" "CCR5" [7] "IFI16" "IFI16" "RRM2" "C1orf29" "AIM2" "CGI-01" [13] "IFI44" "RRM2" "TYMS" "MX1" "C3orf4" "MGC2477" [19] "SMAD7" "PDE4B" > names( CellType.pval[ order(CellType.pval) ] ) [1:20] [1] "CD8A" "CD8B1" "KLRK1" "CD8B1" "CRTAM" "C3F" [7] "HIC" "BZRAP1" "KIAA0182" "PECAM1" "CGI-62" "GPA33" [13] "ChGn" "PCNT2" "TIAM1" "CSDA" "KLRD1" "XCL1" [19] "C9orf28" "KLRD1" > names( Interaction.pval[ order(Interaction.pval) ] ) [1:10] [1] "CHI3L2" "CCR5" "ADCY9" "KIF14" "LOC146909" "PCAF" [7] "ZBTB20" "SDCCAG8" "MAP3K7IP2" "ZNF174"