The aim of this paper is to give some results on the exact density of the I-divergence in the exponential family with gamma distributed observations. It is shown in particular that the I-divergence can be decomposed as a sum of two independent variables with known distributions. Since the considered I-divergence is related to the likelihood ratio statistics, we apply the method to compute the exact distribution of the likelihood ratio tests and discuss the optimality of such exact tests. One of these tests is the exact LR test of the model which is asymptotically optimal in the Bahadur sense. Numerical examples are provided to illustrate the methods discussed.
