The paper presents a procedure based on the EM algorithm for the indirect estimation of the parameters of BiLinear GARCH (BL-GARCH) models. BL-GARCH generalize the class of GARCH models by considering interactions of past shocks and volatilities in the conditional variance equation. In this way...