R语言 线性回归分析实例

y,X1,X2,X3 分别表示第 t 年各项税收收入(亿元),某国生产总值GDP(亿元),财政支出(亿元)和商品零售价格指数(%).

(1) 建立线性模型

 ① 自己编写函数:

> library(openxlsx)
> data = read.xlsx("22_data.xlsx",sheet = 1)
> x = data[,-c(1,2)]
> x = cbind(rep(1,17),x)
> x_mat  = as.matrix(x)
> y =matrix(data[,2],ncol = 1)
> res = solve(t(x_mat)%*%x_mat)%*%t(x_mat)%*%y
> res[,1]
rep(1, 17) 19412.8597818
X1             0.2679605
X2            -0.2874013
X3          -297.3653736

 所以各参数的估计值分别为 

 ② lm函数

> lm(y~x_mat)Call:
lm(formula = y ~ x_mat)Coefficients:(Intercept)  x_matrep(1, 17)          x_matX1  
19412.859781545               NA      0.267960511  x_matX2          x_matX3  -0.287401287   -297.365373557  

 于是各参数的估计值分别为

这两个方法的结果是一样的。

 

(2)要求实验报告中画出矩阵散点图,给出参数的点估计、区间估计、t检验值、判定系数和模型F检验的方差分析表

绘制矩阵散点图。

library(graphics)
pairs(data[,-1]pch = 21,bg = c('red','green3','blue'))
# pch参数是控制点的形状,bg是控制点的颜色

 

下面代码给出参数的点估计,t检验值,判定系数

> summary(lm(y~x_mat+1))Call:
lm(formula = y ~ x_mat + 1)  #调用Residuals:      #残差统计量,残差第一四分位数(1Q)和第三分位数(3Q)有大约相同的幅度,意味着有较对称的钟形分布Min      1Q  Median      3Q     Max  
-4397.9 -1102.4   153.8  1184.4  2934.6 Coefficients: (1 not defined because of singularities)        Estimate Std. Error t value Pr(>|t|)    
(Intercept)      1.941e+04  3.524e+04   0.551    0.591    
x_matrep(1, 17)         NA         NA      NA       NA    
x_matX1          2.680e-01  4.466e-02   6.000 4.45e-05 ***
x_matX2         -2.874e-01  1.668e-01  -1.723    0.109    
x_matX3         -2.974e+02  3.688e+02  -0.806    0.435    
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
#标记为Estimate的列包含由最小二乘法计算出来的估计回归系数。
#标记为Std.Error的列是估计的回归系数的标准误差。
#从理论上说,如果一个变量的系数是0,那么该变量将毫无贡献。然而,这里显示的系数只是估计,它们不会正好为0.
#因此,我们不禁会问:从统计的角度而言,真正的系数为0的可能性有多大?这是t统计量和P值的目的,在汇总中被标记为t value和Pr(>|t|) #P值估计系数不显著的可能性,有较大P值的变量是可以从模型中移除的候选变量

Residual standard error: 2013 on 13 degrees of freedom Multiple R-squared: 0.9982, Adjusted R-squared: 0.9977 F-statistic: 2348 on 3 and 13 DF, p-value: < 2.2e-16

#Residual standard error 表示残差的标准差,F-statistic 表示F的统计量

  区间估计?方差分析表?

 

(3)保留模型中线性关系显著的预测变量确定最后的模型,并利用R软件中的"predict"语句预测2017年的税收收入

 根据回归分析结果,只有变量X1具有显著性。所以模型中仅保留变量X1。

构造模型:

x_mat = cbind(rep(1,17),data[,3]) 
y = data[,2]
res = lm(y~x_mat)
res
> resCall:
lm(formula = y ~ x_mat)Coefficients:
(Intercept)       x_mat1       x_mat2  -6213.0189           NA       0.1915 

 

 该模型为:Y = -6213.0189 + 0.1915 X1

接下来预测2017年的税收收入,先根据数据data对 t 和 y 之间的关系进行回归分析

t = data[,1]
y = data[,2]
res = lm(y~t)
res> resCall:
lm(formula = y ~ t)Coefficients:
(Intercept)            t  -16428607         8213  

 所以 t 与 y 的关系为:y  = -16428607 + 8213 t

 预测 2017 年的税收收入:

> newdata = data.frame(t = 2017)
> pre = predict(res,newdata,interval = "prediction",level = 0.95)
> prefit      lwr      upr
1 136337.8 116018.1 156657.4

  

 

转载于:https://www.cnblogs.com/jiaxinwei/p/11523752.html

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

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

相关文章

【Cef编译】 CefSharp编译失败,检测到“RuntimeLibrary”的不匹配项: 值“MT_StaticRelease”不匹配值“MD_DynamicRelease”...

编译CefSharp生成后一个libcef_dll_wrapper.lib时&#xff0c;供CefSharp使用。结果CefSharp编译的时候报错。遇到以下异常&#xff1a;libcef_dll_wrapper.lib(binary_value_ctocpp.obj) : error LNK2038: 检测到“RuntimeLibrary”的不匹配项: 值“MT_StaticRelease”不匹配值…

CEF编译 执行gn args out\Release_GN_x86异常

gn args out\Debug_GN_x86 用来配置编译参数&#xff0c;执行gn args out\Release_GN_x86时异常&#xff1a; Toolchain is out of date. Run "gclient runhooks" to update the toolchain, or set DEPOT_TOOLS_WIN_TOOLCHAIN0 to use the locally installed toolcha…

拜托!面试请不要再问我 Spring Cloud 底层原理 ...

前些天发现了一个巨牛的人工智能学习网站&#xff0c;通俗易懂&#xff0c;风趣幽默&#xff0c;忍不住分享一下给大家。点击跳转到教程。 >出处&#xff1a; >https://www.fangzhipeng.com > 本文出自[方志朋的博客](http://blog.csdn.net/forezp) 本文为转载文章&…

springCloud - 第1篇 - 服务的注册 Eureka

前些天发现了一个巨牛的人工智能学习网站&#xff0c;通俗易懂&#xff0c;风趣幽默&#xff0c;忍不住分享一下给大家。点击跳转到教程。 PS&#xff1a;这个系列不定时更新&#xff0c;只是个人的学习分享&#xff0c; 内容全程参考书目&#xff1a; 《Spring Cloud 与 Do…

解决:com.mysql.jdbc.MysqlDataTruncation: Data truncation: Data too long for column ‘ip‘ at row 1

前些天发现了一个巨牛的人工智能学习网站&#xff0c;通俗易懂&#xff0c;风趣幽默&#xff0c;忍不住分享一下给大家。点击跳转到教程。 1. 报错&#xff1a; SQL []; Data truncation: Data too long for column ip at row 1; nested exception is com.mysql.jdbc.MysqlDa…

Xampp配置本地域名及常见错误解决

Xampp配置本地域名及常见错误解决 本地域名配置1、计算机-->C盘-->Windows-->System32-->drivers-->etc-->hosts127.0.0.1 localhost//设置你要配置的本地域名2、计算机-->XAMPP-->apache-->conf-->extra-->httpd-vhosts.conf<Virtu…

解决 Idea 卡在 Resolving Maven dependencies ...

前些天发现了一个巨牛的人工智能学习网站&#xff0c;通俗易懂&#xff0c;风趣幽默&#xff0c;忍不住分享一下给大家。点击跳转到教程。 Idea卡在Resolving Maven dependencies的解决方案 在Reimpot All Maven Porjects时, 如果项目过大, maven依赖过多, 会直接卡在Resolvin…

VS Code (visual studio code) VSC 编辑器(微软出品,js开发的编辑器)

一.选择合适的编辑器&#xff0c;提高编程效率 代码编辑器的选择&#xff0c;可以说是开发者社区中一个经久不衰的话题&#xff0c;现今编辑器的数量数不胜数&#xff0c;vim&#xff0c;sublime Text,Emacs,Atom等等&#xff0c;那么对于一个开发者而言&#xff0c;挑选一个合…

Docker 安装 redis 、Redis docker 方式部署

前些天发现了一个巨牛的人工智能学习网站&#xff0c;通俗易懂&#xff0c;风趣幽默&#xff0c;忍不住分享一下给大家。点击跳转到教程。 1. 找镜像&#xff1a; docker search redis2. 拉取镜像&#xff1a; 在这一步可以选择版本&#xff0c;不选择版本&#xff0c;默认为…

xcode windows版安装使用教程

随着iPhone、iPad、Mac等苹果产品越来越火爆&#xff0c;越来越多的初学者想要了解和尝试苹果平台&#xff0c;包括苹果操作系统Mac OS X、苹果演示软件Keynote、苹果开发工具Xcode等。然而&#xff0c;苹果电脑价格昂贵&#xff0c;并不是每个人都可以承受。 因此&#xff0c;…

解决:Unable to open debugger port (127.0.0.1:55017): java.net.SocketException “Socket closed“

项目以前启动正常&#xff0c;突然报错&#xff0c;启动不起来了&#xff0c;报了个Unable to open debugger port (127.0.0.1:55017): java.net.SocketException "Socket closed"这个错。 前些天发现了一个巨牛的人工智能学习网站&#xff0c;通俗易懂&#xff0c;…

html5 如何打包成apk,将H5封装成android应用APK文件的几种方法

直接使用编程软件提供的方法&#xff1a; 1、需要下载安装MyEclipse2014&#xff0c;Android SDK&#xff0c;eclipse(需配置Android开发环境) Java和Android环境安装与配置。 2、打开MyEclipse2014&#xff0c;新建一个HTML5 Mobile Application Project&#xff0c;命名&…

解决 Unmapped Spring configuration files found.Please configure Spring facet.

最近在学习使用IDEA工具&#xff0c;觉得与Eclipse相比&#xff0c;还是有很多的方便之处。 但是&#xff0c;当把自己的一个项目导入IDEA之后&#xff0c;Event Log提示“Unmapped Spring configuration files found.Please configure Spring facet.” 这个提示不影响工程正…

uni-app—从安装到卸载

uni-app实现了一套代码&#xff0c;同时运行到多个平台。支持iOS模拟器、Android模拟器、H5、微信开发者工具、支付宝小程序Studio、百度开发者工具、字节跳动开发者工具 工具安装 开发uni-app需要安装HBuilder X.下载地址。 下载成功后直接解压即可 简单的配置一下开发偏好&am…

jenkins 执行构建 并查看结果

前些天发现了一个巨牛的人工智能学习网站&#xff0c;通俗易懂&#xff0c;风趣幽默&#xff0c;忍不住分享一下给大家。点击跳转到教程。 继完成构建项目配置http://www.cnblogs.com/yajing-zh/p/5111060.html后&#xff0c;则要执行构建。 回到jenkins主页之后&#xff0c;…

使用jquery获取url及url参数的方法

使用jquery获取url以及使用jquery获取url参数是我们经常要用到的操作 1、jquery获取url很简单&#xff0c;代码如下&#xff1a; window.location.href; 其实只是用到了javascript的基础的window对象&#xff0c;并没有用jquery的知识。 2、jquery获取url参数比较复杂&#xff…

js 遍历对象的几种方法

第一种&#xff1a; for......in const obj {id:1,name:zhangsan,age:18}for(let key in obj){console.log(key --- obj[key])} 输出结果&#xff1a; 第二种&#xff1a; 1&#xff09;、Object.keys&#xff08;obj&#xff09; 2&#xff09;、Object.values&#xff08…

简述 maven 命令 package、install、deploy 的区别

前些天发现了一个巨牛的人工智能学习网站&#xff0c;通俗易懂&#xff0c;风趣幽默&#xff0c;忍不住分享一下给大家。点击跳转到教程。 原文有图片说明命令运行过程&#xff0c;我只记录最终结论&#xff1a; 区别即&#xff1a; package &#xff1a;打包&#xff08;到项…

教给孩子的10句“保命金言”(图)

1、平安成长比成功更重要。 教育孩子人人有若干权利&#xff0c;如呼吸权、生命权、隐私权&#xff0c;这些权利任何人不能剥夺。告诉孩子&#xff0c;任何人也无权剥夺他的安全权&#xff0c;安全重于一切。 2、背心裤衩覆盖的地方不许别人摸。 孩子应当知道身体属于自己&a…

Linux 关闭、开启、查看 防火墙命令

前些天发现了一个巨牛的人工智能学习网站&#xff0c;通俗易懂&#xff0c;风趣幽默&#xff0c;忍不住分享一下给大家。点击跳转到教程。 在外部访问CentOS中部署应用时&#xff0c;需要关闭防火墙。 关闭防火墙命令&#xff1a;systemctl stop firewalld.service 开启防火…