4
我想从C中重新编制一些R的代码,所以我尝试使用gsl_fit_linear()
函数来拟合线性回归。从GSL库的C gsl_fit_linear()函数中获得线性回归的p值
在R I会用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值。我是stats和C的新手!