Currently doing this by converting to magnitude-phase form, but lmfit can also support complex data: https://lmfit.github.io/lmfit-py/faq.html#how-can-i-fit-complex-data Might be worth investigating.