文档章节

(C#)直线的最小二乘法线性回归

北风其凉
 北风其凉
发布于 2014/06/05 09:52
字数 450
阅读 400
收藏 0

1.Point结构

在编写C#窗体应用程序时,因为引用了System.Drawing命名空间,其中自带了Point结构,本文中的例子是一个控制台应用程序,因此自己制作了一个Point结构

/// <summary>
/// 二维笛卡尔坐标系坐标
/// </summary>
public struct Point
{
    public double X;
    public double Y;
    public Point(double x = 0, double y = 0)
    {
        X = x;
        Y = y;
    }
}

2.线性回归

/// <summary>
/// 对一组点通过最小二乘法进行线性回归
/// </summary>
/// <param name="parray"></param>
public static void LinearRegression(Point[] parray)
{
    //点数不能小于2
    if (parray.Length < 2)
    {
        Console.WriteLine("点的数量小于2,无法进行线性回归");
        return;
    }

    //求出横纵坐标的平均值
    double averagex = 0, averagey = 0;
    foreach (Point p in parray)
    {
        averagex += p.X;
        averagey += p.Y;
    }
    averagex /= parray.Length;
    averagey /= parray.Length;

    //经验回归系数的分子与分母
    double numerator = 0;
    double denominator = 0;

    foreach (Point p in parray)
    {
        numerator += (p.X - averagex) * (p.Y - averagey);
        denominator += (p.X - averagex) * (p.X - averagex);
    }

    //回归系数b(Regression Coefficient)
    double RCB = numerator / denominator;

    //回归系数a
    double RCA = averagey - RCB * averagex;

    Console.WriteLine("回归系数A: " + RCA.ToString("0.0000"));
    Console.WriteLine("回归系数B: " + RCB.ToString("0.0000"));
    Console.WriteLine(string.Format("方程为: y = {0} + {1} * x",
        RCA.ToString("0.0000"), RCB.ToString("0.0000")));

    //剩余平方和与回归平方和
    double residualSS = 0;   //(Residual Sum of Squares)
    double regressionSS = 0; //(Regression Sum of Squares)

    foreach (Point p in parray)
    {
        residualSS +=
            (p.Y - RCA - RCB * p.X) *
            (p.Y - RCA - RCB * p.X);

        regressionSS +=
            (RCA + RCB * p.X - averagey) *
            (RCA + RCB * p.X - averagey);
    }

    Console.WriteLine("剩余平方和: " + residualSS.ToString("0.0000"));
    Console.WriteLine("回归平方和: " + regressionSS.ToString("0.0000"));
}

3.Main函数调用

static void Main(string[] args)
{
    //设置一个包含9个点的数组
    Point[] array = new Point[9];
    array[0] = new Point(0, 66.7);
    array[1] = new Point(4, 71.0);
    array[2] = new Point(10, 76.3);
    array[3] = new Point(15, 80.6);
    array[4] = new Point(21, 85.7);
    array[5] = new Point(29, 92.9);
    array[6] = new Point(36, 99.4);
    array[7] = new Point(51, 113.6);
    array[8] = new Point(68, 125.1);
    LinearRegression(array);

    Console.Read();
}

4.运行结果

© 著作权归作者所有

共有 人打赏支持
北风其凉

北风其凉

粉丝 114
博文 498
码字总数 463468
作品 4
朝阳
程序员
Machine Learning 2 - 非线性回归算法分析

2017-08-02@erixhao 技术极客TechBooster AI 机器学习第二篇 - 非线形回归分析。我们上文深入本质了解了机器学习基础线性回归算法后,本文继续研究非线性回归。 非线性回归在机器学习中并非热...

erixhao
2017/08/07
0
0
机器学习中的回归理解

机器学习中回归(自变量x和应变量y)之间的关系由两种? 确定性关系(函数关系):如线性相关,能够找到函数或者曲线能够进行拟合。 非确定性关系(相关关系):如正相关,或者负相关。 如何...

qingliangdexiar
2017/05/28
0
0
普通最小二乘法的推导证明(转载)

前言 普通最小二乘法(ordinary least squares, OLS)是线性回归预测问题中一个很重要的概念,在 Introductory Econometrics A Modern Approach (Fourth Edition) 第2章 简单回归模型 中,花...

chensanti234
2017/11/06
0
0
统计-如何读懂线性回归模型的各项指标?

当我们要拟合线性回归模型,有以下三个步骤: 通过最小二乘法,拟合一条直线。 求R^2 求R^2的P值,也就是F值 下面我将一个个分解来讲。 拟合直线:最小二乘法 什么是最小二乘法:每个点(图中...

PriscillaBai
06/09
0
0
机器学习基础篇——最小二乘法

机器学习不是简单的几个算法的累加,如果不能掌握机器学习的一些思考问题的方式以及解决问题的方式,是很难读懂学会这些算法的,也就很难实际去运用这些算法。比如我们前面一直在说的回归问题...

东方神剑
2015/09/12
2.8K
0

没有更多内容

加载失败,请刷新页面

加载更多

70.shell的函数 数组 告警系统需求分析

20.16/20.17 shell中的函数 20.18 shell中的数组 20.19 告警系统需求分析 20.16/20.17 shell中的函数: ~1. 函数就是把一段代码整理到了一个小单元中,并给这个小单元起一个名字,当用到这段...

王鑫linux
今天
2
0
分布式框架spring-session实现session一致性使用问题

前言:项目中使用到spring-session来缓存用户信息,保证服务之间session一致性,但是获取session信息为什么不能再服务层获取? 一、spring-session实现session一致性方式 用户每一次请求都会...

WALK_MAN
今天
5
0
C++ yield()与sleep_for()

C++11 标准库提供了yield()和sleep_for()两个方法。 (1)std::this_thread::yield(): 线程调用该方法时,主动让出CPU,并且不参与CPU的本次调度,从而让其他线程有机会运行。在后续的调度周...

yepanl
今天
4
0
Java并发编程实战(chapter_3)(线程池ThreadPoolExecutor源码分析)

这个系列一直没再写,很多原因,中间经历了换工作,熟悉项目,熟悉新团队等等一系列的事情。并发课题对于Java来说是一个又重要又难的一大块,除非气定神闲、精力满满,否则我本身是不敢随便写...

心中的理想乡
今天
33
0
shell学习之获取用户的输入命令read

在运行脚本的时候,命令行参数是可以传入参数,还有就是在脚本运行过程中需要用户输入参数,比如你想要在脚本运行时问个问题,并等待运行脚本的人来回答。bash shell为此提 供了read命令。 ...

woshixin
今天
4
0

没有更多内容

加载失败,请刷新页面

加载更多

返回顶部
顶部