User: Leonard Goldstein

Reputation:
80
Status:
Trusted
Location:
United States
Last seen:
1 month ago
Joined:
4 years, 6 months ago
Email:
g*******@gene.com

Posts by Leonard Goldstein

<prev • 44 results • page 1 of 5 • next >
1
79
views
1
... 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 4 weeks ago by Leonard Goldstein80
1
79
views
1
... 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 weeks ago by Leonard Goldstein80
4
189
views
4
... 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 5 months ago by Leonard Goldstein80
4
189
views
4
... 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 months ago by Leonard Goldstein80
4
189
views
4
... 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 months ago by Leonard Goldstein80
4
454
views
4
... 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 18 months ago by Leonard Goldstein80
4
454
views
4
... 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 18 months ago by Leonard Goldstein80
1
204
views
1
... 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 18 months ago by Leonard Goldstein80
1
351
views
1
... 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 18 months ago by Leonard Goldstein80
3
388
views
3
... 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 19 months ago by Leonard Goldstein80

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.