Feed the pm.fit results to NUTS sampler

You probably need to input the start dict once for each chain, so try
pm.sample(..., start=[start, start, start, start])