Search
Question: Limma and DEG
0
gravatar for mathavanbioinfo
14 months ago by
mathavanbioinfo0 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 14 months ago by Aaron Lun18k • written 14 months ago by mathavanbioinfo0
1
gravatar for Aaron Lun
14 months ago by
Aaron Lun18k
Cambridge, United Kingdom
Aaron Lun18k 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 14 months ago • written 14 months ago by Aaron Lun18k
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: 247 users visited in the last hour