Search
Question: Limma and DEG
0
gravatar for mathavanbioinfo
23 months ago by
mathavanbioinfo10 wrote:

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")

 

ADD COMMENTlink modified 23 months ago by Aaron Lun21k • written 23 months ago by mathavanbioinfo10
1
gravatar for Aaron Lun
23 months ago by
Aaron Lun21k
Cambridge, United Kingdom
Aaron Lun21k wrote:

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 COMMENTlink modified 23 months ago • written 23 months ago by Aaron Lun21k
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.2.0
Traffic: 373 users visited in the last hour