Hi！
Thanks for such a great tool and it's very helpful for my project. I followed the tutorial and it runs well from step1 to step5 of part1. 725 candidate regions are found. But for the part2, I want to found the difference between repressed regions(states 1,2,7,8) and active TSS regions(states 4,5,6) using the following commond.
sciddo.py --debug dump --data-file testdata.h5 --support-file test_run_penem.h5 --output test_run_penem_split-DCD.enh-to-pcr.tsv --data-type dynamics --split-segments --from-states 1 2 7 8 --to-states 4 5 6 --add-inverse --limit-bed-output
It gives the error:
Here is the design matrix with no replicates.
It runs well with the SCIDDO testdata. Is this caused by the data features or other?
Another question is can the SCIDDO identify the regions which are in the same states for the two groups but with significant different signal strength? It seems that the 725 candidate regions are the switching between the states.
Many thanks!
Nan
Hi！
Thanks for such a great tool and it's very helpful for my project. I followed the tutorial and it runs well from step1 to step5 of part1. 725 candidate regions are found. But for the part2, I want to found the difference between repressed regions(states 1,2,7,8) and active TSS regions(states 4,5,6) using the following commond.
sciddo.py --debug dump --data-file testdata.h5 --support-file test_run_penem.h5 --output test_run_penem_split-DCD.enh-to-pcr.tsv --data-type dynamics --split-segments --from-states 1 2 7 8 --to-states 4 5 6 --add-inverse --limit-bed-output
It gives the error:
Here is the design matrix with no replicates.
It runs well with the SCIDDO testdata. Is this caused by the data features or other?
Another question is can the SCIDDO identify the regions which are in the same states for the two groups but with significant different signal strength? It seems that the 725 candidate regions are the switching between the states.
Many thanks!
Nan