Random remark, are you sure you want to model i/r as normal, and allow it to be negative? Maybe you want a half-normal.