Closed jvanzand closed 1 year ago
Update: I changed to the basis
fitting_basis = 'per tc secosw sesinw logk'
and got a similar error later on in rvsearch.periodogram.Periodogram.per_bic():
Traceback (most recent call last):
File "/scr/vanzandt/miniconda3/envs/radvel_fix/bin/rvsearch", line 8, in <module>
sys.exit(main())
File "/scr/vanzandt/miniconda3/envs/radvel_fix/lib/python3.7/site-packages/rvsearch/cli.py", line 158, in main
args.func(args)
File "/scr/vanzandt/miniconda3/envs/radvel_fix/lib/python3.7/site-packages/rvsearch/driver.py", line 62, in run_search
searcher.run_search(outdir=args.output_dir)
File "/scr/vanzandt/miniconda3/envs/radvel_fix/lib/python3.7/site-packages/rvsearch/search.py", line 437, in run_search
perioder.per_bic()
File "/scr/vanzandt/miniconda3/envs/radvel_fix/lib/python3.7/site-packages/rvsearch/periodogram.py", line 203, in per_bic
self.post.params['k{}'.format(self.num_known_planets+1)].vary = True
KeyError: 'k6'
Looks like there's another assumption that 'k' is in the basis, whereas my basis contains 'logk' instead.
I think there is only one rvsearch-compatible basis
I ran a basic search using the following basis in my setup file:
This gave an error message:
In rvsearch.utils.initialize_default_pars(), the params object is converted from anybasis to fitting_basis, which is then assumed to contain the "per1" param:
Because my basis instead uses "logper," the search fails here.