{
    Double_t data[3];
    TGraphErrors g("data.dat");
    TRobustEstimator rob(g.GetN(), 3);

    for (Int_t i = 0; i < g.GetN(); i++)
    {
        data[0] = g.GetX()[i];
        data[1] = g.GetY()[i];
        data[2] = g.GetEX()[i];
        rob.AddRow(data);
    }

    rob.Evaluate();
    rob.GetMean()->Print();
    rob.GetCovariance()->Print();
}