Skip to content

Simulations guide ChIPsim

Jump to bottom
Vaso Triantafyllidou edited this page Feb 8, 2021 · 3 revisions

ChIPsim Guide

  • download the ChIPsim code

  • All the chromosome files must be in the same directory, FASTA format and have the chr?.fa extension, in order to be accepted. The chr? part will be used to create the chr?.fastq and control_chr?.fastq files in the same directory.

  • It is required to specify the path to the directory with the chromosomes, and a number for the seed while it is optional to specify a number for the cores that will be used (default 1).

.\chipsim.py --dir /path-to-your-dir/ -s INT -n INT

usage: ChIPsim [-h] -d DIR -s S [-n N]

Simulated ChIPseq data for treatment and control using R package "ChIPsim"

optional arguments:
  -h, --help         show this help message and exit
  -n N               number of cores

required arguments:
  -d DIR, --dir DIR  directory with chromosome .fa files
  -s S               seed number for the simulation
  • The python code does not perform any additional check if the files provided are FASTA format, if there are nucleotides other than ATGCN etc, all these are handled by the ChIPsim package. Also we should note that although the user selects only one seed numbed, in reality (#chromosomes * 5 scripts) = # seed will be generated, and this is the reason that a dictionary assigning a number to each chromosome file is created before calling the ChIPsim package.

Clone this wiki locally