popgenmethods / smcpp

SMC++ infers population history from whole-genome sequence data.
GNU General Public License v3.0
149 stars 32 forks source link

assertion error #245

Open ndussex opened 1 year ago

ndussex commented 1 year ago

Hi,

I suddenly get an error when trying to estimate the split between 2 populations and I didn't have this error before when using the following command:

smc++ split -o split/NovZ_NE NovZ/NovZ.final.json NE/NE.final.json results_divergence_NovZ_NE/.smcgz -c 100000

I have this error using either v1.15.2 or 1.15.4

Here's the error below. Has anyone seen this before?

Thanks a lot!

My error:

Traceback (most recent call last): File "/sw/bioinfo/SMC++/1.15.2/rackham/bin/smc++", line 8, in sys.exit(main()) File "/sw/bioinfo/SMC++/1.15.2/rackham/lib/python3.8/site-packages/smcpp/frontend/console.py", line 26, in main cmds[args.command].main(args) File "/sw/bioinfo/SMC++/1.15.2/rackham/lib/python3.8/site-packages/smcpp/commands/split.py", line 44, in main analysis = SplitAnalysis(args.data, args) File "/sw/bioinfo/SMC++/1.15.2/rackham/lib/python3.8/site-packages/smcpp/analysis/split.py", line 23, in init self._init_inference_manager(args.polarization_error, hs) File "/sw/bioinfo/SMC++/1.15.2/rackham/lib/python3.8/site-packages/smcpp/analysis/base.py", line 108, in _init_inference_manager im.model = self._model File "smcpp/_smcpp.pyx", line 198, in smcpp._smcpp.PyInferenceManager.model.set File "/sw/bioinfo/SMC++/1.15.2/rackham/lib/python3.8/site-packages/smcpp/observe.py", line 19, in f(self, *args, **kwargs) File "smcpp/_smcpp.pyx", line 328, in smcpp._smcpp.PyOnePopInferenceManager.update File "/sw/bioinfo/SMC++/1.15.2/rackham/lib/python3.8/site-packages/smcpp/model.py", line 297, in for_pop assert self.model1.N0 == self.model2.N0 AssertionError