User: benformatics

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

Posts by benformatics

<prev • 21 results • page 1 of 3 • next >
1
vote
1
answer
112
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 28 days ago by benformatics150
0
votes
2
answers
1.5k
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 5 weeks 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 6 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 6 months ago by benformatics150
3
votes
2
answers
455
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 10 months ago by benformatics150 • updated 10 months ago by michael.ante1.6k
1
vote
3
answers
3.8k
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 12 months ago by benformatics150
1
vote
1
answer
499
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 12 months ago by benformatics150
0
votes
5
answers
754
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 12 months ago by benformatics150
1
vote
1
answer
733
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 12 months ago by benformatics150
1
vote
1
answer
553
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 12 months ago by benformatics150

Latest awards to benformatics

Scholar 24 days ago, created an answer that has been accepted. For A: get mRNAs by TxDb
Appreciated 4 weeks ago, created a post with more than 5 votes. For A: ATAC-seq peak calling with MACS
Good Answer 4 weeks ago, created an answer that was upvoted at least 5 times. For A: ATAC-seq peak calling with MACS
Teacher 12 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: 925 users visited in the last hour