<t>Late answer but anyway:you found the right answer, the parameter R(t) is used to recover exactly the initial yield curve.Indeed the sde dr(t) gives you r(t) in terms of R(t), and then you can calculate the zero coupon bond bythe affine formula P(t, T) = exp( A(t, T) - B(t, T) r(t) ),with A(t, T) ...