User: Leonard Goldstein
Leonard Goldstein • 110
- Reputation:
- 110
- Status:
- Trusted
- Location:
- United States
- Last seen:
- 1 week, 1 day ago
- Joined:
- 5 years, 2 months ago
- Email:
- g*******@gene.com
Posts by Leonard Goldstein
<prev
• 49 results •
page 1 of 5 •
next >
0
votes
1
answer
44
views
1
answers
... Please read the help page for getBamInfo and SGSeq vignette section 2. Note that your BAM files must contain the custom tag 'XS' for spliced reads. If you are looking for help with a particular problem or error you need to provide more information (error message, sessionInfo() output, aligner inform ...
written 8 days ago by
Leonard Goldstein • 110
1
vote
1
answer
53
views
1
answers
Comment:
C: Custom gene annotations for SGSeq
... You can add junctions manually for a TxFeatures object. But you have to be careful that coordinates are consistent with exons. Junction starts and ends have to coincide with exon ends and start, respectively. Here is an example for TxFeatures object txf_ann that comes with the package:
new_junc ...
written 16 days ago by
Leonard Goldstein • 110
2
votes
1
answer
53
views
1
answers
... Yes you can work with a combination of annotated and predicted gene models. The easiest way to do this is to use analyzeFeatures() for annotated gene models as explained in the vignette but with predict = TRUE. ...
written 16 days ago by
Leonard Goldstein • 110
0
votes
1
answer
76
views
1
answers
... Actually looking at the code in more detail, BamFile details from the sample information currently do not propagate. The hard-coded behavior is that only primary alignments in the BAM file are considered for analysis. In other words, each read should be analyzed only once, but the analysis may inclu ...
written 8 weeks ago by
Leonard Goldstein • 110
0
votes
1
answer
76
views
1
answers
... You can influence which reads are analyzed when setting up the data frame with sample information. To do this column file_bam needs to be a BamFileList and sample information must be stored in a DataFrame object. Please see the manual page for Rsamtools::BamFileList for more details. ...
written 8 weeks ago by
Leonard Goldstein • 110
0
votes
1
answer
192
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 8 months ago by
Leonard Goldstein • 110
1
vote
1
answer
192
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 8 months ago by
Leonard Goldstein • 110
0
votes
4
answers
311
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 13 months ago by
Leonard Goldstein • 110
0
votes
4
answers
311
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 13 months ago by
Leonard Goldstein • 110
0
votes
4
answers
311
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 13 months ago by
Leonard Goldstein • 110
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: 147 users visited in the last hour