Hi Bijita,
The quantum regression theorem can only be used to calculate correlators with a specific type of time ordering (see Gardiner and Zoller quantum noise), which I believe is due to the inability to evolve the master equation in reverse. Unlike the Schrodinger equation, the master equation is not reversible due to a loss of information to the environment. However, in the steady-state, you can still calculate the negative time correlator you want by switching the operators between the two times and calculating a positive correlation.
I suppose it would make sense to change the function so it just does this when called, but one needs to be more careful when it's not the steady-state correlation.
Regards,
Kevin