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