аюсь переписать некоторый код из R в C, поэтому я пытаюсь подогнать линейную регрессию, используяgsl_fit_linear() функция.

В R я бы использовал функцию lm (), которая возвращает p-значение для подгонки, используя этот код:

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

Я понятия не имею, как перейти от вывода C к p-значению, мой код выглядит примерно так:

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

Кажется, это правильно рассчитывает наклон и перехват, но я не знаю, как получить значение р. Я новичок в статистике и C!