User: benformatics

gravatar for benformatics
benformatics150
Reputation:
150
Status:
Trusted
Location:
Last seen:
5 days, 5 hours ago
Joined:
1 year, 9 months ago
Email:
s*************@gmail.com

Posts by benformatics

<prev • 21 results • page 1 of 3 • next >
1
vote
1
answer
149
views
1
answers
Answer: A: get mRNAs by TxDb
... The *transcripts()* function will return coordinates for the entire range of a given transcript including introns. To get only the exonic regions use *exonsBy()* ## known exons for each transcript exonsBy(TxDb.Hsapiens.UCSC.hg38.knownGene,'tx') ## known exons for each gene exon ...
written 11 weeks ago by benformatics150
0
votes
2
answers
1.9k
views
2
answers
Comment: C: COSMIC vcf file compatibility for Mutect2
... If anyone wants to do this in R, it a lot easier: library(VariantAnnotation) ## import cosmic vcf file cos.ens <- readVcf('./data/CosmicCodingMuts.vcf') ## check current chr names seqinfo(cos.ens) ## set UCSC compatible chromosome names seqlevelsStyle(cos.e ...
written 3 months ago by benformatics150
1
vote
1
answer
269
views
1
answers
Comment: C: smallRNA aligner and pipeline
... BWA Small RNA: https://www.biostars.org/p/3399/ There are a lot of small RNA profiling papers here are some examples: 1. https://www.ncbi.nlm.nih.gov/pmc/articles/PMC3334700/ 2. https://www.ncbi.nlm.nih.gov/pmc/articles/PMC4079977/ ...
written 8 months ago by benformatics150
0
votes
1
answer
269
views
1
answers
Answer: A: smallRNA aligner and pipeline
... Many publications involving smallRNA use standard aligners (e.g. bowtie, BWA, etc..). In general, they align their reads to the genome and then subset their alignments based on metrics such as size and whether or not they match the locations of previously annotated RNA features. The major obstacle ...
written 8 months ago by benformatics150
3
votes
2
answers
513
views
2
answers
Check my work? Hyper geometric test on enrichment between iCLIP sites across experiments
... Hi I was hoping for some more feedback on this test I am trying to perform to check for enrichment at specific sites throughout the genome and my statistics background isn't great. _____ Here is the current setup: Given: - genomic coordinates of iCLIP binding sites (single nucleotide position corre ...
enrichment hypergeometric written 12 months ago by benformatics150 • updated 12 months ago by michael.ante1.8k
1
vote
3
answers
4.0k
views
3
answers
Answer: A: How To Get The Gene Name Interactively In A R Plot
... This is also super easy if you are already good with ggplot as [plotly][1] has added a library/function to do just this. All you need to do is make sure that your 'text' aesthetic is the text of interest. library(ggplot2) library(plotly) ## given test.df your data.frame with inter ...
written 13 months ago by benformatics150
1
vote
1
answer
560
views
1
answers
Answer: A: Duplicate miRNAs with different read counts in miRDeep2
... I believe this was answered previously: https://www.biostars.org/p/90159/ ...
written 13 months ago by benformatics150
0
votes
5
answers
925
views
5
answers
Answer: A: How Many Reads In A Bam File Are Aligned To a Specific Region
... An R code version using GenomicRanges. Not super fast but will do the job and let you keep working with the information in R. library(GenomicAlignments) library(rtracklayer) ## import bam file bamfile <- readGAlignments('sample.bam') ## assuming single end otherwise use read ...
written 13 months ago by benformatics150
1
vote
1
answer
869
views
1
answers
Answer: A: How to merge (or perform left joins) on GRange objects?
... I think the simplest way would be to do as you suggested and use the peakNames string. Your main problem overall however is the following: ![example case][1] You need to have a plan to address the above. And that is only a simple case you could have multiple chains like the above. Are you going ...
written 13 months ago by benformatics150
1
vote
1
answer
608
views
1
answers
Answer: A: Calculation of read counts in RNA-seq
... Given that they did this in R using `GenomicRanges`/`GenomicFeatures`, my assumption would be that the last section refers to the exact ranges that were used in the `summarizeOverlaps()` call. In this case they probably used their alignments (BAM files) and overlapped them with a `GRangesList` that ...
written 13 months ago by benformatics150

Latest awards to benformatics

Scholar 11 weeks ago, created an answer that has been accepted. For A: get mRNAs by TxDb
Appreciated 12 weeks ago, created a post with more than 5 votes. For A: ATAC-seq peak calling with MACS
Good Answer 12 weeks ago, created an answer that was upvoted at least 5 times. For A: ATAC-seq peak calling with MACS
Teacher 14 months ago, created an answer with at least 3 up-votes. For A: ATAC-seq peak calling with MACS

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 895 users visited in the last hour