Limma and DEG
1
0
Entering edit mode
@mathavanbioinfo-11955
Last seen 2 days 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 • 805 views
ADD COMMENT
1
Entering edit mode
Aaron Lun ★ 28k
@alun
Last seen 1 hour 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: 891 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