User: Leonard Goldstein
Leonard Goldstein • 80
- Reputation:
- 80
- Status:
- Trusted
- Location:
- United States
- Last seen:
- 3 months, 4 weeks ago
- Joined:
- 4 years, 9 months ago
- Email:
- g*******@gene.com
Posts by Leonard Goldstein
<prev
• 44 results •
page 1 of 5 •
next >
0
votes
1
answer
133
views
1
answers
... 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 3 months ago by
Leonard Goldstein • 80
1
vote
1
answer
133
views
1
answers
... 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 4 months ago by
Leonard Goldstein • 80
0
votes
4
answers
230
views
4
answers
... 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 8 months ago by
Leonard Goldstein • 80
0
votes
4
answers
230
views
4
answers
... 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 8 months ago by
Leonard Goldstein • 80
0
votes
4
answers
230
views
4
answers
... 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 8 months ago by
Leonard Goldstein • 80
0
votes
4
answers
502
views
4
answers
... 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 21 months ago by
Leonard Goldstein • 80
0
votes
4
answers
502
views
4
answers
... 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 21 months ago by
Leonard Goldstein • 80
0
votes
1
answer
231
views
1
answers
... 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 21 months ago by
Leonard Goldstein • 80
1
vote
1
answer
399
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 21 months ago by
Leonard Goldstein • 80
1
vote
3
answers
435
views
3
answers
... 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 22 months ago by
Leonard Goldstein • 80
Latest awards to Leonard Goldstein
No awards yet. Soon to come :-)
Use of this site constitutes acceptance of our User
Agreement
and Privacy
Policy.
Powered by Biostar
version 16.09
Traffic: 352 users visited in the last hour