从GSL库获取C gsl_fit_linear()函数中的线性回归的p值
我正试图从C中重新定义R中的一些代码,所以我试图使用
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);
}
这似乎正确计算斜率和截距,但我不知道如何获得p值。我是统计数据和C的新手!
没有找到相关结果
已邀请:
1 个回复
氮顺
这使 :
R给出: