回归分析(2) 一元回归模型

news2025/1/27 12:35:35

        如上所述,为了易于确定回归函数μ(x)中的未知参数,我们来讨论变量Y与x之间存在着线性相关关系的情形

 散布在某一条直线的周围.于是,我们可以用线性回归方程

 来描述Y与x之间的相关关系,并假设相应的误差(称为随机误差)

 服从正态分布N(0,\sigma^2),其中a,b,\sigma^2 都是不依赖于x的未知参数,参数b 称为回归系数

        对Y作这样的假设, 相当于假设对于x (在某个区间内)的每一个值,都有

 于是对x 与Y的每一组独立试验(或样本)

 也有

 式(5. 8)或式(5.10)称为一元(正态)线性回归模型.若将试验观测值(5.6)代入,则
可获得随机误差\varepsilon _i的一组观测值

 对于一元线性回归模型,需要研究下面几个问题:

(1)根据样本(5. 9)估计参数a,b,\sigma^2,从而建立线性回归模型;

(2)线性相关关系的显著性检验;

(3)利用所获得的线性回归模型对变量Y的取值进行预测.

1.参数a,b 的无偏估计及其分布

对于一元线性回归模型(5. 10),这时的离差平方和(5.4)为

为了使s取得最小值,按最小二乘法,分别求s对a及b的偏导,有

 令这两个偏导数为零,得

 

 整理,得

 或者

 其中

 分别为样本的各种矩,方程组(5. 13) 或方程组(5. 14) 称为正规方程(组).解此正规
方程可得a与b的估计值(分别记作\widehat{a}\widehat{b})

 其中

x_1,x_2,...,x_n的二阶中心矩.

``       与式(5.16)相对应的估计量可表示为

 

 作为正规方程的解,估计量a与b自然应满足正规方程(5. 13)和方程组(5.14),即

 或

 这时残差平方和(即s的最小值)可表示为

 将通过式(5. 16)计算得到的a 与b 的值代入方程(5.7),就可得到Y对x 的(经验)线性回归方程

 和相应的回归值(或拟合值)

 方程(5.22)的图形称为(经验)线性回归直线.

 

 该式表明,经验回归直线过由点集(5.6)构成的散点图的几何中心(\overline{x},\overline{y})

        下面讨论估计量\widehat{b}的分布.由估计量(5.18)的第一式知

所以\widehat{b}是随机变量Y_1,Y_2,..,Y_n的线性组合.再注意到x_i(i = 1,2,...,n)的非随机
性及Y1 ,Y2, …,Y"是独立的正态随机变量,且

 

 则可知\widehat{b}也服从正态分布,且有

 故\widehat{b}是b的无偏估计.同样,我们可算得

 

 这样,我们已经证明并获得了

定理5.2.1 在上述记号下,估计量\widehat{b}是b的无偏估计,且

 其中

 同理可证

是a的无偏估计(见习题5.1,由读者自己证明),并且有下述定理.

定理5.2.2 估计量\widehat{a}是a的无偏估计,且

5.2.2 参数\sigma^2的无偏估计及其分布

        为了获得\sigma^2的无偏估计,我们作离差分解

 

 由式(5. 19)知上式的最后两个和式均为0,故有

 其中右端第一项正是式(5. 21)定义的残差平方和

于是我们得到

 其中右端三项的自由度分别为又由于

 且相互独立,故由定理1. 2. 2可得如下的定理.

 定理5.2.3 在上述记号下, S_E,\overline{Y},\widehat{b}相互独立且

 由于卡方分布的期望是它的自由度,故由式(5. 32)

从而

 因此

 是\sigma^2的无偏估计.在此记号下,定理5.2. 3又可改写成如下形式.

2. 一元线性回归参数的计算

为了计算无偏估计\sigma^2,将

 代入式(5.33),得

 再注意到式(5.18)的第一式,有

 并引进计算器容易获得其值的二阶中心矩的记号

 则有

 式(5.38)中的\sigma^2,\widehat{a},\widehat{b},也可利用表5.1及计算器的统计功能方便地获得.

例题

 

 

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.coloradmin.cn/o/39950.html

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈,一经查实,立即删除!

相关文章

Go-zero框架学习+xorm+gorm配置

Go-zero框架学习xormgorm配置 文章目录Go-zero框架学习xormgorm配置一、框架介绍二、go-zero快速搭建1.下载go-zero2.安装goctl3.项目初始化4.测试5.项目结构6.快速生成rpc服务7.快速添加api接口8.快速生成model服务8.快速生成Dockerfile9.快速生成K8s部署文件三.golang的ORM框…

[go学习笔记.第十七章.redis的使用] 1.redis的使用

1.redis基本介绍 (1).Redis 是 NoSQL 数据库,不是传统的关系型数据库,官网: https://redis.io/ 和http://redis.cn/ (2).Redis: REmote Dlctionary Sever(远程字典服务器), Redis 性能非常高,单机能够达到 15w qps,通常适合做缓存…

刷爆力扣之盛最多水的容器

刷爆力扣之盛最多水的容器 HELLO,各位看官大大好,我是阿呆 🙈🙈🙈 今天阿呆继续记录下力扣刷题过程,收录在专栏算法中 😜😜😜 该专栏按照不同类别标签进行刷题&#xff…

安卓的分区一点有用知识:super、lpunpack、lpdump

我们知道这个安卓的镜像分区有很多个。 那么这个文章要介绍什么呢? 三个点: 一是现在的android支持动态分区,很多的东西都被放到super分区里面了,这个应该是可以配置的。然后super里面有比如system、vendor这种比较大的分区。那…

教务排课系统毕业设计,大学排课系统设计与实现,排课系统论文作品参考

功能清单 【后台管理员功能】 录入分院:录入分院名称,简介,详情 分院管理:管理已经录入分院,支持修改和删除 老师录入:录入老师姓名、联系方式、头像、老师简介 老师管理:管理所有已经录入老师…

基于Matlab模拟用于海况海洋学研究的 X 波段雷达系统(附源码)

目录 一、定义雷达系统参数 二、对海面进行建模 三、配置雷达收发器 四、生成数据多维数据集 五、处理海面回波 六、总结 七、程序 海事雷达系统在充满挑战的动态环境中运行。为了改进对感兴趣目标的检测并评估系统性能,必须了解海面返回的性质。 在本例中&a…

【操作系统】2.4 死锁

这一节也非常重要 2.4.1 死锁的概念 2.4.1 死锁的概念_StudyWinter的博客-CSDN博客 在并发环境下,各种进程因竞争资源而造成的一种互相等待对方手里的资源,导致各进程都阻塞,都无法向前推进的现象。这就是死锁,死锁发生后&#…

Jest API使用方法

如上面的知识图谱所示,一个常见的测试框架通常需要实现这些功能: ● before/after 钩子函数: 如beforeEach,afterEach, ● Mock方法: 函数Mock,时间mock等。 ● 断言: 判断一个描述是否正确,在Jest中常为 e…

你心心念念的RabbitMQ个人实践来了来了它来了

前言 MQ(Message Queue)就是消息队列,其有点有很多:解耦、异步、削峰等等,本文来聊一下RabbitMQ的一些概念以及使用。 RabbitMq 案例 Springboot整合RabbitMQ简单案例 基本概念 Exchange:消息交换机&a…

云原生系列 【基于CCE Kubernetes编排实战二】

✅作者简介: CSDN内容合伙人,全栈领域新星创作者,阿里云专家博主,阿里云问答板块版主,华为云享专家博主,掘金后端评审团成员 💕前言: 最近云原生领域热火朝天,那么云原生…

Hystirx限流:信号量隔离和线程池隔离

背景: 最近工作中要处理服务高并发的问题,大流量场景下限流熔断降级可以说是必不可少的,打算对限流做一次改造,所以要先了解一下hytrix相关内容,比如了解一下线程池隔离和信号量隔离的区别。 **信号量:**信…

[网络工程师]-应用层协议-DHCP

BOOTP是最早的主机配置协议,动态主机配置协议(Dynamic Host Configuration Protocol,DHCP)则是在其基础上进行了改良的协议,是一种用于简化主机IP配置管理的IP管理标准。通过DHCP协议,DHCP服务器为DHCP客户…

集合学习笔记——Collection 全家桶

Collection是我们日常开发中使用频率非常高的集合,它的主要实现有List和Set,区别是List是有序的,元素可以重复;Set是无序的,元素不可以重复,我们简单看下继承关系: List的实现类主要线程不安全的ArrayList和LinkedList…

推挽输出和开漏输出-三极管-mos管

一、推挽输出 1.1推挽输出的概念 推挽(push-pull)输出是由两个MOS或者三极管组成,两个管子始终保持一个导通,另一个截止的状态。 图1 推挽电路示意图 当输入高电平时,叫做推; 上管Q1导通,下管…

【目标检测】Faster R-CNN论文的讲解

目录:Faster R-CNN论文的讲解一、前言二、回顾Fast R-CNN三、引入Faster R-CNN四、Faster R-CNN的介绍4.1 框架结构4.2 RPN如何产生候选区域的4.3 损失函数4.4 训练候选框提取网络4.5 RPN和Fast R-CNN共享特征的方法4.5.1 交替训练法4.5.2 近似联合训练法一、前言 …

C语言——学生信息管理系统

目录 功能展示 界面展示 所有功能模块: 功能1:菜单模块(显示功能菜单) 功能2:增加学生信息 功能3:输出学生信息(查看所有学习信息) 功能4:修改学生信息 功能5&a…

python3-GUI概述及应用

目录一、什么是GUI二、Python GUIPySimpleGUI概述一、PySimpleGUI简介二、PySimpleGUI特征三、输出设备hello,world猜数字一、玩家猜数字二、电脑猜数字21点游戏一、21点游戏简介二、程序代码一、什么是GUI 图形用户界面(Graphical User Interface,简称…

十六、CANdelaStudio深入-CDD与CDDT的差异(新建自定义服务)

本专栏将由浅入深的展开诊断实际开发与测试的数据库编辑,包含大量实际开发过程中的步骤、使用技巧与少量对Autosar标准的解读。希望能对大家有所帮助,与大家共同成长,早日成为一名车载诊断、通信全栈工程师。 本文介绍CANdelaStudio的CDD与CDDT的差异与新建自定义服务,欢迎…

数字图像处理(一)——什么是数字图像

一、什么是数字图像处理? 一副图像可以被定义为一个二维函数f(x,y),其中x和y是空间平面坐标,而对任意一对空间坐标(x,y)处幅值f称为图像在该点的强度或者灰度。当x和y以及灰度值f是有限的离散数值时,我们称该图像为数字图像。像素…

排序算法简述

一、概述 常见的排序算法有冒泡排序、插入排序、选择排序、快速排序、归并排序、桶排序、基数排序,这些排序各自有各自的特点。按照时间时间复杂度可以分为 O(n^2):冒泡、插入、选择排序;O(nlogn):归并、快速排序;O(n):桶排序、计数排序、基…