User: Gary

gravatar for Gary
Gary10
Reputation:
10
Status:
New User
Location:
Last seen:
2 weeks, 1 day ago
Joined:
3 years, 10 months ago
Email:
y**********@gmail.com

Posts by Gary

<prev • 64 results • page 1 of 7 • next >
0
votes
1
answer
84
views
1
answers
Comment: C: How to produce the fraction of reads in peaks (FRiP) using DiffBind
... Hi Rory, Many thanks. Here you are. [http://68.181.92.180/~Gary/temporal/20190325_DiffBind_ATACseq_FRiP/dba_atac_peaks.RData][1] [http://68.181.92.180/~Gary/temporal/20190325_DiffBind_ATACseq_FRiP/dba_atac_counts.RData][2] By the way, I find that DiffBind can estimate FRiP values based on bam fi ...
written 25 days ago by Gary10
0
votes
1
answer
84
views
1
answers
Comment: C: How to produce the fraction of reads in peaks (FRiP) using DiffBind
... Hi Rory, Thank you very much. Below are their bam, bai, narrowPeak, and CSV files. Is it OK? Could you teach me the command lines how to save the DBA object (atac) before and after the call to dba, count()? After that, I can provide the files you need. Thank you very much. [http://68.181.92.180/~G ...
written 26 days ago by Gary10
0
votes
0
answers
42
views
0
answers
Enough cell number or not based on library complexity evaluation by ATACseqQC?
... Hi, I use ATACseqQC to estimate the library complexity of one ATAC-Seq data. When I don't trim any reads, it seems that our library complexity (cell number) is not enough because of a saturation plateau (fig. 1). However, it appears that the library complexity (cell number) is enough because of an ...
0
votes
1
answer
84
views
1
answer
How to produce the fraction of reads in peaks (FRiP) using DiffBind
... Hi, We have three mouse ATAC-Seq data. However, there is an error message when I would like to estimate their FRiP values using DiffBind below (i.e., Error in if (colnames(res)[i] == "Peak caller") { : argument is of length zero). Could you help me? Many thanks. Best, Gary > sessionInfo() ...
diffbind dba.count atac-seq frip dba.show written 27 days ago by Gary10 • updated 27 days ago by Rory Stark2.8k
0
votes
1
answer
398
views
1
answers
Comment: C: How about our ATAC-Seq data quality analyzed by ATACseqQC
... Dear Haibo, Now I can run your command helped by my local friend (the detail below). However, we still don't understand the sub-command: awk 'BEGIN{FS=OFS="\t"} function abs(v) {return v < 0 ? -v : v}; /^@/. Does it mean that it removes duplicate reads in a BAM file? Many thanks for your kindly ...
written 6 weeks ago by Gary10
0
votes
1
answer
398
views
1
answers
Comment: C: How about our ATAC-Seq data quality analyzed by ATACseqQC
... Dear Haibo, I really appreciate your help in resolving our ATAC-Seq problem. The information you provided is very helpful. I totally agree with you that I have to learn more about the Samtools, awk, and Unix commands. Hopefully, I can use your command correctly in the near future. Thank you once ag ...
written 7 weeks ago by Gary10
0
votes
1
answer
398
views
1
answers
Comment: C: How about our ATAC-Seq data quality analyzed by ATACseqQC
... Dear Julie and Haibo, Thank you for your very useful help. I am not sure whether you saw my responses. I integrate my additional questions below. Could you help me again? Many thanks. After removing mitochondrial reads, duplicated reads, and other alignments of no interest (the partial bamQC resul ...
written 7 weeks ago by Gary10
0
votes
1
answer
398
views
1
answers
Comment: C: How about our ATAC-Seq data quality analyzed by ATACseqQC
... Dear Haibo, After removing mitochondrial reads, duplicated reads, and other alignments of no interest (the partial bamQC results below), I reproduce the fragment size distribution and library complexity plots below. The library complexity profile is very different. Based on them, could you help to ...
written 7 weeks ago by Gary10
0
votes
1
answer
398
views
1
answers
Comment: C: How about our ATAC-Seq data quality analyzed by ATACseqQC
... Dear Haibo, Thank you for your very useful information. I have not removed duplicated reads in the fragment size distribution and library complexity figures. We didn't make a size selection but used a Qiagen MinElute PCR Purification Kit for the purification as suggested by Buenrostro, et al. (201 ...
written 7 weeks ago by Gary10
0
votes
1
answer
398
views
1
answers
Comment: C: How about our ATAC-Seq data quality analyzed by ATACseqQC
... Dear Haibo, Many thanks for your suggestion and information. I have tried blast on the web, but it doesn't work (the figures below). Maybe the file size is too large, and I have to try a standalone blast. But it needs more time for me to learn how to install and run a standalone blast. ![enter ima ...
written 7 weeks ago by Gary10

Latest awards to Gary

Popular Question 9 months ago, created a question with more than 1,000 views. For Could not find function "makeTranscriptDbFromGFF"
Popular Question 9 months ago, created a question with more than 1,000 views. For DiffBind installation doesn't work
Supporter 19 months ago, voted at least 25 times.
Popular Question 2.3 years ago, created a question with more than 1,000 views. For How to run the rpkm function of edgeR
Popular Question 2.3 years ago, created a question with more than 1,000 views. For DiffBind installation doesn't work

Help
Access

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