example_ra_fitting.cpp

 
#include <Pds/Ra>
#include <Pds/Ml>
 
double Func(const Pds::Vector &x,const Pds::Vector &c)
{
    if(c.IsEmpty()) return 0;
    return c[0]+c[1]*x[0];
}
 
Pds::Vector DcFunc(const Pds::Vector &x,const Pds::Vector &c)
{
    if(c.IsEmpty()) return Pds::Vector();
    
    Pds::Vector Dc(c.Nel());
 
    Dc.SetRaw(0,1.0);
    Dc.SetRaw(1,x[0]);
    return Dc;
}
 
 
int main(void)
{
    Pds::IterationConf Conf;
    unsigned int L=1000;
    Pds::Matrix X;
    Pds::Vector Y;
    Pds::DataSet::Regression::LoadDataLine(L,X,Y);
    Pds::Vector Yres;
    
    Pds::Vector C="-0.3 0.1";
    
    C.Print("C:\n");
    C=Pds::Regression::FitModel(Conf,X,Y,Func,DcFunc,C);
    C.Print("C:\n");
 
    Yres=X.OperateRows(Func,C);
    
    Pds::Octave::Plot::CurveXYXY(X,Y,X,Yres,"testando.m","testando.png");
 
}