/*UVOS*/

/* This file is part of MAClient copyright © 2021 Carl Philipp Klemm.
 *
 * MAClient is free software: you can redistribute it and/or modify
 * it under the terms of the GNU General Public License (GPL) version 
 * 3 as published by the Free Software Foundation.
 *
 * MAClient is distributed in the hope that it will be useful,
 * but WITHOUT ANY WARRANTY; without even the implied warranty of
 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
 * GNU General Public License for more details.
 *
 * You should have received a copy of the GNU General Public License
 * along with MAClient.  If not, see <http://www.gnu.org/licenses/>.
*/
#include "regessioncalculator.h"
#include <stdexcept>
#include <math.h>
#include <QDebug>

RegessionCalculator::RegessionCalculator(const std::vector<double>& xValues, const std::vector<double>& yValues)
{
    if(xValues.size() == yValues.size())
    {
        this->xValues = xValues;
        this->yValues = yValues;

        double sumY = 0;
        double sumX = 0;
        double sumXTimesY = 0;
        double sumSquaredX = 0;
        double sumSquaredY = 0;
        for (unsigned i = 0; i < xValues.size(); i++)
        {
            sumY += yValues[i];
            sumX += xValues[i];
            sumSquaredX += xValues[i]*xValues[i];
            sumSquaredY += yValues[i]*yValues[i];
            sumXTimesY += xValues[i]*yValues[i];
        }

        slope = (sumXTimesY - (sumX*sumY)/xValues.size()) / (sumSquaredX - (sumX*sumX)/xValues.size());
        offset = sumY/xValues.size() - slope*(sumX/xValues.size());

        double error = (xValues.size()*sumSquaredY-sumY*sumY-slope*slope*(xValues.size()*sumSquaredX-sumX*sumX)) / (xValues.size()*(xValues.size()-2));
        stdError = sqrt( (error*error*xValues.size() ) / ( xValues.size()*sumSquaredX-sumX*sumX));
    }
    else throw std::invalid_argument("xValues and yValues need to be the same size");
}