## User: story.benjamin

Reputation:
10
Status:
New User
Location:
Last seen:
22 hours ago
Joined:
2 years, 8 months ago
Email:
s*************@gmail.com

#### Posts by story.benjamin

<prev • 9 results • page 1 of 1 • next >
1
60
views
1
... [Here][1] is a *hacky* fix that I came up with (stealing the code from disjointExons)... but maybe something more robust would be better. [1]: https://www.biostars.org/p/379856/ ...
written 28 days ago by story.benjamin10
1
60
views
1
... Here would be reproducible code: library(GenomicRanges) library(GenomicFeatures) library(GenomicAlignments) hse = makeTxDbFromBiomart( biomart="ensembl", dataset="hsapiens_gene_ensembl", host="grch37.ensembl.org" ) exonicParts = exonicParts( hse, linked.to.single.gene.only ...
written 28 days ago by story.benjamin10
1
60
views
1
... The current DEXSeq manual suggests using: exonicParts = exonicParts( txdb, linked.to.single.gene.only = TRUE ) To generate an "exonic parts" feature set (GRanges object). However, the current exonicParts call in GenomicFeatures omits a step, that was **previously present** for `disjointExon ...
written 28 days ago by story.benjamin10 • updated 27 days ago by Hervé Pagès ♦♦ 14k
1
71
views
1
... I keep all my R packages in a specific location that is not my home directory. For some reason when ExperimentHub was downloading (through biocLite call) it decided to push files to ~/.ExperimentHub... but my server doesn't allow users to surpass a certain amount of space (in their homes) so that p ...
written 10 weeks ago by story.benjamin10 • updated 10 weeks ago by shepherl ♦♦ 1.4k
1
172
views
1
... When I run the program I get the following column names in the output from bam2R    A    T    C    G  -   N  INS DEL HEAD TAIL   QUAL   The "-" represents gaps (so it may be present for spliced alignments; e.g. introns).  My guess is you converted the output to a data.frame which renamed that co ...
written 6 months ago by story.benjamin10
0
212
views
0
... So I think the problem has something to do with the command: exns = stack(exns) which is made by the make_exongr() function, made by make_gbrecord(), which is made by the original readGenBank() command I hacked together a workaround for my purposes but I'm not sure if this is a bug or a 1-off sc ...
written 10 months ago by story.benjamin10
0
212
views
0
... I am trying to pull a specific GenBank entry off of the NCBI database using the genbankr library. My hope was to have all the feature information in a more "R friendly" (e.g. GRanges) format so I could incorporate it into my analysis. Other IDs seem to work but this one was giving me trouble. This ...
written 10 months ago by story.benjamin10
1
553
views
1
... Hi, A colleague approached me in hopes of figuring out a strange export problem they were having with rtracklayer. For some reason it seems that when I have over 255-256 chromosomes (manual tested this) in a genome, that certain chromosome(s) are lost during the export of a bigwig. This was notice ...
written 2.6 years ago by story.benjamin10 • updated 2.6 years ago by Michael Lawrence11k
0
539
views
0
... I am also getting this error for specific datasets. It seems somewhat random. SRP051830 is an example of one that doesn't work.   EDIT: What did work is to set the fileType = 'sra'  ... so perhaps it's an issue with fastq availability. Downside is obviously that you have to do the sra --> fastq ...
written 2.7 years ago by story.benjamin10

#### Latest awards to story.benjamin

No awards yet. Soon to come :-)

Content
Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.