gymreklab / chips

Simulation tool for ChIP- and other -seq experiments
GNU General Public License v3.0
12 stars 3 forks source link

Running WCE on high read lengths #46

Closed mlamkin7 closed 4 years ago

mlamkin7 commented 4 years ago

Below was the command ran which encountered what seemed to be an infinite loop. However, it finished to completion when the readlength was 250.

/storage/mlamkin/projects/tulip/src/tulip simreads -t wce --model params/GM12878_CTCF_ENCFF124LLI_ENCFF256QBB-1.9.json -f /storage/resources/dbase/human/hg19/hg19.fa -o fastqs/readlen500_wce --scale-outliers -c 7 --numcopies 1000 --numreads 10000 --readlen 500 --region chr22:1-51304566

GM12878_CTCF_ENCFF124LLI_ENCFF256QBB-1.9.json

{
    "frag": {
        "k": 21,
        "theta": 9.5
    },
    "pcr_rate": 0.80,
    "pulldown": {
        "f": 0.0003432990051805973,
        "s": 0.3999766409397125
    }
}