7

Tengo una función que me calcula la regresión lineal, pero falta que calcule el p valor. ¿Alguien sabría de manera sería?

function linearRegression(y,x)
{
    var lr = {};
    var n = y.length;
    var sum_x = 0;
    var sum_y = 0;
    var sum_xy = 0;
    var sum_xx = 0;
    var sum_yy = 0;

    for (var i = 0; i < y.length; i++) 
    {
        sum_x += parseFloat(x[i]);
        sum_y += parseFloat(y[i]);
        sum_xy += (parseFloat(x[i])*parseFloat(y[i]));
        sum_xx += (parseFloat(x[i])*parseFloat(x[i]));
        sum_yy += (parseFloat(y[i])*parseFloat(y[i]));
    } 

    lr['slope'] = (n*sum_xy - sum_x*sum_y) / (n*sum_xx - sum_x*sum_x);
    lr['intercept'] = (sum_y - lr.slope*sum_x)/n;
    lr['r2'] = Math.pow((n*sum_xy - sum_x*sum_y)/Math.sqrt((n*sum_xx-sum_x*sum_x)*(n*sum_yy - sum_y*sum_y)),2).toFixed(2);

    return lr;
}
5
  • 1
    Hola Aaron, ¿es un problema matemático o de programación?
    – César
    Commented el 24 may. 2016 a las 12:12
  • Buenas tardes, pues lo que necesito es un código parecido que realice el calculo de p valor porque desde esta página realiza el cálculo de todos los valores estadísticos: home.ubalt.edu/ntsbarsh/business-stat/otherapplets/…
    – Aaron VC
    Commented el 25 may. 2016 a las 15:51
  • 5
    Creo que iría bien que intentases hacer tu el calculo del p valor y si te falla algo nosotros te ayudamos.. A la mayoría de nosotros la época de estudiante nos pilla bastante lejos y no recordamos ciertos conceptos matemáticos.Aquí estas pidiendo que hagamos el código de un concepto matemático no trivial. Commented el 24 jun. 2016 a las 7:37
  • Tal vez te sirva esto... Regresión lineal simple, y viendo el código fuente de esa página podrás ver cómo funciona. Quizás hasta puedas simplificar el código.
    – user13558
    Commented el 13 sept. 2016 a las 19:38
  • 1
    Estaría muy bien si además del código explicaras un poco lo que significa el p valor y la regresión lineal. O al menos los datos de entrada y lo que supuestamente te tendría que devolver. Commented el 1 nov. 2016 a las 14:10

1 respuesta 1

4
+50

Aqui te dejo el codigo extraido de esta web, adaptado a tu función.

function linearRegression(y,x)
{
    var lr = {};
    var n = y.length;
    var sum_x = 0;
    var sum_y = 0;
    var sum_xy = 0;
    var sum_xx = 0;
    var sum_yy = 0;
    var cantidad = y.length;
    var varianzaX = 0.0;
    var varianzaY = 0.0;
    var varianza = 0.0
    var promedioX = 0.0;
    var promedioY = 0.0;

    for (var i = 0; i < n; i++) {
        sum_x += parseFloat(x[i]);
        sum_y += parseFloat(y[i]);
        sum_xy += (parseFloat(x[i])*parseFloat(y[i]));
        sum_xx += (parseFloat(x[i])*parseFloat(x[i]));
        sum_yy += (parseFloat(y[i])*parseFloat(y[i]));
    }

    lr['slope'] = (n*sum_xy - sum_x*sum_y) / (n*sum_xx - sum_x*sum_x);
    lr['intercept'] = (sum_y - lr.slope*sum_x)/n;
    lr['r2'] = Math.pow((n*sum_xy - sum_x*sum_y)/Math.sqrt((n*sum_xx-sum_x*sum_x)*(n*sum_yy - sum_y*sum_y)),2).toFixed(2);
    lr['p'] = 0;
    
    //
    promedioX = sum_x / cantidad;
    promedioY = sum_y / cantidad;
    
    // Calculo de varianzas
    for (i = 0; i < cantidad; i++) {
      varianzaX = varianzaX + Math.pow(x[i] - promedioX, 2);
      varianzaY = varianzaY + Math.pow(y[i] - promedioY, 2);
      varianza += ((x[i] - promedioX) * (y[i] - promedioY));
    }
    
    // P-value
    var SLP = varianza / varianzaX;
    var SSEF = (varianzaY - SLP * SLP * varianzaX);
    var SSRF = SLP * SLP * varianzaX;
    var MSEF = SSEF / (cantidad - 2);
    var FVAL1 = SSRF/MSEF;
    var Fp = FVAL1 + "";
    var fDesde = 1;
    var fHasta = cantidad - 2;
    var p = Fmt(FishF(FVAL1,fDesde,fHasta));
    if ((FVAL1 + "").indexOf("Infinity") == -1) {
      if ((p + "").indexOf("e") != -1) {
        lr['p'] = "Almost Zero";
      }
      else {  
        lr['p'] = p;
      }
    } 

    return lr;
}

function FishF(f, n1, n2) {
    var x = n2 / (n1 * f + n2)
    if ((n1 % 2) == 0) {
        return StatCom(1 - x, n2, n1 + n2 - 4, n2 - 2) * Math.pow(x, n2 / 2);
    }
    if ((n2 % 2) == 0) {
        return 1 - StatCom(x, n1, n1 + n2 - 4, n1 - 2) * Math.pow(1 - x, n1 / 2);
    }
    var th = Math.atan(Math.sqrt(n1 * f / n2));
    var a = th / (Math.PI / 2);
    var sth = Math.sin(th);
    var cth = Math.cos(th);
    if (n2 > 1) {
        a = a + sth * cth * StatCom(cth * cth, 2, n2 - 3, -1) / (Math.PI / 2);
    }
    if (n1 == 1) {
        return 1 - a;
    }
    var c = 4 * StatCom(sth * sth, n2 + 1, n1 + n2 - 4, n2 - 2) * sth * Math.pow(cth, n2) / Pi;
    if (n2 == 1) {
        return 1 - a + c / 2;
    }
    var k = 2;
    while (k <= (n2 - 1) / 2) {
        c = c * k / (k - .5);
        k = k + 1;
    }
    return 1 - a + c;
}

function StatCom(q, i, j, b) {
    var zz = 1;
    var z = zz;
    var k = i;
    while (k <= j) {
        zz = zz * q * k / (k - b);
        z = z + zz;
        k = k + 2;
    }
    return z
}

/* FMT: Redondea a 4 decimales */
function Fmt(x) {
    var v;
    if (x >= 0) {
      v = '' + (x + 0.00005);
    } else {
      v = '' + (x - 0.00005);
    }
    return v.substring(0, v.indexOf('.') + 5);
}

x = [1,1,0,1,0,1,0,1];
y = [1,1,0,1,1,1,1,1];
document.getElementById('resultado').innerHTML = JSON.stringify(linearRegression(y,x));
<div id="resultado"></div>

6
  • 3
    Excelente respuesta. Por cierto, el término estadístico es varianza (no variante o variación :-).
    – Mariano
    Commented el 3 nov. 2016 a las 5:48
  • 2
    Gracias Mariano, muy acertada tu observación! :)
    – Marcos
    Commented el 3 nov. 2016 a las 11:23
  • 4
    Para mi, la respuesta parece funcionar, pero esta no es Excelente como se dijo, debería incluir una buena explicación de las correcciones hechas y los errores cometidos. El asunto de fondo es aprender, no solamente desbloquearse de los problemas. Commented el 6 nov. 2016 a las 3:58
  • Tenés razón @rnd. En realidad quise decir "muy bueno el código". Es cierto que me llevó tiempo leerlo y entenderlo, cosa que se podría facilitar describiendo qué hace cada fragmento, por lo menos comentándolo en las partes más importantes... y quizás agregando también una imagen con una curva de distribución, mostrando distancias.
    – Mariano
    Commented el 7 nov. 2016 a las 10:21
  • 1
    Coincido en que mi respuesta solo brinda una solución al problema y carece de una explicación detallada. Siéntanse libres de tomar el código y documentarlo. Esta semana estoy demasiado cargado como para dedicarle el tiempo que se merece a esta respuesta. Saludos!
    – Marcos
    Commented el 7 nov. 2016 a las 12:29

Tu Respuesta

By clicking “Publica tu respuesta”, you agree to our terms of service and acknowledge you have read our privacy policy.

¿No es la respuesta que buscas? Examina otras preguntas con la etiqueta o formula tu propia pregunta.