Hi,
I am doing goseq on RNA-seq data of S. pombe. Below I am showing the last few lines of calculating the enriched GOs for my up-regulated genes. In the GO.hyper_up I don't get two last columns "term" and "ontology" as it is shown in goseq manual. Does any one have any suggestions? Let me know if you need more info
I am doing goseq on RNA-seq data of S. pombe. Below I am showing the last few lines of calculating the enriched GOs for my up-regulated genes. In the GO.hyper_up I don't get two last columns "term" and "ontology" as it is shown in goseq manual. Does any one have any suggestions? Let me know if you need more info
Code:
> Spo.genes.Go = read.table('gene_association.pombase', sep='\t', comment='!', stringsAsFactors = F, quote="\"") > gene2GO = data.frame(Spo.genes.Go$V2, Spo.genes.Go$V5) > colnames(gene2GO)=c("id", "GO") > head(gene2GO) id GO 1 SPAC212.11 GO:0005634 2 SPAC212.11 GO:0043140 3 SPAC212.11 GO:0000722 4 SPAC212.08c GO:0031362 5 SPAC212.04c GO:0009897 6 SPAC212.01c GO:0009897 > GO.hyper_up=goseq(pwf_up, gene2cat = gene2GO, method="Hypergeometric")