楼主: xuelida
8049 15

[分享]我自己编的小程序:matlab的DF-GLS的单位根检验功效模拟程序 [推广有奖]

  • 1关注
  • 94粉丝

教师

计量经济学

已卖:5445份资源

学科带头人

50%

还不是VIP/贵宾

-

威望
1
论坛币
58865 个
通用积分
30.5491
学术水平
394 点
热心指数
256 点
信用等级
169 点
经验
47973 点
帖子
1506
精华
7
在线时间
1083 小时
注册时间
2005-12-5
最后登录
2025-12-2

楼主
xuelida 在职认证  发表于 2007-7-3 15:59:00 |AI写论文

+2 论坛币
k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群
赵安豆老师微信:zhaoandou666

经管之家联合CDA

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

感谢您参与论坛问题回答

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

+2 论坛币

一、数据生成过程:DGP1、yt=βyt-1t, DGP2、yt= ut,ut= βut-1+ εt ,其中,α、δ = 0.5,β = 0.95下同。

这时共4个程序, 程序Untitled53.m的DGP1,GLS去势p=1,a=-7时,模拟的样本可以自己选定,根据不同的样本计算不同的概率,这时的概率计算的临界值由我给的上一个下载里的https://bbs.pinggu.org/thread-198766-1-1.html&page=3

“matlab的DF-GLS的分位数为5%时的单侧临界值模拟程序”

程序Untitled54.m的DGP1,GLS去势p=2,a=-13时。

程序Untitled73.m的DGP2,GLS去势p=1,a=-7时。

程序Untitled74.m的DGP2,GLS去势p=2,a=-13时。

二、数据生成过程:DGP1、yt = α +β yt-1+ εt DGP2、yt = α + ut,ut = βut-1+ εt

这时共4个程序, 程序Untitled91.m的DGP1,GLS去势p=1,a=-7时,模拟的样本可以自己选定,根据不同的样本计算不同的概率,这时的概率计算的临界值由我给的上一个下载里的https://bbs.pinggu.org/thread-198766-1-1.html&page=3

“matlab的DF-GLS的分位数为5%时的单侧临界值模拟程序”

程序Untitled92.m的DGP1,GLS去势p=2,a=-13时。

程序Untitled93.m的DGP2,GLS去势p=1,a=-7时。

程序Untitled94.m的DGP2,GLS去势p=2,a=-13时。

三、数据生成过程:DGP1、yt = α + δ t+βyt-1t, DGP2、yt = α + δ t +ut,ut = βut-1+ εt

这时共4个程序, 程序Untitled41.m的DGP1,GLS去势p=1,a=-7时,模拟的样本可以自己选定,根据不同的样本计算不同的概率,这时的概率计算的临界值由我给的上一个下载里的https://bbs.pinggu.org/thread-198766-1-1.html&page=3

“matlab的DF-GLS的分位数为5%时的单侧临界值模拟程序”

程序Untitled42.m的DGP1,GLS去势p=2,a=-13时。

程序Untitled43.m的DGP2,GLS去势p=1,a=-7时。

程序Untitled44.m的DGP2,GLS去势p=2,a=-13时。

132564.rar (6.42 KB, 需要: 6 个论坛币)
二维码

扫码加我 拉你入群

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

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

关键词:MATLAB atlab 单位根检验 matla 模拟程序 MATLAB 检验 单位 模拟程序 功效

已有 2 人评分经验 论坛币 学术水平 收起 理由
Xaero + 20 + 20 + 1 我很赞同
xuehe + 30 tk

总评分: 经验 + 20  论坛币 + 50  学术水平 + 1   查看全部评分

沙发
hhssbb(真实交易用户) 发表于 2007-8-7 22:08:00

xuelida你好,想请教一下你程序中的几个小问题:(以untitled42.m为例)

1、第52行 k(i)=T*(q-1) %这句命令起到什么作用呢?q相当于et=q*et-1+u中的系数,T为样本容量。

2、第54-56行 if tau(i)<-2.9121
astau5=astau5+1; %why?
end

3、检验的功效 最终由哪个公式计算得出?功效值在程序中如何体现呢?

以上是我不大明白的地方。请指教。

藤椅
xuelida(未真实交易用户) 在职认证  发表于 2007-8-8 01:27:00

q不就是单位根检验回归参数值,这个就是DF单位根检验统计量呀,k(i)=T*(q-1)

如果tau(i)<-2.9121
把astau5加起来:astau5=astau5+1;

然后astau5/n就是就算的检验功效了

板凳
hhssbb(真实交易用户) 发表于 2007-8-13 10:53:00

明白了,谢谢你的指点!!!

报纸
hhssbb(真实交易用户) 发表于 2007-8-13 11:27:00

untitled42.m给出的功效值是0啊。

另外,个人觉得将临界值-2.9121直接表达为相应分布函数的反函数更直观。

地板
xuelida(未真实交易用户) 在职认证  发表于 2007-8-13 11:48:00

这个看什么形式更直观,或者计算速度更快,更好理解,

我也不太会编程,有些很繁琐,希望各位同仁批评修改

7
zhaomn200145(未真实交易用户) 发表于 2007-8-13 11:49:00
楼主啊,你有Gauss版的单位根检验功效模拟程序吗?我的机器跑不了matlab,只能跑gauss啊,我现在正在学习这个检验功效模拟程序,卖贵一点我也可以啊,多谢多谢。

8
lyrasover(未真实交易用户) 发表于 2007-8-13 12:43:00
不错的

9
xuelida(未真实交易用户) 在职认证  发表于 2007-8-14 13:30:00
gauss版本的还没有写出,等有时间我会贴出

10
chen251213(真实交易用户) 发表于 2007-10-30 17:41:00
呵呵,学习学习!

您需要登录后才可以回帖 登录 | 我要注册

本版微信群
加好友,备注cda
拉您进交流群
GMT+8, 2026-1-18 02:28