User: colin.kern

gravatar for colin.kern
colin.kern650
Reputation:
650
Status:
Trusted
Location:
United States
Last seen:
11 hours ago
Joined:
6 years, 10 months ago
Email:
c*********@gmail.com

about me

Posts by colin.kern

<prev • 100 results • page 1 of 10 • next >
0
votes
1
answer
73
views
1
answers
Comment: C: Why is the pValue peaks very high in EpigenomeRoadMap data.
... You can see the format specifications here: https://genome.ucsc.edu/FAQ/FAQformat.html#format13 The only difference is the narrowPeak format has one additional column that indicates the peak summit. ...
written 4 days ago by colin.kern650
0
votes
0
answers
69
views
0
answers
Comment: C: GC bias correction for Chip-SEQ
... The way current sequencing technology works introduces biases for GC-rich areas of the genome, so they tend to get more reads than other areas. This is hard to account for in ChIP-seq because the areas of interest (promoters, enhancers, etc.) often are GC-rich themselves, so the bias is expected bio ...
written 5 days ago by colin.kern650
0
votes
2
answers
76
views
2
answers
Answer: A: Chip-Seq macs2 peak calling
... The output files are explained in the README. Just go here: https://github.com/taoliu/MACS and look near the bottom of the page for the "Output Files" section. The narrowPeak file is an extension of the BED format, which is why it's used for bedtools. You are going to have be a lot more specific a ...
written 5 days ago by colin.kern650
2
votes
1
answer
73
views
1
answers
Answer: A: Why is the pValue peaks very high in EpigenomeRoadMap data.
... A raw p-value is in the range 0 to 1, so 19.2 is clearly not a p-value. The narrowPeak and broadPeak formats, which are standard formats for peak calls that are extended from the BED format, represent p- and q-values as -log10 transformed values. The Roadmap data structure is probably derived from t ...
written 5 days ago by colin.kern650
5
votes
2
answers
106
views
2
answers
Answer: A: Effect of gene size on WGCNA
... Like DEG analysis, WGCNA is looking at each gene independently of one another. All that matters is how the expression of the same gene compares across samples. That comparison will be proportionally the same whether you normalize by gene length or not, because all the values within each comparison a ...
written 6 days ago by colin.kern650
0
votes
3
answers
100
views
3
answers
Answer: A: Different number of peaks between biological replicates?
... Is this data from cell lines or frozen tissue? ChIP-seq can be a very hard assay to get consistent data from, especially with tissue, because there can be very large differences in the signal-to-noise ratio even when following the exact same protocol. Seeing that it is consistently the second replic ...
written 8 days ago by colin.kern650
2
votes
1
answer
88
views
1
answers
Answer: A: Detection of ncRNA in RNAseq transcript expression Data: lncRNA vs lincRNA
... LincRNA means "long intergenic non-coding RNA", whereas lncRNA is just "long non-coding RNA". Many people use them interchangeably. The only scenario I can think of where there might be a lncRNA that's not a lincRNA would be a lncRNA that overlaps a coding gene on the antisense strand, but I have a ...
written 8 days ago by colin.kern650
1
vote
0
answers
113
views
0
answers
Comment: C: MACS2 "total fragments in treatment" gives less number of fragments than expecte
... Samtools flagstat reports individual reads, not read pairs. The fragments in Macs2 output is referring to pairs of reads, so that's why you're getting half the number. The 370,935 difference may be Macs2 filtering read pair alignments that are considered "properly paired" in samtools but not by Macs ...
written 8 days ago by colin.kern650
1
vote
3
answers
129
views
3
answers
Answer: A: Why you shouldn't normalize count data using gene-length AND library size for co
... It depends on the analysis you're doing. Normalization for gene length isn't done in edgeR or DESeq2 because for DEG analysis you are only comparing counts of the same gene across samples, so normalizing by gene length won't change the proportions of the numbers of the same gene across samples. If y ...
written 8 days ago by colin.kern650
0
votes
1
answer
170
views
1
answers
Comment: C: Which one to use for ChipSeq replicate comparison: Spearman or Pearson
... Yes, skipZeros removes bins with values of zero in all samples, so it will remove bins in unmappable regions, but if you have just one read in one sample, then that bin is included, so it doesn't help to remove any of your background noise. TPM's benefit over FPKM/RPKM is well documented in RNA-seq ...
written 11 days ago by colin.kern650

Latest awards to colin.kern

Centurion 4 days ago, created 100 posts.
Appreciated 6 days ago, created a post with more than 5 votes. For A: Effect of gene size on WGCNA
Scholar 6 days ago, created an answer that has been accepted. For A: How do I interpret genotype likelihood in the VCF format?
Teacher 6 days ago, created an answer with at least 3 up-votes. For A: Quick Way to Annotate a Bed File
Scholar 7 days ago, created an answer that has been accepted. For A: How do I interpret genotype likelihood in the VCF format?
Popular Question 6 weeks ago, created a question with more than 1,000 views. For ChIP-seq Genome Coverage
Teacher 7 weeks ago, created an answer with at least 3 up-votes. For A: Quick Way to Annotate a Bed File
Scholar 4 months ago, created an answer that has been accepted. For A: How do I interpret genotype likelihood in the VCF format?
Scholar 4 months ago, created an answer that has been accepted. For A: How do I interpret genotype likelihood in the VCF format?
Teacher 4 months ago, created an answer with at least 3 up-votes. For A: Quick Way to Annotate a Bed File
Popular Question 9 months ago, created a question with more than 1,000 views. For Tophat align_summary.txt and samtools flagstat accepted_hits.bam disagree
Popular Question 9 months ago, created a question with more than 1,000 views. For ChIP-seq Genome Coverage
Teacher 9 months ago, created an answer with at least 3 up-votes. For A: Quick Way to Annotate a Bed File
Scholar 9 months ago, created an answer that has been accepted. For A: Quick Way to Annotate a Bed File
Popular Question 12 months ago, created a question with more than 1,000 views. For Tophat align_summary.txt and samtools flagstat accepted_hits.bam disagree
Popular Question 12 months ago, created a question with more than 1,000 views. For ChIP-seq Genome Coverage
Popular Question 12 months ago, created a question with more than 1,000 views. For Scaffold name mapping between NCBI and Ensembl
Popular Question 15 months ago, created a question with more than 1,000 views. For Tophat align_summary.txt and samtools flagstat accepted_hits.bam disagree
Scholar 16 months ago, created an answer that has been accepted. For A: Quick Way to Annotate a Bed File
Great Question 20 months ago, created a question with more than 5,000 views. For Can't convert paired end BAM to bed using bedtools
Popular Question 20 months ago, created a question with more than 1,000 views. For vcf-consensus Error: Fasta sequence does not match REF
Popular Question 20 months ago, created a question with more than 1,000 views. For Tophat align_summary.txt and samtools flagstat accepted_hits.bam disagree
Popular Question 20 months ago, created a question with more than 1,000 views. For ChIP-seq Genome Coverage
Popular Question 20 months ago, created a question with more than 1,000 views. For Get # of unique mapped reads from Tophat output

Help
Access

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