EmoryUniversityTheoreticalBiophysics / SirIsaac

Automated dynamical systems inference
MIT License
38 stars 11 forks source link

overflow error #26

Open gengyuncong opened 4 years ago

gengyuncong commented 4 years ago

Hi, I tested it with my own dataset and I got some errors. I used the same parameter as in the simpleExample.ipynb:

complexityStepsize = 2 # increase complexity with steps of size 2
complexityMax = 25 # don't try models with complexity > 25
complexityList = range(0,complexityMax,complexityStepsize) 
totalSteps = 1e3
keepSteps = 10
seeds = (1,1) # use a fixed random seed
ensTemperature = 100.
ensGen = fittingProblem.EnsembleGenerator( totalSteps, keepSteps,
    temperature=ensTemperature, seeds=seeds )
avegtol = 1e-2
maxiter = 100
priorSigma = 3.
numprocs = 3
N = len(indepParamsList)

Here is how my data looks like:

indepParamsList[:10 ]
[[1.0], [1.0], [1.0], [1.0], [1.0], [1.0], [1.0], [1.0], [1.0], [2.0]]

sirIsaacData[:3]
[{'x1': {3.0: (1.2346039866581013, 0.19381314496763835)},
  'x2': {3.0: (5.223171191926012, 0.5260055984594344)},
  'x3': {3.0: (7.072974151579377, 0.704560417875636)}},
 {'x1': {3.5: (1.9176442151604471, 0.7198853888751312)},
  'x2': {3.5: (3.871928103727898, 0.41154764318844306)},
  'x3': {3.5: (6.751740435788017, 0.842283495315975)}},
 {'x1': {4.5: (1.8585518592890489, 0.9657200975255507)},
  'x2': {4.5: (4.711265598482662, 0.5967139128320695)},
  'x3': {4.5: (9.877799185687941, 1.1240765611335537)}}]

And here is my error message:

error.txt