User: Leonard Goldstein

Reputation:
70
Status:
Trusted
Location:
United States
Last seen:
3 weeks, 2 days ago
Joined:
4 years, 2 months ago
Email:
g*******@gene.com

Posts by Leonard Goldstein

<prev • 42 results • page 1 of 5 • next >
0
votes
4
answers
104
views
4
answers
Comment: C: Which sequencing parameters are required for SGSeq analysis of splice events?
... There is no specific requirement for fragment or read length. I would follow standard recommendations for library preparation. The sequencing specs and read length sound fine. ...
written 23 days ago by Leonard Goldstein70
0
votes
4
answers
104
views
4
answers
Answer: A: Which sequencing parameters are required for SGSeq analysis of splice events?
... Hi Andrea,  If you have the option to do strand-specific sequencing, this is still a good idea for transcriptome analysis (just not required for SGSeq). Read depth and replicates sound reasonable.  Best,  Leonard ...
written 5 weeks ago by Leonard Goldstein70
0
votes
4
answers
104
views
4
answers
Answer: A: Which sequencing parameters are required for SGSeq analysis of splice events?
... Hi Andrea,  Actually strand-specific sequencing is not a requirement (the package was written a few years ago when most available data sets were not strand-specific). Since SGSeq relies on split reads, longer reads are helpful. There is no requirement for a particular sequencing depth. More reads a ...
written 5 weeks ago by Leonard Goldstein70
0
votes
4
answers
373
views
4
answers
Answer: A: Cannot extract TxFeatures for a defined GRanges object
... exonsBy seems to be dropping a lot of genes, including MYC: > tx <- makeTxDbFromBiomart() > txs_by_gene <- transcriptsBy(tx, "gene") > exs_by_gene <- exonsBy(tx, "gene") > length(txs_by_gene) [1] 63967 > length(exs_by_gene) [1] 36751 > subsetByOverlaps(txs_by_gene, GRang ...
written 13 months ago by Leonard Goldstein70
0
votes
4
answers
373
views
4
answers
Answer: A: Cannot extract TxFeatures for a defined GRanges object
... It looks like MYC is missing from the TxDb object you created, so this is a question about the makeTxDbFromBiomart function in the GenomicFeatures package. Maybe repost your question and add GenomicFeatures as one of the tags, and hopefully someone can help. Alternatively you can try using one of th ...
written 13 months ago by Leonard Goldstein70
0
votes
1
answers
160
views
1
answers
Answer: A: SGFeatureCounts function in SGSeq package failed to generate splice graphs
... Hi Magda, Can you please provide more details -- what commands were used to generate the SGFeatureCounts object and what package versions were used i.e. output of sessionInfo(). For larger analyses I strongly recommend not to use wrapper functions like analyzeFeatures() and instead run individual a ...
written 13 months ago by Leonard Goldstein70
0
votes
1
answers
266
views
1
answers
Answer: A: SGSeq error on getBamInfo
... The error message should list the problematic files but it returns an empty string. Looks like something went wrong that was not caught by the internal checks (probably unrelated to the XS tag). Are your BAM file paths correct and are index files present? You can try running the following and see if ...
written 14 months ago by Leonard Goldstein70
0
votes
3
answers
306
views
3
answers
Answer: A: R session crash with SGSeq::convertToTxFeatures
... Hi, since I can't reproduce the problem it is difficult to troubleshoot. As an alternative to TxDb objects, the convertToTxFeatures function also accepts a GRangesList of exons grouped by transcript. So you could try something like this > tx <- exonsBy(hsaEnsembl, "tx", use.names = TRUE) > ...
written 14 months ago by Leonard Goldstein70
0
votes
3
answers
306
views
3
answers
Answer: A: R session crash with SGSeq::convertToTxFeatures
... I could run this without any issues. If the problem persists maybe you can try reinstalling R and Bioconductor?  > library(GenomicFeatures) > library(SGSeq) > hsaEnsembl <- makeTxDbFromBiomart(dataset="hsapiens_gene_ensembl") > txf <- convertToTxFeatures(hsaEnsembl) > > se ...
written 14 months ago by Leonard Goldstein70
0
votes
1
answers
486
views
1
answers
Comment: C: SGSeq: moving toward diffex from SGSeq analysis
... 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 ...
written 15 months ago by Leonard Goldstein70

Latest awards to Leonard Goldstein

No awards yet. Soon to come :-)

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.2.0
Traffic: 233 users visited in the last hour