Obtención del valor de p para la regresión lineal en la función gsl_fit_linear () de la biblioteca GSL

Estoy tratando de reportar algo de código desde R en C, así que estoy tratando de ajustar una regresión lineal usando la función gsl_fit_linear() .

En R usaría la función lm (), que devuelve un valor de p para el ajuste usando este código:

 lmAvgs<- lm( c(1.23, 11.432, 14.653, 21.6534) ~ c(1970, 1980, 1990, 2000) ) summary(lmAvgs) 

No tengo idea de cómo pasar de la salida de C a un valor de p, mi código se ve algo así hasta ahora:

 int main(void) { int i, n = 4; double x[4] = { 1970, 1980, 1990, 2000 }; double y[4] = {1.23, 11.432, 14.653, 21.6534}; double c0, c1, cov00, cov01, cov11, sumsq; gsl_fit_linear (x, 1, y, 1, n, &c0, &c1, &cov00, &cov01, &cov11, &sumsq); } 

Esto parece calcular correctamente la pendiente y la intersección, pero no sé cómo obtener un valor de p. Soy novato en stats y C!

Todo está en: http://en.wikipedia.org/wiki/Ordinary_least_squares . pero aquí hay un fragmento de código que muestra una salida similar al resumen (lmAvgs) en R. Para ejecutar esto, necesita la Biblioteca GSL :

 int n = 4; double x[4] = { 1970, 1980, 1990, 2000}; double y[4] = {1.23, 11.432, 14.653, 21.6534}; double c0, c1, cov00, cov01, cov11, sumsq; gsl_fit_linear (x, 1, y, 1, n, &c0, &c1, &cov00, &cov01, &cov11, &sumsq); cout<<"Coefficients\tEstimate\tStd. Error\tt value\tPr(>|t|)"< 

Lo que da :

 Coefficients Estimate Std. Error t value Pr(>|t|) Intercept -1267.91 181.409 -6.98922 0.0198633 x 0.644912 0.0913886 7.05681 0.0194956 Multiple R-squared: 0.961389, Adjusted R-squared: 0.942083 F-statistic: 49.7986 on 1 and 2 DF, p-value: 0.0194956 

R da:

 Coefficients: Estimate Std. Error t value Pr(>|t|) (Intercept) -1.268e+03 1.814e+02 -6.989 0.0199 * c(1970, 1980, 1990, 2000) 6.449e-01 9.139e-02 7.057 0.0195 * --- Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1 Residual standard error: 2.044 on 2 degrees of freedom Multiple R-squared: 0.9614, Adjusted R-squared: 0.9421 F-statistic: 49.8 on 1 and 2 DF, p-value: 0.01950