September 15th, 2015, 12:39 am
QuoteOriginally posted by: fmfreshmanLet [$]R_{t}[$] be the CIR process at time [$]t[$]. I want to calculate the expectation $$E(R_{t}|R_{T})$$ conditional on [$]R_{T}[$] and [$]R_{0}[$] , where [$]t<T[$]. Then, is there a good way to calculate it?As pointed out, your [$]R_t[$] is a diffusion bridge process. It has a known time-inhomogeneous SDE, say eqn (68) here.The bridge SDE drift is known in terms of the standard (unbridged) SDE coefs and the standard (CIR) transition density, [$]p(t, R_t|R_0)[$], which is known. So,(*) [$]dR_t = \gamma(t,R_t; T, R_T) dt + \sigma \sqrt{R_t} dW_t [$],where [$]\gamma(\cdot)[$] is the coefficient of dt at (68).Integrating,(**) [$]E[R_t] - R_0 = E[\int_0^t \gamma(s,R_s;T, R_T) \, ds | R_0] = \int_0^{\infty} \int_0^t \gamma(s,y;T, R_T) \, p(s,y|R_0) \, ds \, dy[$].Note the last term of the integrand is, again, the known CIR transition density.So, unless I have made a mistake, that's it. [The expectation in the l.h.s. of (**) is what you seek].In Mathematica, for example, just give that double integral to NIntegrate and you are done.You could check my formula by doing a Monte Carlo with the bridged SDE process (*).=================================================================================Correction.Thinking about it some more, the first equality in (**) seems correct, but not the second.In the second, the last term of the integral should really be the transition density for the bridged process.Since we don't know that one, that seems to leave either a Monte Carlo or numerical PDE solution for [$]E[R_t][$].Alternatively, what I posted should be approximately correct if [$]t \ll T[$], since in that case the motion shouldbe practically unrestricted. Finally, googling "simulating diffusion bridges" will turn up some things.
Last edited by
Alan on September 14th, 2015, 10:00 pm, edited 1 time in total.