Question: How do I get groseq signal for each gene from groSeq in groHMM
0
gravatar for jsalsager
12 months ago by
jsalsager40
jsalsager40 wrote:

hello!

I worked my way through the tutorial for groHMM but I'm not sure how to get signal (RPKM, FPKM etc.) for each gene.  The tutorial just shows how to make graphs.

Thanks!

grohmm groseq • 276 views
ADD COMMENTlink modified 12 months ago by anusha.nagari30 • written 12 months ago by jsalsager40
Answer: How do I get groseq signal for each gene from groSeq in groHMM
1
gravatar for anusha.nagari
12 months ago by
UT SOUTHWESTERN MEDICAL CENTER
anusha.nagari30 wrote:

Hello,

You can do that using the following extension to section 3.8 in the tutorial:

d <- calcNormFactors(d)
d_rpkm <- data.frame(rpkm(d,width(kgChr7),normalized.lib.sizes=TRUE))
gene_rpkm <- data.frame(kgChr7,kgLimit, d_rpkm)
head(gene_rpkm)

 

 

ADD COMMENTlink written 12 months ago by anusha.nagari30

Thank you again, I really appreciate your help!

ADD REPLYlink written 12 months ago by jsalsager40

Can I bug you with one final question?

If I wanted to export a similar file but using groHMM called transcripts rather than knownGenes, how would I change that code? I tried switching kgLimit with txLimit and ended up getting errors that I was combining a date.frame with different lengths.

Thanks so much for all your help!

Edit:

Would this work?

d_rpkm <- data.frame(rpkm(d,width(txFinal),normalized.lib.sizes=TRUE))
transcript_rpkm <- data.frame(txFinal, d_rpkm)
head(transcript_rpkm)
ADD REPLYlink modified 12 months ago • written 12 months ago by jsalsager40
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: 180 users visited in the last hour