Comparing alternative estimation methods of variance-covariance matrices in the class of Fractionally Integrated GARCH models