source: branches/rate/sample/estimation/alumina.mso @ 502

#using "types";

FlowSheet Alumina_Flow
        VARIABLES
        y as Real;
        x1 as Real;
        x2 as Real;

        PARAMETERS
        b0 as Real(Default=1);
        b1 as Real(Default=1);
        b2 as Real(Default=1);

        SPECIFY
        x1 = 0.647;
        x2 = 77.1;
        
        EQUATIONS
        y = b0 + b1*x1 + b2*x2;
        
        OPTIONS
        Dynamic = false;
        #NLASolver(AbsoluteAccuracy=1e-13);
end

Estimation Alumina_Est as Alumina_Flow

        ESTIMATE
        # PAR   START  LOWER  UPPER  UNIT
        b0     -1.00    -5  5;
        b1      2.00    -5  5;
        b2     -0.00    -5  5;

        EXPERIMENTS
        # FILE                WEIGTH    TYPE
# This sequence works correctly, but the simulation results go all to the first level of variables (fix me)
        "alumina_est.dat"               1               "fit";
        "alumina2_est.dat"              1               "test";

# This sequence does not show all data (fix me)
#       "alumina2_est.dat"              1               "fit";
#       "alumina_est.dat"      1                "test";

# This sequence does not compute statistics (fix me)
#       "alumina2_est.dat"              1               "test";
#       "alumina_est.dat"               1               "fit";

# This sequence has both problems above (fix me)
#       "alumina_est.dat"               1               "test";
#       "alumina2_est.dat"              1               "fit";

        OPTIONS
        Statistics(
                Fit=true,
                Parameter=true,
                Prediction=false
        );
        Significance = 0.95;
        BiLateral = true;
        NumJac = true;
        
        NLPSolver(
                AbsoluteAccuracy=1e-9,
                #File = "complex"
                File = "ipopt_emso"
        );
        Dynamic = false;
end