The design matrix looks correct. 

Can you set P.wb=0 in your gPPI script. This should fix the issues. If it is already set to 0, then please let me know and I'll try to identify another solution.

This will revert the code back to slice-by-slice computations, rather than reading in the whole volume.

Best,
Donald