c语言gsl,从GSL库获取C gsl_fit_linear()函数中的线性回归的p值
一切都在:
http://en.wikipedia.org/wiki/Ordinary_least_squares.但是这里有一段代码,在R中显示类似于summary(lmAvgs)的输出.要运行它,你需要
GSL Library:
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<|t|)"<
double stdev0=sqrt(cov00);
double t0=c0/stdev0;
double pv0=t0<0?2*(1-gsl_cdf_tdist_P(-t0,n-2)):2*(1-gsl_cdf_tdist_P(t0,n-2));//This is the p-value of the constant term
cout<
double stdev1=sqrt(cov11);
double t1=c1/stdev1;
double pv1=t1<0?2*(1-gsl_cdf_tdist_P(-t1,n-2)):2*(1-gsl_cdf_tdist_P(t1,n-2));//This is the p-value of the linear term
cout<
double dl=n-2;//degrees of liberty
double ym=0.25*(y[0]+y[1]+y[2]+y[3]); //Average of vector y
double sct=pow(y[0]-ym,2)+pow(y[1]-ym,2)+pow(y[2]-ym,2)+pow(y[3]-ym,2); // sct = sum of total squares
double R2=1-sumsq/sct;
cout<
double F=R2*dl/(1-R2);
double p_value=1-gsl_cdf_fdist_P(F,1,dl);
cout<
这使 :
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给出:
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
c语言gsl,从GSL库获取C gsl_fit_linear()函数中的线性回归的p值相关推荐
- PHP获取二维数组中某一列的值集合
PHP还是比较常用的,于是我研究了一下PHP二维数组,下面通过本文给大家介绍PHP获取二维数组中某一列的值集合,对php数组二维数组的值相关知识感兴趣的朋友一起学习吧 PHP还是比较常用的,于是我研究 ...
- 如何获取Html或Jsp中select框内的值(JS原生与非原生)
JS原生 第一步 首先我们先给select定义一个Id号,我们现在给它定义为aaa <select id="aaa"> <option>选项一</op ...
- java获取插入数据库表中的自增列值的一个简单方法
PreparedStatement ps = conn.prepareStatement(sql,Statement.RETURN_GENERATED_KEYS); // 第一步 ps.execute ...
- python中os.system.获取输出信息_python中os.system()的返回值
[python中os.system()的返回值] 如果第三方程序返回的是布尔型返回值,os.system会将true转为1,false转为0进行返回. 问题: /bin/xxx.py是一个返回码为1的 ...
- JavaScript封装异步函数 —— 【异步编程】 —— 如何获取一个函数中异步操作的结果?
我们在写项目的时候有没有遇到这样的一个情况:一个函数里,嵌套着另一个异步函数,那如何在外面获取这个异步函数的结果呢?(return值) 就比如这样一个例子: function fn() ...
- C语言-----函数中的变量
目录 1.局部变量与全局变量 1.1 局部变量 1.2 全局变量 2.变量的存储类别 2.1 自动类型 (auto) 2.2 寄存器类型 (register) 2.3 静态类型 (static) 2. ...
- c语言 算法库 设计 gsl,GSL科学计算库、随机变量的Erlang分布与Weibull分布
更新:本文升级版"GSL科学计算库.随机变量的Erlang分布与Weibull分布"已经迁移至我的新博客http://gnailuy.com/.新文章对已有内容做了修改,并新增关于 ...
- gsl科学计算库文档,翻译了索引,凑合看看。
1.介绍 2.库的使用 本章描述如何编译使用GSL的程序,介绍GSL的一般用法. 2.1例子程序 2.2编译和链接 2.3共享库 2.4与ANSI C兼容性 2.5inline函数 2.6长双精度lo ...
- GSL科学计算库——计算高斯-勒让德积分
相关文章: Windows系统Qt5配置GSL科学计算库 Windows系统下Cygwin+Dev C ++ 配置GSL科学计算库 假设计算下列积分: ∫0πexcos(x)dx\int_0^\pi ...
最新文章
- mysql 删除 修改密码_Mysql数据库root密码忘记了,如何在不删除Mysql的情况下修改密码...
- 3.Spring Security 自定义用户认证
- 使用Azure portal Create Virtual Machine
- matlab计算函数区域面积,matlab 中函数regionprops 求解区域面积的原理
- protel99se 问题汇总(不定期更新)
- JavaScript 复习之 事件模型 和 Event对象
- (二)注册服务提供者
- 决策者根据什么曲线做出决策_如何在开放社区中做出品牌决策
- Android属性动画赏析,Android源码分析—属性动画的工作原理
- 在Filfter中使用注解了spring 的bean,报javax.naming.NamingException: Cannot create resource instance...
- 项目如何开始:怎样和客户一起搞定需求
- JavaWeb 利用jsp 实现分页查询
- pc电脑版java模拟器下载安装_PC电脑版java模拟器
- 在Visual Studio中使用Visual Basic教程:创建一个控制台应用
- PCWORLD评出的十大杀毒软件
- 《从零开始学架构:照着做,你也能成为架构师》李运华 读后感
- nvidia控制面板点了没反应win7_nvidia控制面板点击没反应 - 卡饭网
- 【Unity3D开发小游戏】《青蛙过河》Unity开发教程
- C++开发斗地主(QT)第三篇之动画发牌与位置计算
- 神操作!“AI程序员”走进偶像剧,网友:我服!