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")
If you want the entire list of DEGs, set number=Inf.
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:
genelist
topTable
fit$genes
fit$genes <- data.frame(GeneId=paste0("Gene", 1:nrow(fit)))
Login before adding your answer.
Use of this site constitutes acceptance of our User Agreement and Privacy Policy.