Limma and DEG
1
0
Entering edit mode
@mathavanbioinfo-11955
Last seen 12 weeks ago
India

hi

I am using LIMMA to identify DEGs. I got the top 5 DEGs. But i am not able to get the whole list of DEGs with proper id and genes. 

Here i used the code 

write.table(topTable(fit, coef=1, adjust="fdr", sort.by="B", number=500), file="limma_complete.xls", row.names=F, sep="\t")

 

limma • 845 views
ADD COMMENT
1
Entering edit mode
Aaron Lun ★ 28k
@alun
Last seen 7 hours ago
The city by the bay

If you want the entire list of DEGs, set number=Inf.

If you want the gene names in the output, supply a data frame or vector with the relevant information to genelist when calling topTable. Or, store it in fit$genes, as shown in the example below:

fit$genes <- data.frame(GeneId=paste0("Gene", 1:nrow(fit)))
ADD COMMENT

Login before adding your answer.

Traffic: 1388 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