Extract the transition matrix of jump problem


In DiffEqJump, I would like to extract the transition matrix assiociated to the jump process. Can you tell me if this is implemented please?

For example, when defining

sir_model = @reaction_network SIR begin
           c1, s + i --> 2i
           c2, i --> r
       end c1 c2

I would like to be able to get this matrix for the jump process sir_model.

Thank you for your help,