Estou tentando repordar algum código de R em C, então estou tentando ajustar uma regressão linear usando ogsl_fit_linear() funçã

In R eu usaria a função lm (), que retorna um valor p para o ajuste usando este código:

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

Não tenho idéia de como ir da saída C para um valor p, meu código se parece com isso até agora:

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);
}

Isto parece calcular corretamente a inclinação e a interceptação, mas não sei como obter um valor-p. Sou iniciante em estatísticas e C!