Got 24 samples, and got: 25 data fields. Header: uninfected_01 infected_01 infected_02 infected_03 infected_04 infected_05 infected_06 uninfected_02 uninfected_03 infected_07 infected_08 infected_09 infected_10 infected_11 infected_12 infected_13 uninfected_04 uninfected_05 infected_14 infected_15 infected_16 infected_17 infected_18 uninfected_06 Next: TRINITY_DN0_c0_g1 1.87 2.00 311.14 1749.31 1307.99 5342.68 1991.08 0.00 0.00 0.00 0.00 0.00 3.39 1.89 2.47 0.00 0.00 0.00 0.00 2.63 0.00 6.13 1.91 1.50 $VAR1 = { 'uninfected_06' => 24, 'infected_16' => 21, 'infected_11' => 14, 'infected_17' => 22, 'infected_06' => 7, 'infected_02' => 3, 'infected_10' => 13, 'infected_08' => 11, 'infected_07' => 10, 'uninfected_01' => 1, 'infected_05' => 6, 'infected_04' => 5, 'infected_15' => 20, 'infected_03' => 4, 'uninfected_03' => 9, 'infected_09' => 12, 'infected_18' => 23, 'infected_12' => 15, 'infected_14' => 19, 'uninfected_05' => 18, 'uninfected_02' => 8, 'infected_01' => 2, 'uninfected_04' => 17, 'infected_13' => 16 }; Contrasts to perform are: $VAR1 = [ [ 'infected', 'uninfected' ] ]; Loading required package: edgeR Loading required package: limma