Problem while modelling

[AngieD9110]

Hi,
I am trying to test modelling with following command:
tadbit model --workdir model_Ctrl_test --input_matrix test_Ctrl_matrix.f.bed --noX --optimize --beg 0 --end 20000000 --reso 50000 --maxdist 2000:4000:100 --nmodels 20 --nkeep 20 --exp_type Hi-C --cpu 32 --force
However, an error occurred with the minimum maxdist value:

No matter how I adjust the minimum maxdist , the error keeps showing with only the minimum maxdist:(here I changed to 1000:4000:100)

Could you tell me how this happen? Many thanks!

[david-castillo]

Hi,

Are you using the last version of TADbit (just type tadbit version)?
I remember a change back in May to keep compatibility with IMP where get_particle() was deprecated and that was changed to get().

Regards

David

[AngieD9110]

Thanks so much for your response. It seems conda can only update to version 1.0. I updated the package using GitHub and it works. However, as I try to do downstream analysis(e.g. correlation real/models) with the following command:
tadbit model --workdir model_Ctrl_test --input_matrix test_Ctrl_matrix.f.bed --noX --optimize --beg 0 --end 20000000 --reso 50000 --maxdist 1000:1100:100 --nmodels 20 --nkeep 20 --exp_type Hi-C --force --model --analyze --analyze_list 2
another problem occured:

[david-castillo]

Hi,

Yes, sorry, that's a bit unclear in the documentation but you should run three different commands: one for the optimization where you infer the optimal parameters, one for the creation of the final ensemble of models and one to create the analysis plots.
You can do several optimizations and tadbit will just add the combinations of parameters that haven't been tried in the previous runs. The idea is to maximize the correlation by trying different combinations of maxdist, lowfreq, upfreq and dcutoff. When you find the best correlation issue the tadbit model with --model.

tadbit model --optimize --workdir model_Ctrl_test --input_matrix test_Ctrl_matrix.f.bed --noX --beg 0 --end 20000000 --reso 50000 --maxdist 1000:1100:100 --nmodels 20 --nkeep 20 --exp_type Hi-C

tadbit model --model --workdir model_Ctrl_test --input_matrix test_Ctrl_matrix.f.bed --noX --beg 0 --end 20000000 --reso 50000 --nmodels 200 --nkeep 200 --exp_type Hi-C

tadbit model --analyze --workdir model_Ctrl_test

Regards

David