OK - the only correlation appears to be the IMU/magnetic yaw discrepancy, which is a clear function of yaw. That means, I think, that there is a component of the magnetic field of the aircraft that is not being subtracted out by calibration, perhaps because it is too large, or that the compass is faulty. You could try demagnetizing next.