User: bernatgel

gravatar for bernatgel
bernatgel230
Reputation:
230
Status:
Trusted
Location:
Barcelona, Spain
Twitter:
@bernatgel
Last seen:
2 hours ago
Joined:
3 years, 1 month ago
Email:
b********@gmail.com

about me

Posts by bernatgel

<prev • 23 results • page 1 of 3 • next >
1
vote
1
answer
92
views
1
answers
Comment: C: What does BAF mean?
... Actually, the definition of what is allele A and B is a [bit more complex][1] tha assigning calling allele A to the reference allele. What is allele A and allele B is sequence dependent and it's not related to the population allele frequency (which would be the case if the reference allele was allel ...
written 2 hours ago by bernatgel230
0
votes
2
answers
171
views
2
answers
Comment: C: Circos histogram from SNP.VCf file
... Circos only accepts data in its simple but strict tabular format and does not accept VCF as an input. You will need to either use another tool or preprocess the data until you have the exact data you want to plot in the format Circo needs. ...
written 5 days ago by bernatgel230
0
votes
1
answer
160
views
1
answers
Answer: A: SNP Array data with 0 copies for a segment
... If there's no normal contamination you should see a deep drop in LRR and a kind of a smear in BAF. BAF is computed as the part of the total signal that comes from the B allele. Since with an homozygous loss both alleles should have a signal 0, all signal you see will be background random noise. T ...
written 9 weeks ago by bernatgel230
1
vote
2
answers
420
views
2
answers
Comment: C: How to plot density between 2 bed files?
... I agree with Alex Reynolds that the small differences might come from reads counted multiple times. In the example code I've provided the windows are not overlapping but there's no check preventing a read to be counted twice in a similar way to bedops. In any case, if you need to ensure that every r ...
written 10 weeks ago by bernatgel230
0
votes
2
answers
420
views
2
answers
Comment: C: How to plot density between 2 bed files?
... Sure, you can specify a different genome when creating the karyoplot with `plotKaryotype(genome="hg38")` or `plotKaryotype(genome="mm10")`, for example. It has a small cache with a few frequently used organisms ans it the cache misses, it will use the available BSGenomes installed and download the c ...
written 11 weeks ago by bernatgel230
0
votes
2
answers
420
views
2
answers
Comment: C: How to plot density between 2 bed files?
... As far as I can see in GSM there are 4117 regions between positions 58819693 and 59363443 and in STEM there are 6958 regions between 58819362 and 59363565, so these peaks are real. I think the difference with your plots is due to the differences between GRCh37 and GRCh38. ChrY is shorter in GRCh38 ...
written 11 weeks ago by bernatgel230
3
votes
2
answers
420
views
2
answers
Answer: A: How to plot density between 2 bed files?
... If you prefer a pure R solution, you can user karyoploteR. The package is still in Bioconductor devel, so you will need R-devel to use it until mid April. You can compute the density and create the plot in a few lines: library(karyoploteR) window.size <- 1e5 #read the data with re ...
written 11 weeks ago by bernatgel230
0
votes
2
answers
177
views
2
answers
Comment: C: bedtools shuffle vs. random
... Yes, I forgot to add this to my answer. Strand specific randomization is in the pipeline but not ready yet. It is possible to do it in a strand specific way right now by defining a custom randomization function that internally randomizes separately according to strand. If you think you strand speci ...
written 11 weeks ago by bernatgel230
0
votes
2
answers
177
views
2
answers
Answer: A: bedtools shuffle vs. random
... If you can use R, you could user [regioneR][1] to test if there is a significant overlap between your dataset and the one from the publication. It will perform the whole process explained by @fridaymeetssunday: randomization a number of times (1000), computing the overlaps, the mean and standard dev ...
written 11 weeks ago by bernatgel230
3
votes
1
answer
606
views
1
answers
Answer: A: Testing for over-representation of chip-peaks
... You can use the R/Bioconductor package [regioneR][1] for this. It implements a statistical test for the association of genomic regions (such as chip peaks and exons) based on random permutations. In this case I think the best aproach would be to "flip" the question and ask whether alternatively spl ...
written 12 weeks ago by bernatgel230

Latest awards to bernatgel

Teacher 11 weeks ago, created an answer with at least 3 up-votes. For A: Сalculating fold-enrichment of ChIP-seq peaks intersecting with promoters (vs. g
Scholar 12 weeks ago, created an answer that has been accepted. For A: Permutation test with gene size and GC contents
Teacher 12 weeks ago, created an answer with at least 3 up-votes. For A: Сalculating fold-enrichment of ChIP-seq peaks intersecting with promoters (vs. g
Teacher 3 months ago, created an answer with at least 3 up-votes. For A: Сalculating fold-enrichment of ChIP-seq peaks intersecting with promoters (vs. g
Supporter 11 months ago, voted at least 25 times.
Good Answer 22 months ago, created an answer that was upvoted at least 5 times. For A: Permutation test with gene size and GC contents
Scholar 22 months ago, created an answer that has been accepted. For A: Permutation test with gene size and GC contents
Appreciated 22 months ago, created a post with more than 5 votes. For A: Permutation test with gene size and GC contents
Teacher 22 months ago, created an answer with at least 3 up-votes. For A: Сalculating fold-enrichment of ChIP-seq peaks intersecting with promoters (vs. g
Teacher 24 months ago, created an answer with at least 3 up-votes. For A: Сalculating fold-enrichment of ChIP-seq peaks intersecting with promoters (vs. g
Teacher 24 months ago, created an answer with at least 3 up-votes. For A: Association between bed files - statistical significance

Help
Access

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