It may be that we should make sure the fitting is robust to there being NaN values in the PRF. At the end of the day, if the PRF has NaN values, we should just ignore them with the fitting and fit only where there are valid model and data values?