R语言数据分析案例29-基于ARIMA模型的武汉市房价趋势与预测研究
一、选题背景
房地产行业对于国民经济和社会及居民的发展和生活具有很大的影响,而房价能够体现经济运转的好坏,因而房价的波动牵动着开发商和购房者的关注,城市房价预测是一个研究的热点问题,研究房价对民生问题具有重要意义。
本文首先介绍了房地产行业及房价的背景 ,并整理了国内外的相关文献。。。。。
(一)研究背景
房地产行业对于国民经济和社会及居民的发展和生活具有很大的影响,而房价能够体现经济运转的好坏。首先,房地产行业与大众的生活息息相关,比如房产买卖和房屋租赁等等,其发展直接关系到人们的居住条件。
(二)文献综述
国外有许多学者研究了房地产行业的问题。Hekman在1979年综合一些经济因素对房价进行了分析,结果表明经济因素对于房价有显著的影响;Clayton研究了基于波动的理性预期对于房地产价格的影响,然而研究结论与理想的结果并不符合;Normanm Liang对美国两百多个城市的房地产销售价格进行了研究,发现经济条件不同其波动规律不同,并且房价的变化速率的时间段不一样也会存在很大差异。。。。
二、方案论证(设计理念)
RIMA模型被称为自回归移动平均模型,通常当数据序列不是平稳序列时应用,之后通过差分、季节分解等一系列方法后将原本不平稳的序列变为平稳序列,之后利用平稳序列建模方法进行建模。
ARIMA (p ,d, q)模型的特征形式如下:

季节模型
ARIMA可以分为简单的季节模型和乘法季节模型。它是根据季节效应的提取的难易程度进行区分的。当季节效应提取较为容易时,就是简单季节模型,当季节效应提取较为困难时就是乘法季节模型。。。。
三、过程论述
本文所用数据来源于房价官网以及国家统计局官网。对于部分缺失值,采用平均值替代法。
数据集和代码
报告代码数据
数据的描述性统计如表2所示。
表 2 数据的描述性统计
data<-read.table("dataw.csv",header=TRUE,sep=",")
data
price=data$wuhan
summary(price)
| Min. | 1stQu. | Median | Mean | 3rdQu. | Max. | |
| price | 10001 | 15779 | 16492 | 15892 | 17141 | 18581 |
序列的时序图、自相关图及单根检验等可以用来检验序列的平稳性。
应用R作2016年1月-2022年5月武汉市房价的时序图如图1所示。
price<-ts(price,start=c(2016,1),frequency=12)
plot(price,main="2016年1月-2022年5月武汉市房价",xlab="年份",ylab="房价")

从时序图图1中可以清楚的看到该序列蕴涵曲线增长的长期趋势,为非平稳序列。
接下来进行自相关检验,2016年1月-2022年5月武汉市房价自相关图如图2所示。

四、结果分析
根据上面结果可以认为2016年1月到2022年5月武汉市房价序列为非平稳时间序列,不能直接构建ARIMA模型,需要进行差分处理。

#绘制差分后序列自相关图和偏自相关图
acf(price.dif)
pacf(price.dif)
武汉市房价二阶差分时序图和自相关图分别如图5和图6所示,观察到序列已经非常平稳,所有数据均在相同的高度轻微波动。
对二阶差分后的数据进行ADF平稳性检验,结果如表3所示。
表 3 延迟2期平稳性检验
| 检验形式 | no drift no trend | with drift no trend | with drift and trend | |
| Price | ADF统计量 | -6.78 | -6.74 | -6.69 |
| 对应P值 | 0.01 | 0.01 | 0.01 |
#序列的白噪声检验
for(i in 1:2) print(Box.test(price.dif2,lag=6*i))
表 4 白噪声检验
| 滞后期数 | 卡方统计量 | P值 | |
| Price | 滞后6期P值 | 29.972 | 3.979e-05 |
| 滞后12期P值 | 41.721 | 3.71e-05 |
按照上面自相关图和偏自相关图的内容,以及对武汉市房价数据序列进行了二次差分,并结合自动定阶的函数,计算得到模型应该采用ARIMA(1,2,2),拟合得到模型系数如图7所示。
#自动定阶
auto.arima(price.dif2)#模型拟合
price.fit<-arima(price,order=c(1,2,2))
price.fit

进一步观察残差图。
plot(price.fit$residuals,main = "price模型残差图",xlab = "日期",ylab="残差")

最后利用ARIMA(1,2,2)模型对武汉市的房价进行预测,预测6期,即未来六个月房价数据,得到的整体拟合和预测图如下图所示:
从图中可以看到武汉市房价在2022年5月后的六个月呈现出较为平稳的趋势。
五、课程设计总结
本文首先介绍了研究房地产行业及房价的背景,并对国内外相关文献进行了整理。接着利用武汉市2016年1月-2022年五月的房价月度数据,基于时间序列模型,对武汉市房价的动态数据做了时序的分析,并且预测了武汉市六期的房价,并利用三次平滑指数法预测进行对比。。。。。
创作不易,希望大家多点赞关注评论!!!(类似代码或报告定制可以私信)
相关文章:
R语言数据分析案例29-基于ARIMA模型的武汉市房价趋势与预测研究
一、选题背景 房地产行业对于国民经济和社会及居民的发展和生活具有很大的影响,而房价能够体现经济运转的好坏,因而房价的波动牵动着开发商和购房者的关注,城市房价预测是一个研究的热点问题,研究房价对民生问题具有重要意义。 …...
面试-NLP八股文
机器学习 交叉熵损失: L − ( y l o g ( y ^ ) ( 1 − y ) l o g ( 1 − ( y ^ ) ) L-(ylog(\hat{y}) (1-y)log(1-(\hat{y})) L−(ylog(y^)(1−y)log(1−(y^))均方误差: L 1 n ∑ i 1 n ( y i − y ^ i ) 2 L \frac{1}{n}\sum\limits_{i1}^{n}…...
数据仓库之离线数仓
离线数据仓库(Offline Data Warehouse)是一种以批处理方式为主的数据仓库系统,旨在收集、存储和分析大量历史数据。离线数据仓库通常用于定期(如每日、每周、每月)更新数据,以支持各种业务分析、报表生成和…...
Mybatis源码解析
MybatisAutoConfiguration或者MybatisPlusAutoConfiguration核心作用是初始化工厂类SqlSessionFactory,其中包含属性interceptors、MapperLocations、TypeAliasesPackage、TypeEnumsPackage、TypeHandlers等。 MybatisAutoConfiguration自动装配类是由依赖…...
前端学习CSS之神奇的块浮动
在盒子模型的基础上就可以对网页进行设计 不知道盒子模型的可以看前面关于盒子模型的内容 而普通的网页设计具有一定的原始规律,这个原始规律就是文档流 文档流 标签在网页二维平面内默认的一种排序方式,块级标签不管怎么设置都会占一行,而同一行不能放置两个块级标签 行级…...
【Java】内部类、枚举、泛型
目录 1.内部类1.1概述1.2分类1.3匿名内部类(重点) 2.枚举2.1一般枚举2.2抽象枚举2.3应用1:用枚举写单例2.4应用2:标识常量 3.泛型3.1泛型认识3.2泛型原理3.3泛型的定义泛型类泛型接口泛型方法 3.4泛型的注意事项 1.内部类 1.1概述 内部类:指…...
LabVIEW电子类实验虚拟仿真系统
开发了基于LabVIEW开发的电子类实验虚拟仿真实验系统。该系统通过图形化编程方式,实现了复杂电子实验操作的虚拟化,不仅提高了学生的操作熟练度和学习兴趣,而且通过智能评价模块提供即时反馈,促进教学和学习的互动。 项目背景 在…...
SVM支持向量机
SVM的由来和概念 间隔最大化是找最近的那个点的距离’ 之前我们学习的都是线性超平面,现在我们要将超平面变成圈 对于非线性问题升维来解决 对于下图很难处理,我们可以将棍子立起来,然后说不定red跑到左边了,green跑到右边了(可能增加了某种筛选条件导致两个豆子分离)(只是一种…...
【Unity】RPG2D龙城纷争(二)关卡、地块
更新日期:2024年6月12日。 项目源码:后续章节发布 索引 简介地块(Block)一、定义地块类二、地块类型三、地块渲染四、地块索引 关卡(Level)一、定义关卡类二、关卡基础属性三、地块集合四、关卡初始化五、关…...
mediamtx流媒体服务器测试
MediaMTX简介 在web页面中直接播放rtsp视频流,重点推荐:mediamtx,不仅仅是rtsp-CSDN博客 mediamtx github MediaMTX(以前的rtsp-simple-server)是一个现成的和零依赖的实时媒体服务器和媒体代理,允许发布,读取&…...
C# 循环
C# 循环 在编程中,循环是一种控制结构,它允许我们重复执行一段代码多次。C# 提供了几种循环机制,以适应不同的编程需求。本文将详细介绍 C# 中常用的几种循环类型,包括 for 循环、while 循环、do-while 循环和 foreach 循环&…...
PHP杂货铺家庭在线记账理财管理系统源码
家庭在线记帐理财系统,让你对自己的开支了如指掌,图形化界面操作更简单,非常适合家庭理财、记账,系统界面简洁优美,操作直观简单,非常容易上手。 安装说明: 1、上传到网站根目录 2、用phpMyad…...
机器学习中的神经网络重难点!纯干货(上篇)
. . . . . . . . .纯干货 . . . . . . 目录 前馈神经网络 基本原理 公式解释 一个示例 卷积神经网络 基本原理 公式解释 一个示例 循环神经网络 基本原理 公式解释 一个案例 长短时记忆网络 基本原理 公式解释 一个示例 自注意力模型 基本原理…...
[DDR4] DDR1 ~ DDR4 发展史导论
依公知及经验整理,原创保护,禁止转载。 专栏 《深入理解DDR4》 内存和硬盘是电脑的左膀右臂, 挑起存储的大梁。因为内存的存取速度超凡地快, 但内存上的数据掉电又会丢失,一直其中缓存的作用,就像是我们的工…...
享元和代理模式
文章目录 享元模式1.引出享元模式1.展示网站项目需求2.传统方案解决3.问题分析 2.享元模式1.基本介绍2.原理类图3.外部状态和内部状态4.类图5.代码实现1.AbsWebSite.java 抽象的网站2.ConcreteWebSite.java 具体的网站,type属性是内部状态3.WebSiteFactory.java 网站…...
[英语单词] ellipsize,动词化后缀 -ize
openvswitch manual里的一句话:里面有使用ellipsize,但是查字典是没有这个单词,这就是创造出来的动词。将单词ellipsis,加动词化后缀,-ize。 Often we ellipsize arguments not important to the discussion, e.g.: &…...
自然资源-测绘地信专业术语,值得收藏!
自然资源-测绘地信专业术语,值得收藏! 1、1954年北京坐标系 1954年我国决定采用的国家大地坐标系,实质上是由原苏联普尔科沃为原点的1942年坐标系的延伸。 2、1956年黄海高程系统 根据青岛验潮站1950年一1956年的验潮资料计算确定的平均海面…...
如何在小程序中实现页面之间的返回
在小程序中实现页面之间的返回,通常有以下几种方法,这些方法各有特点,适用于不同的场景: 1. 使用wx.navigateBack方法 描述:wx.navigateBack是微信小程序中用于关闭当前页面,返回上一页面或多级页面的API…...
深入解析数据结构之B树:平衡树中的王者
在计算机科学中,数据结构是算法和程序设计的基础。而在众多数据结构中,B树作为一种平衡树,在数据库和文件系统中有着广泛应用。本文将详细介绍B树的概念、特点、操作、优缺点及其应用场景,帮助读者深入理解这一重要的数据结构。 …...
18. 第十八章 继承
18. 继承 和面向对象编程最常相关的语言特性就是继承(inheritance). 继承值得是根据一个现有的类型, 定义一个修改版本的新类的能力. 本章中我会使用几个类来表达扑克牌, 牌组以及扑克牌性, 用于展示继承特性.如果你不玩扑克, 可以在http://wikipedia.org/wiki/Poker里阅读相关…...
Linux应用开发之网络套接字编程(实例篇)
服务端与客户端单连接 服务端代码 #include <sys/socket.h> #include <sys/types.h> #include <netinet/in.h> #include <stdio.h> #include <stdlib.h> #include <string.h> #include <arpa/inet.h> #include <pthread.h> …...
深入浅出Asp.Net Core MVC应用开发系列-AspNetCore中的日志记录
ASP.NET Core 是一个跨平台的开源框架,用于在 Windows、macOS 或 Linux 上生成基于云的新式 Web 应用。 ASP.NET Core 中的日志记录 .NET 通过 ILogger API 支持高性能结构化日志记录,以帮助监视应用程序行为和诊断问题。 可以通过配置不同的记录提供程…...
Android Wi-Fi 连接失败日志分析
1. Android wifi 关键日志总结 (1) Wi-Fi 断开 (CTRL-EVENT-DISCONNECTED reason3) 日志相关部分: 06-05 10:48:40.987 943 943 I wpa_supplicant: wlan0: CTRL-EVENT-DISCONNECTED bssid44:9b:c1:57:a8:90 reason3 locally_generated1解析: CTR…...
Java 语言特性(面试系列2)
一、SQL 基础 1. 复杂查询 (1)连接查询(JOIN) 内连接(INNER JOIN):返回两表匹配的记录。 SELECT e.name, d.dept_name FROM employees e INNER JOIN departments d ON e.dept_id d.dept_id; 左…...
k8s从入门到放弃之Ingress七层负载
k8s从入门到放弃之Ingress七层负载 在Kubernetes(简称K8s)中,Ingress是一个API对象,它允许你定义如何从集群外部访问集群内部的服务。Ingress可以提供负载均衡、SSL终结和基于名称的虚拟主机等功能。通过Ingress,你可…...
IGP(Interior Gateway Protocol,内部网关协议)
IGP(Interior Gateway Protocol,内部网关协议) 是一种用于在一个自治系统(AS)内部传递路由信息的路由协议,主要用于在一个组织或机构的内部网络中决定数据包的最佳路径。与用于自治系统之间通信的 EGP&…...
深入理解JavaScript设计模式之单例模式
目录 什么是单例模式为什么需要单例模式常见应用场景包括 单例模式实现透明单例模式实现不透明单例模式用代理实现单例模式javaScript中的单例模式使用命名空间使用闭包封装私有变量 惰性单例通用的惰性单例 结语 什么是单例模式 单例模式(Singleton Pattern&#…...
Cinnamon修改面板小工具图标
Cinnamon开始菜单-CSDN博客 设置模块都是做好的,比GNOME简单得多! 在 applet.js 里增加 const Settings imports.ui.settings;this.settings new Settings.AppletSettings(this, HTYMenusonichy, instance_id); this.settings.bind(menu-icon, menu…...
MySQL 8.0 OCP 英文题库解析(十三)
Oracle 为庆祝 MySQL 30 周年,截止到 2025.07.31 之前。所有人均可以免费考取原价245美元的MySQL OCP 认证。 从今天开始,将英文题库免费公布出来,并进行解析,帮助大家在一个月之内轻松通过OCP认证。 本期公布试题111~120 试题1…...
Python Ovito统计金刚石结构数量
大家好,我是小马老师。 本文介绍python ovito方法统计金刚石结构的方法。 Ovito Identify diamond structure命令可以识别和统计金刚石结构,但是无法直接输出结构的变化情况。 本文使用python调用ovito包的方法,可以持续统计各步的金刚石结构,具体代码如下: from ovito…...
