The present work deals with a Gompertz-type diffusion process, which includes in the drift term a time-dependent function C(t) representing the effect of a therapy able to modify the dynamics of the underlying process. However, in experimental studies is not immediate to deduce the functional form of C(t) from a treatment protocol. So a statistical approach is proposed in order to estimate this function when a control group and one or more treated groups are observed. In order to validate the proposed strategy a simulation study for several interesting functional forms of C(t) has been carried out. Finally, an application to infer the net effect of cisplatin and doxorubicin+cyclophosphamide in actual murine models is presented