Utilizando o mesmo algoritmo da resposta do @Marcos Banikda resposta do @Marcos Banik, escrevi um pequeno algoritmo em C baseado no método dos mínimos quadrados:
void lms(double *x, double *y, int n, double *m, double *b)
{
int i;
double sumYX = 0.;
double sumX = 0.;
double sumY = 0.;
double sumX2 = 0.;
double sum2X = 0.;
for(i = 0; i < n; i++) {
sumYX += x[i] * y[i];
sumX += x[i];
sumY += y[i];
sumX2 += x[i] * x[i];
}
sum2X = sumX * sumX;
*m = (sumYX - (sumX * sumY) / (double)n) / (sumX2 - sum2X / (double)n);
*b = sumY / (double)n - *m * sumX / (double)n;
}
O algoritmo é exatamente o mesmo utilizado na resposta dele. O laço for
inicial é feito para calcular todos as somatórias e posteriormente seus resultados são utilizados para achar m
e b
. Acredito que isso seja facilmente traduzível para sua linguagem de preferência. :)
Com m
e b
você tem a equação da reta. Para calcular qual o valor de y, basta fazer uma função:
double trend(double m, double b, double x)
{
return m*x + b;
}
Para chamar as funções:
int main(void)
{
double x[] = {1., 2., 3., 4.};
double y[] = {10., 20., 30., 40.};
double m, b;
double nx = 5.;
double ny;
lms(x, y, 4, &m, &b);
ny = trend(m, b, nx);
printf("m: %lf \nb: %lf \nnx: %lf \nny: %lf\n", m, b, nx, ny);
return 0;
}