User: Leonard Goldstein

gravatar for Leonard Goldstein
Reputation:
80
Status:
Trusted
Location:
United States
Last seen:
6 months, 3 weeks ago
Joined:
5 years ago
Email:
g*******@gene.com

Posts by Leonard Goldstein

<prev • 44 results • page 1 of 5 • next >
0
votes
1
answer
165
views
1
answers
Comment: C: I can't import my gff file by `importTransprits()`,when I worked with SGSeq.
... Both are possible. Please see vignette section 6 for an example on how to analyze a particular region of the genome using the `which` argument. If no region is specified the analysis is genome-wide. Depending on your data set, genome-wide predictions can be computationally intensive. In case you run ...
written 6 months ago by Leonard Goldstein • 80
1
vote
1
answer
165
views
1
answers
Answer: A: I can't import my gff file by `importTransprits()`,when I worked with SGSeq.
... When posting a question about a software package, please always tag your post with the package name. This triggers an automatic email from the system to the package maintainer. Regarding your question -- the `tag_tx` and `tag_gene` arguments are for specifying names of the relevant GFF tags. Pleas ...
written 6 months ago by Leonard Goldstein • 80
0
votes
4
answers
278
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 11 months ago by Leonard Goldstein • 80
0
votes
4
answers
278
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 11 months ago by Leonard Goldstein • 80
0
votes
4
answers
278
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 11 months ago by Leonard Goldstein • 80
0
votes
4
answers
531
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 23 months ago by Leonard Goldstein • 80
0
votes
4
answers
531
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 23 months ago by Leonard Goldstein • 80
0
votes
1
answer
287
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 23 months ago by Leonard Goldstein • 80
1
vote
1
answer
446
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 2.0 years ago by Leonard Goldstein • 80
1
vote
3
answers
482
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 2.0 years ago by Leonard Goldstein • 80

Latest awards to Leonard Goldstein

No awards yet. Soon to come :-)
Content
Help
Access
Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 16.09
Traffic: 274 users visited in the last hour