[Questions about the results of the enrichGo()]More genes are involved but less biological activity.
0
0
Entering edit mode
414457044 • 0
@414457044-15408
Last seen 6.7 years ago

Hi:

I downloaded the gene mutation table for breast cancer patients and ranked them according to the mutation frequency of the genes.

The first 885 genes in the database (with a mutation frequency greater than or equal to 3) were analyzed with enrichGO() to obtain 144 results for participating in biological activities.

The first 348 genes in the database (with a mutation frequency greater than or equal to 4) were analyzed with enrichGO() to obtain the results of 173 biological activities.

But the first 885 genes contained the first 348 genes. Why did they get less results than using the first 348 genes?

The following code:
```
gene <- read.table("morethan3.txt",colClasses = "character")
gene<-as.character(gene[,1])
ego <- enrichGO(gene = gene, OrgDb = org.Hs.eg.db,
                ont = "BP", pAdjustMethod = "BH",
                pvalueCutoff = 0.05, qvalueCutoff = 0.05)
ego_10<-as.data.frame(ego)
write.csv(ego_10,file = "./ego_3_0.05_0.05.csv", row.names = F)

gene <- read.table("morethan4.txt",colClasses = "character")
gene<-as.character(gene[,1])
ego <- enrichGO(gene = gene, OrgDb = org.Hs.eg.db,
                ont = "BP", pAdjustMethod = "BH",
                pvalueCutoff = 0.05, qvalueCutoff = 0.05)
ego_10<-as.data.frame(ego)
write.csv(ego_10,file = "./ego_4_0.05_0.05.csv", row.names = F)
```

in  file morethan3.txt,the data is:
[morethan3.txt](https://github.com/chriszhuge/enrich/blob/master/morethan3.txt)
and the result is :
[ego_3_0.05_0.05.csv](https://github.com/chriszhuge/enrich/blob/master/ego_3_0.05_0.05.csv)

in  file morethan4.txt,the data is:
[morethan4.txt](https://github.com/chriszhuge/enrich/blob/master/morethan4.txt)
and the result is :
[ego_4_0.05_0.05.csv](https://github.com/chriszhuge/enrich/blob/master/ego_4_0.05_0.05.csv)

enrichGO() results • 900 views
ADD COMMENT

Login before adding your answer.

Traffic: 464 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6