Question: SGSeq: moving toward diffex from SGSeq analysis
0
gravatar for Sylvain Foisy
2.0 years ago by
Canada
Sylvain Foisy30 wrote:

Hi,

I have use SGSeq mostly to look at exon usage for genes of interest in my RNASeq datasets. I would like to move toward exon diffex analysis. I followed the recipes in the HTML guide to get counts, exons and gene ID but I am perplexed about what I got... I have an SGFeatureCounts object with created from rowRanges on the SGVariantCounts object previously created that has all the infos like geneName, etc. However, after performing the getSGVariantCounts method on it and using the eventID method, I only get integers, that I am guessing are the positions of the genes in the previous object. How can I get geneName instead or possibly add geneName data on the objects?

I am guessing that it would always be possible to go back with integer values to explore the SGFeatureCounts object?

Thanks in advance

S

rnaseq sgseq • 587 views
ADD COMMENTlink modified 2.0 years ago by Leonard Goldstein80 • written 2.0 years ago by Sylvain Foisy30
Answer: SGSeq: moving toward diffex from SGSeq analysis
0
gravatar for Leonard Goldstein
2.0 years ago by
United States
Leonard Goldstein80 wrote:

Hi Sylvain, 

Not sure I completely follow your question. You seem to be confusing SGFeatureCounts and SGVariantCounts objects. The idea behind section 11 of the vignette is to test for differential splice variant usage within a splice event, rather than differential exon usage within a gene. This approach will detect differential splice events (which are identified by event IDs) and you can look them up in your original SGVariants/Counts object. 

Leonard

ADD COMMENTlink written 2.0 years ago by Leonard Goldstein80

Hi Leonard,

No problem ;-) Ok, as per section 11, I created objects for counts, exon ID and gene ID with the methods in the vignette:

# My R stuff:
# sgfc.denovo.var is my equivalent to sgvc_pred
> sgv.ranges<-rowRanges(sgfc.denovo.var)
> sgv.counts<-getSGVariantCounts(sgv.ranges,sample_info=data.r)
> x <-counts(sgv.counts)
> vid<-variantID(sgv.counts> 
> eid<-eventID(sgv.counts)

When I look at vid or eid, I only get matrices of integer values, so I am guessing that they are corresponding indexes in the previous  sgv.ranges object. Am I right? If this is the case, what is the easiest way to make out the gene symbol from the geneName column in sgv.ranges?

Thanks in advance

S

ADD REPLYlink modified 24 months ago • written 24 months ago by Sylvain Foisy30

Hi Sylvain, 

Yes variantID() and eventID() are accessor functions for metadata columns in the SGVariant/Counts objects. You can see all available columns by typing mcols(sgv.counts). To extract gene names you can simply use geneNames(sgv.counts). Note this will return an empty CharacterList unless you previously annotated your predictions. I hope this helps. 

Leonard

ADD REPLYlink written 24 months ago by Leonard Goldstein80

Hi Leonard,

Thanks for the info; I'll be playing some more and if needed, may be bug you again ;-)

Best regards

S

ADD REPLYlink written 24 months ago by Sylvain Foisy30
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: 127 users visited in the last hour