Hi , I have tried to run GOseq on my differential expressed gene list, but when I try to run the following command: GO.wall <-goseq(DEGgenes, pwf, "bosTau4", "ensGene", method= "Wallenius") I get the following error: Error in pwf$DEgenes : $ operator is invalid for atomic vectors Does anyone has experience in solving this issue? Thanks, Fernando