Question: How to extract genes with greatest BCV?
0
gravatar for Mich
12 months ago by
Mich0
Singapore
Mich0 wrote:

Hi,

I am analyzing RNAseq data using the EdgeR package and have made the MDS plot using the following code.

y<-DGEList(counts=data, group=Source)

keep<- rowSums(cpm(y)>1) >= 1

y<- y[keep, , keep.lib.sizes=FALSE]

y<-calcNormFactors(y)

plotMDS(y, method="bcv",dim= c(1,2), top=1000)

I am wondering if it is possible to obtain the gene names of the top 1000 genes with the greatest biological coefficient of variation.

Many thanks! 

edger • 226 views
ADD COMMENTlink modified 12 months ago by Aaron Lun24k • written 12 months ago by Mich0
Answer: How to extract genes with greatest BCV?
1
gravatar for Aaron Lun
12 months ago by
Aaron Lun24k
Cambridge, United Kingdom
Aaron Lun24k wrote:

The BCV is just the square root of the tagwise dispersion, so just take the top 1000 genes with the largest tagwise dispersions:

order(y$tagwise.dispersion, decreasing=TRUE)[1:1000]

... after you've run estimateDisp.

ADD COMMENTlink written 12 months ago by Aaron Lun24k

Thanks Aaron!

ADD REPLYlink written 12 months ago by Mich0
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 16.09
Traffic: 305 users visited in the last hour