R软件实现灰色预测GM(1,1)模型

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

灰色预测，即GM模型，是邓聚龙先生在上世纪80年代发明的一种预测方法，突出优点是适用于小样本情况。
因为在小样本条件下，复杂、高端的统计方法往往派不上用场。
但是，GM模型没有被一般的统计软件收录，如SPSS、EViews等等，无法用菜单操作，这给广大使用者带来了一定的困难。
本贴提供了用R软件进行GM(1,1)模型的做法，简单方便，0基础的也容易上手。
这是我三年前给河北省统计系统建模大赛培训时的课件，首次公开。
课件样本：





GM(1,1)课件: GM(1,1)灰色预测.pptx (8.93 MB, 需要: RMB 9 元)

2017-10-4 11:06:22 上传

GM(1,1)课件（R软件操作）
需要: RMB 9 元  [购买]

案例文件： GM(1,1)模型在煤炭需求预测中的应用.pdf (286.29 KB)

2017-10-4 11:07:41 上传

案例文件（免费）

GM(1,1)的R脚本文件： GM11.rar (1.58 KB) 本附件包括：

• GM11.R

2017-10-4 11:08:53 上传

GM(1,1)的R脚本文件

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏14 回帖

关键词：灰色预测 r软件 EVIEWS Eview Views

谢谢楼主分享！

朋友们，代码已经更新，在这里，修正了原程序的一个小bug

GM11.zip (1.65 KB) 本附件包括：

• GM11.R

2022-1-5 23:15:22 上传

完善后的GM11代码

多谢“307662276”朋友支持，问好！

耕耘使者 发表于 2017-10-4 11:10
灰色预测，即GM模型，是邓聚龙先生在上世纪80年代发明的一种预测方法，突出优点是适用于小样本情况。
因为 ...

您好，感谢楼主分享，想知道怎么获取资源呢？我点击链接，显示的是本板块暂不支持直接下载，谢谢！

谁能帅过陈小良1 发表于 2017-11-17 16:00
您好，感谢楼主分享，想知道怎么获取资源呢？我点击链接，显示的是本板块暂不支持直接下载，谢谢！

人大论坛各位朋友好，程序包已经更新，修正了一个小bug,新包见3楼回复中的链接。
如果看不到，请把下面部分代码代替原来的相关部分：
if(p > 0.95){     
cat("小误差概率p>0.95, GM(1,1)预测精度等级为：好",'\n','\n')
}else{
if(p>0.8){
cat("小误差概率p属于(0.8,0.95], GM(1,1)模型预测精度等级为：合格",'\n','\n')
}else{
if(p>0.7){
cat("小误差概率p属于(0.7,0.8], GM(1,1)模型预测精度等级为：勉强合格",'\n','\n')
}else{
cat("小误差概率p<=0.7, GM(1,1)模型预测精度等级为：不合格",'\n','\n')
}
更改说明：上面的两个红色的p是更新的，原代码中是cv.用R软件打开后，手动更改也很方便。

人大论坛各位朋友好，程序包已经更新，修正了一个小bug,新包见3楼回复中的链接。
如果看不到，请把下面部分代码代替原来的相关部分：
if(p > 0.95){     
cat("小误差概率p>0.95, GM(1,1)预测精度等级为：好",'\n','\n')
}else{
if(p>0.8){
cat("小误差概率p属于(0.8,0.95], GM(1,1)模型预测精度等级为：合格",'\n','\n')
}else{
if(p>0.7){
cat("小误差概率p属于(0.7,0.8], GM(1,1)模型预测精度等级为：勉强合格",'\n','\n')
}else{
cat("小误差概率p<=0.7, GM(1,1)模型预测精度等级为：不合格",'\n','\n')
}
更改说明：上面的两个红色的p是更新的，原代码中是cv.用R软件打开后，手动更改也很方便。

人大论坛各位朋友好，程序包已经更新，修正了一个小bug,新包见3楼回复中的链接。
如果看不到，请把下面部分代码代替原来的相关部分：
if(p > 0.95){     
cat("小误差概率p>0.95, GM(1,1)预测精度等级为：好",'\n','\n')
}else{
if(p>0.8){
cat("小误差概率p属于(0.8,0.95], GM(1,1)模型预测精度等级为：合格",'\n','\n')
}else{
if(p>0.7){
cat("小误差概率p属于(0.7,0.8], GM(1,1)模型预测精度等级为：勉强合格",'\n','\n')
}else{
cat("小误差概率p<=0.7, GM(1,1)模型预测精度等级为：不合格",'\n','\n')
}
更改说明：上面的两个红色的p是更新的，原代码中是cv.用R软件打开后，手动更改也很方便。

感谢“伏立特8”朋友支持，新春快乐！

感谢“欧拉的宝宝”朋友支持，春节快乐！