I would appreciate any help or comment if anybody experienced same problem.
I am trying to analyze ChIPseq data with PeakSeq v1.1 on MacOS.
I provided mapability file from http://archive.gersteinlab.org/proj/...Map/H.sapiens/ and I created chromosome list file by putting chromosome names on separate line (chr1, chr2, chr3,..)
Processing SAM file with -preprocess option seemed work fine and .bin files were created for each chromosome.
The problem comes in next step when I try to run it with -peak_select option.
It returns weird kind of error message:
Anybody has an idea what can cause this?
here is the config.dat file I am using:
I am trying to analyze ChIPseq data with PeakSeq v1.1 on MacOS.
I provided mapability file from http://archive.gersteinlab.org/proj/...Map/H.sapiens/ and I created chromosome list file by putting chromosome names on separate line (chr1, chr2, chr3,..)
Processing SAM file with -preprocess option seemed work fine and .bin files were created for each chromosome.
The problem comes in next step when I try to run it with -peak_select option.
It returns weird kind of error message:
Code:
PeakSeq -peak_select config.dat for reading.n Chromosome_list_HG.txt
here is the config.dat file I am using:
Code:
# Experiment id is used as a prefix to the output file name. Experiment_id IP_vs_Inp # Chromosome ID list file, is used to generate data file paths. chromosome_list_file Chromosome_list_HG.txt # Enrichment fragment length For tag extension, this is the value of average fragment length. Enrichment_fragment_length 200 # Target FDR in the simulations. target_FDR 0.05 # Number of simulations performed while estimating the putative peaks. N_Simulations 10 # Minimum distance between consecutive peaks Minimum_interpeak_distance 200 # Mappability file that includes the uniquely mappable number of nucleotides per window for each chromosome. Mappability_map_file Mapability_HG.txt # The directory that contains the preprocessed ChIP-Seq reads, can specify multiple directories to pool reads from multiple source (e.g. replicates) ChIP_Seq_reads_data_dirs PeakSeq_IP # The directory that contains the preprocessed Input (control) experiment reads. (Multiple directories allowed) Input_reads_data_dirs PeakSeq_Inp # Seed for pseudo-random number generator. This is necessary for simulated background option (specified below). #Simulation_seed 1234567 # Q-value threshold applied on the final set of peaks. max_Qvalue 0.05 # There are currently two models for simulating the background for threshold selection # Simulated background is the simulation based method that is explained in the PeakSeq paper. # Poisson background uses a simple Poisson background with mean estimated from the read statistics. This option is still experimental but it is much faster than the simulated background option. # Background_model Poisson Background_model Simulated