全部版块 我的主页
论坛 数据科学与人工智能 数据分析与数据科学 R语言论坛
2012-4-15 14:53:51
雁茗轩 发表于 2012-4-15 14:48
真是不好意思,太麻烦你了老师!
老师您太厉害了,我这回可以估计出来统计量了!真的太感激您了。您不愧学术精湛,而且毫不吝啬赐教,真的深受感动!
二维码

扫码加我 拉你入群

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

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

2012-4-17 19:11:27
雁茗轩 发表于 2012-4-15 14:53
老师您太厉害了,我这回可以估计出来统计量了!真的太感激您了。您不愧学术精湛,而且毫不吝啬赐教,真的 ...
epoh老师您好,我想再对所建的模型进行一个脉冲响应分析,不知道是我理解的不到位还是怎么回事,我估计的时候老出错,出来了所需要的图的标题,还出不来曲线图,这是数据和程序,您有时间的话帮我看看,谢谢了!
附件列表

d.txt

大小:2.29 KB

 马上下载

二维码

扫码加我 拉你入群

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

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

2012-4-17 20:44:14
雁茗轩 发表于 2012-4-17 19:11
epoh老师您好,我想再对所建的模型进行一个脉冲响应分析,不知道是我理解的不到位还是怎么回事,我估计的 ...

确认一下thVar的范围
yan=read.table("d.txt",header=TRUE)
y=yan$CPI
x=yan$do1
summary(y)
#   Min. 1st Qu.  Median    Mean 3rd Qu.    Max.
#-1.3000 -0.2000  0.2000  0.2369  0.6875  2.6000
summary(x)
#    Min.  1st Qu.   Median     Mean  3rd Qu.     Max.
#-27.5000  -1.4830   1.1250   0.6761   3.6520  14.2800

二维码

扫码加我 拉你入群

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

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

2012-4-17 21:26:12
雁茗轩 发表于 2012-4-17 19:11
epoh老师您好,我想再对所建的模型进行一个脉冲响应分析,不知道是我理解的不到位还是怎么回事,我估计的 ...
thVar=yan$do1,图形不对
若改为:mod=lstar(y,m=3,d=1,control=list(maxit=3000))
图形如下:
   yan.bmp
二维码

扫码加我 拉你入群

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

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

2012-4-17 22:12:53
epoh 发表于 2012-4-17 21:26
thVar=yan$do1,图形不对
若改为:mod=lstar(y,m=3,d=1,control=list(maxit=3000))
图形如下:
好的,太感谢您了,我再重新做一下!
二维码

扫码加我 拉你入群

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

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

2012-4-17 23:19:42
epoh 发表于 2012-4-17 21:26
thVar=yan$do1,图形不对
若改为:mod=lstar(y,m=3,d=1,control=list(maxit=3000))
图形如下:
老师不好意思,老麻烦您,还是有问题,我确认了thVar的范围,像您183楼所说,但是结果还是那样子;
若把模型改为:mod=lstar(y,m=3,d=1,control=list(maxit=3000))这个的话,是不是相当于只是Y自己回归的转换函数的脉冲了呢?谢谢您了啊!
二维码

扫码加我 拉你入群

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

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

2012-4-18 10:16:35
雁茗轩 发表于 2012-4-17 23:19
老师不好意思,老麻烦您,还是有问题,我确认了thVar的范围,像您183楼所说,但是结果还是那样子;
若把 ...
你确认了thVar没错,
模型画出的 transition functions,也可以
   tf.bmp
可是girf图形却是....,这个头疼,得花时间查查

二维码

扫码加我 拉你入群

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

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

2012-4-18 12:01:18
epoh 发表于 2012-4-18 10:16
你确认了thVar没错,
模型画出的 transition functions,也可以
是啊,不知道哪里出了问题……
二维码

扫码加我 拉你入群

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

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

2012-4-18 14:21:12
雁茗轩 发表于 2012-4-18 12:01
是啊,不知道哪里出了问题……
呵呵太久没用忘了
这个方法maybe不适用external threshold variable

二维码

扫码加我 拉你入群

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

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

2012-4-19 22:17:05
epoh 发表于 2012-4-18 14:21
呵呵太久没用忘了
这个方法不适用external threshold variable
因为bootstrap 过程,threshold variable ...
哦,好的,那就不做了!真心谢谢epoh老师了,老麻烦您!
二维码

扫码加我 拉你入群

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

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

2012-4-21 23:34:45
雁茗轩 发表于 2012-4-19 22:17
哦,好的,那就不做了!真心谢谢epoh老师了,老麻烦您!
epoh老师:
         您好!
         还是关于那个LSTAR模型中的脉冲响应函数的,那我那个模型中,已经估计出来了伽马和c,是不是可以按照泰勒三阶展开了呢?然后对这个线性的模型进行脉冲影响估计?这个能实现吗或者说有意义没?我给您看个文章,这个也是在LSTAR模型中做了脉冲,但是她分的两种环境,我不知道这个是具体怎么分析方法,您看能对解决我的问题有帮助没?谢谢您了,不好意思麻烦您,但是真的很需要您的帮助,谢谢啦!
附件列表
二维码

扫码加我 拉你入群

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

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

2012-4-22 08:29:08
雁茗轩 发表于 2012-4-21 23:34
epoh老师:
         您好!
         还是关于那个LSTAR模型中的脉冲响应函数的,那我那个模型中,已经 ...
这篇文献是LSTVAR(4)
双变量: 产出缺口,通货膨胀      
转换变量(St),选滞后二期通胀

你要尝试?
二维码

扫码加我 拉你入群

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

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

2012-4-22 23:03:50
epoh 发表于 2012-4-22 08:29
这篇文献是LSTVAR(4)
双变量: 产出缺口,通货膨胀      
转换变量(St),选滞后二期通胀
呵呵,我就是看到文章的后面用的G=0和G=1两种环境下分别做的脉冲,那相当于是分段了?分段的话其实就可以写成线性的模型了,那是不是就可以直接像传统的那样做脉冲了呢?
二维码

扫码加我 拉你入群

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

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

2012-4-23 10:25:43
雁茗轩 发表于 2012-4-22 23:03
呵呵,我就是看到文章的后面用的G=0和G=1两种环境下分别做的脉冲,那相当于是分段了?分段的话其实就可以 ...
我说过这篇文献是LSTVAR(4), 2 regimes(通货紧缩,通货膨胀),
                双变量: 产出缺口,通货膨胀      
                转换变量(St),选滞后二期通胀
由于转换变量(St),是选滞后二期通胀(非external transition variable)
所以可以做GIRF

如果你要GIRF,就需改变你的transition variable
二维码

扫码加我 拉你入群

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

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

2012-4-23 14:41:44
epoh 发表于 2012-4-23 10:25
我说过这篇文献是LSTVAR(4), 2 regimes(通货紧缩,通货膨胀),
                双变量: 产出缺口,通货膨胀 ...
赞epoh
二维码

扫码加我 拉你入群

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

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

2012-4-23 23:48:58
epoh 发表于 2012-4-23 10:25
我说过这篇文献是LSTVAR(4), 2 regimes(通货紧缩,通货膨胀),
                双变量: 产出缺口,通货膨胀 ...
恩,是,看来我的还是不好做
二维码

扫码加我 拉你入群

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

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

2012-4-25 11:30:48
雁茗轩 发表于 2012-4-23 23:48
恩,是,看来我的还是不好做
问题可能出在
你的thVar的取值范围
我再确认一下再贴上来
二维码

扫码加我 拉你入群

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

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

2012-4-25 11:50:39
觉得很好。谢谢了
二维码

扫码加我 拉你入群

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

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

2012-4-25 14:35:57
雁茗轩 发表于 2012-4-23 23:48
恩,是,看来我的还是不好做
yan=read.table("d.txt",header=TRUE)
y=yan$CPI
x=yan$do1           # [0.00 1.00 3.81 1.65....10.84 1.40 1.71 1.93]
x1=yan$do1[2:120]   # [1.00 3.81 1.65....10.84 1.40 ]
x2=yan$do1[3:121]   # [3.81 1.65....10.84 1.40 1.71 ]
x3=yan$do1[4:122]   # [1.65....10.84 1.40 1.71 1.93]
mod1=lstar(y,m=3,d=1,thVar=x1,control=list(maxit=3000))
mod1
LSTAR model
Coefficients:
Low regime:
     const1      phi1.1      phi1.2      phi1.3
0.27269434  0.72507134 -0.03027709  0.14946705

High regime:
     const2      phi2.1      phi2.2      phi2.3
-0.11229556 -0.48641219 -0.08182624 -0.21365960

Smoothing parameter: gamma = 100

Threshold
Variable: external
Value: -0.04647
x1.jpeg
###########
mod2=lstar(y,m=3,d=1,thVar=x2,control=list(maxit=3000))
mod2
LSTAR model
Coefficients:
Low regime:
    const1     phi1.1     phi1.2     phi1.3
0.1246322  0.3213363 -0.2095328  0.1370120

High regime:
    const2     phi2.1     phi2.2     phi2.3
0.2047646  0.2967180  0.3461295 -0.3922440

Smoothing parameter: gamma = 100

Threshold
Variable: external
Value: 2.993
x2.jpeg
##########
mod3=lstar(y,m=3,d=1,thVar=x3,control=list(maxit=3000))
mod3
LSTAR model
Coefficients:
Low regime:
     const1      phi1.1      phi1.2      phi1.3
0.20220828  0.56478303 -0.08529851  0.02761796

High regime:
     const2      phi2.1      phi2.2      phi2.3
-0.14253275 -0.70277458  0.07884266  0.03247553

Smoothing parameter: gamma = 110.5

Threshold
Variable: external
Value: 4.786
x3.jpeg
二维码

扫码加我 拉你入群

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

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

2012-4-25 21:04:17
很好很热心和强大
二维码

扫码加我 拉你入群

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

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

2012-4-27 07:58:08
二维码

扫码加我 拉你入群

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

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

2012-4-27 08:00:03
二维码

扫码加我 拉你入群

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

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

2012-5-3 00:38:59
epoh 发表于 2012-4-25 14:35
yan=read.table("d.txt",header=TRUE)
y=yan$CPI
x=yan$do1           # [0.00 1.00 3.81 1.65....10 ...
但是epoh老师,这个是选取了do1的不同取值分别作为转换变量来你和模型,接着做了脉冲分析吗?那岂不是x1做转换变量得到的脉冲就可以作为我那个模型的了吗?老麻烦您,真不好意思!谢谢您了……
二维码

扫码加我 拉你入群

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

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

2012-5-3 10:49:45
雁茗轩 发表于 2012-5-3 00:38
但是epoh老师,这个是选取了do1的不同取值分别作为转换变量来你和模型,接着做了脉冲分析吗?那岂不是x1做 ...
哈哈,休假完看到了.
你说对了,x1,x2,x3做转换变量可以得到脉冲(你的模型)
不过一阶差分,smooth效果不佳,图形如下:
do1_tf.bmp

我改以thVar=log(oilprice)[3:121]
oil=cpo$Oil
oilln=log(oil)
thVar=oilln[3:121]

smooth效果较好,fitted.values也较好,你参考看看
请注意短信息
logoil_girf.bmp
logoil_tf.bmp
二维码

扫码加我 拉你入群

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

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

2012-5-4 08:39:20
epoh老师,您这是我的数据,谢谢您了。
附件列表

od.xls

大小:88.5 KB

 马上下载

二维码

扫码加我 拉你入群

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

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

2012-5-4 12:04:08
老师请问这个结果是什么意思

Using default threshold variable: thDelay=0
Testing linearity...   p-Value =  0.2999844
The series is linear. Using linear model instead.

Non linear autoregressive model

AR model
Coefficients:
     const      phi.1      phi.2
0.70289952 0.08038689 0.22178193
二维码

扫码加我 拉你入群

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

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

2012-5-5 23:15:32
怀想的季节 发表于 2012-5-4 12:04
老师请问这个结果是什么意思

Using default threshold variable: thDelay=0
说明该模型不具有非线性性,建立线性模型就可以了!
二维码

扫码加我 拉你入群

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

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

2012-5-22 23:00:56
epoh 发表于 2012-5-3 10:49
哈哈,休假完看到了.
你说对了,x1,x2,x3做转换变量可以得到脉冲(你的模型)
不过一阶差分,smooth效果不 ...
老师,我想请教您一下,咱们估计出来的的LSTAR模型的前后两个线性部分的形式是一样,即是同阶的,那可不可以说前面的线性部分是3阶的,后面乘以转换函数的线性部分不是三阶呢?比如是滞后一期?这种可以估计出来吗?谢谢老师……
二维码

扫码加我 拉你入群

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

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

2012-5-23 09:32:55
雁茗轩 发表于 2012-5-22 23:00
老师,我想请教您一下,咱们估计出来的的LSTAR模型的前后两个线性部分的形式是一样,即是同阶的,那可不可 ...
lstar(x, m, d=1, steps=d, mL, mH, mTh,control=list(),....)
  mL : autoregressive order for 'low' regime (dafult: m). Must be <=m
  mH : autoregressive order for 'high' regime (default: m). Must be <=m

可以利用mL,mH 来调整
mod1=lstar(y,m=3,d=1,thVar=oilln[3:121],mL=3,mH=1,control=list(maxit=3000))
> mod1

Non linear autoregressive model

LSTAR model
Coefficients:
Low regime:
      const1       phi1.1       phi1.2       phi1.3
0.190190160  0.721915078 -0.084351533 -0.002905174

High regime:
     const2      phi2.1
0.02355544 -0.50692537

Smoothing parameter: gamma = 100

Threshold
Variable: external
Value: 3.507
二维码

扫码加我 拉你入群

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

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

2012-5-23 12:15:26
epoh 发表于 2012-5-23 09:32
lstar(x, m, d=1, steps=d, mL, mH, mTh,control=list(),....)
  mL : autoregressive order for 'low ...
太感谢您了,好的,我试试!
二维码

扫码加我 拉你入群

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

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

栏目导航
热门文章
推荐文章

说点什么

分享

扫码加好友,拉您进群
各岗位、行业、专业交流群