SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
mouse ENCODE transcription factor ChIP-seq data gene_x Bioinformatics 1 02-04-2013 09:22 PM
Correlation of CHiP-Seq with ENCODE transcription factor data Tjcooper Bioinformatics 0 11-23-2012 04:34 AM
ENCODE ChIP-seq data hwu Bioinformatics 5 10-04-2012 01:32 PM

Reply
 
Thread Tools
Old 02-08-2013, 06:35 PM   #1
apredeus
Senior Member
 
Location: Bioinformatics Institute, SPb

Join Date: Jul 2012
Posts: 151
Default ENCODE ChIP-seq Metrics (NSC,RSC)

Hello all,

I was wondering if there's a software package that calculates the metrics that are used by ENCODE project (described in http://genome.cshlp.org/content/22/9/1813 )

I'm interested in NSC, RSC and data for "cross-correlation - strand shift" plots.

Anybody else uses those for their experiments?

Thank you in advance.
apredeus is offline   Reply With Quote
Old 04-15-2013, 07:34 AM   #2
phansen
Junior Member
 
Location: berlin

Join Date: Feb 2011
Posts: 1
Default SPP package

Hi,

download and in install the spp R-package from here

http://compbio.med.harvard.edu/Supplements/ChIP-seq/

This package has a function 'get.binding.characteristics' which will do the job.

Peter
phansen is offline   Reply With Quote
Old 04-15-2013, 10:44 PM   #3
apredeus
Senior Member
 
Location: Bioinformatics Institute, SPb

Join Date: Jul 2012
Posts: 151
Default

Awesome, thank you for the advice!
apredeus is offline   Reply With Quote
Old 05-29-2013, 03:11 PM   #4
mthornton
Junior Member
 
Location: LA

Join Date: May 2013
Posts: 2
Default

Hello,

I would use this script

https://code.google.com/p/phantompeakqualtools/

It uses SPP as well. It's convenient.
mthornton is offline   Reply With Quote
Old 10-05-2016, 01:42 AM   #5
zgene
Junior Member
 
Location: Asia

Join Date: Sep 2012
Posts: 7
Default

Quote:
Originally Posted by mthornton View Post
Hello,

I would use this script

https://code.google.com/p/phantompeakqualtools/

It uses SPP as well. It's convenient.
Dear All,

I am a newbie in ChIP-seq analysis and now I am trying to assess these values using "phantompeakqualtools", but unfortunately I got an error and cant find any similar posts for that.
Would you please help me !

I used this command to determine strand cross-correlation peak:

Rscript /data/NGS/phantompeakqualtools/run_spp.R -c=ChIP-WT1-R2-BWA_sorted_unique.bam -savp -out=R2unique-phantom

and I got the following error:

awk: line 2: function and never defined
Error converting BAM to tagalign file: /tmp/RtmpbkYntW/ChIP-WT1-R2-BWA_sorted_unique.bam6f036a95b6e

I assumed that the error was due to the input file and tried the following command to first convert BAM to tagalign

samtools view -F 0x0204 -o - ChIP-WT1-R1_sorted_unique.bam | awk 'BEGIN{OFS="\t"}{if (and($2,16) > 0) {print $3,($4-1),($4-1+length($10)),"N","1000","-"} else {print $3,($4-1),($4-1+length($10)),"N","1000","+"} }' | gzip -c > ChIP-WT1-R1_sorted_unique.tagAlign.gz

and again I got the following error;
awk: line 2: function and never define

So, now I know that the problem is due to awk. But I can't understand what it is exactly!

I would be happy to hear your comments.

Thanks in advance.
zgene is offline   Reply With Quote
Reply

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off




All times are GMT -8. The time now is 10:13 PM.


Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2020, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO