如何做出LSTAR和estar的散点图？

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

如题，请教大虾们，不胜感激！

扫码加我 拉你入群

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

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

分享0 收藏4 回帖

关键词：ESTAR LSTAR STAR TAR 散点图 不胜感激 散点图 如何

回复你的短信息.
#####
library(tsDyn)
mod <- lstar(log10(lynx), m=2, mTh=c(0,1), control=list(maxit=3000))

phi1=mod$model.specific$coefficients[1:3]
phi2=mod$model.specific$coefficients[4:6]
gamma=mod$model.specific$coefficients[7]
th=mod$model.specific$coefficients[8]

z=mod$model.specific$thVar
G <- function(y, g, th)  plogis(y, th, 1/g)
tf= G(z,gamma,th)
plot(z,tf)
   

epoh 发表于 2012-4-30 15:47
回复你的短信息.
#####
library(tsDyn)

EPOH：你好！
    为你的学识及为人所感动，通前期的有关你的贴子而学习，让我深入了解到了用R到处理STAR模型的技术；通过学习后翻阅资料，让我对SATA有更深的了解。
    既然STAR是一个较前沿的问题，那么为什么不发动起来以你为首开设一个专题呢：省得学子们到处找资料和学习，那怕是我们学子们交些费用都行！
    好从STAR基础开始，还如一元的STAR或多元的SATR的R办法都更好，省的你一个一个解答。你也很辛苦的。数据吗，就用网友提供的数据，这样更有针对性和实效性。

心若灿烂 发表于 2012-5-1 18:54
EPOH：你好！
    为你的学识及为人所感动，通前期的有关你的贴子而学习，让我深入了解到了用R到处理STA ...

哈哈,老兄您言重了.
个人所学,也仅止于能跟各位切磋而已.

epoh 发表于 2012-5-1 20:06
哈哈,老兄您言重了.
个人所学,也仅止于能跟各位切磋而已.

epoh兄就是这么低调又富有热情
赞一个~

epoh 发表于 2012-5-1 20:06
哈哈,老兄您言重了.
个人所学,也仅止于能跟各位切磋而已.

考虑一下，开一个专题如何！！

epoh 发表于 2012-4-30 15:47
回复你的短信息.
#####
library(tsDyn)

就用你上面的例子得到了的结果，也存在不显著的系数：那么要不要进行处理呢：
> data=(lynx)
> mod=lstar(log10(lynx), m=2, mTh=c(0,1), control=list(maxit=3000))
> phi1=mod$model.specific$coefficients[1:3]
> phi2=mod$model.specific$coefficients[4:6]
> gamma=mod$model.specific$coefficients[7]
> th=mod$model.specific$coefficients[8]
> z=mod$model.specific$thVar
> G=function(y,g,th) plogis(y, th, 1/g)
> tf=G(z,gamma,th)
> plot(z,tf)
> mod
> summary(mod)

Non linear autoregressive model

LSTAR model
Coefficients:
Low regime:
    const1     phi1.1     phi1.2
0.4891014  1.2465399 -0.3664328
High regime:
    const2     phi2.1     phi2.2
-1.0240758  0.4232669 -0.2546088
Smoothing parameter: gamma = 11.15
Threshold
Variable: Z(t) = + (0) X(t) + (1) X(t-1)
Value: 3.339
Residuals:
      Min        1Q    Median        3Q       Max
-0.594820 -0.107360  0.014309  0.111098  0.510342
Fit:
residuals variance = 0.03805,  AIC = -357, MAPE = 5.58%

Coefficient(s):
        Estimate  Std. Error  t value  Pr(>|z|)   
const1  0.489101    0.204914   2.3869 0.0169929 *  
phi1.1  1.246540    0.067871  18.3663 < 2.2e-16 ***
phi1.2 -0.366433    0.104301  -3.5132 0.0004427 ***
const2 -1.024076    2.430066  -0.4214 0.6734492   
phi2.1  0.423267    0.172146   2.4588 0.0139415 *  
phi2.2 -0.254609    0.585416  -0.4349 0.6636207   
gamma  11.153834   10.004728   1.1149 0.2649120   
th      3.339199    0.092748  36.0030 < 2.2e-16 ***
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
Non-linearity test of full-order LSTAR model against full-order AR model
F = 12.446 ; p-value = 1.3815e-05
Threshold
Variable: Z(t) = + (0) X(t) + (1) X(t-1)
       在结果中，有不显著的系数：phi2.2 （其P值   0.6636207）；
       如何去掉这个不显著的phi2.2呢？一直被这个问题纠结！
谢谢！

epoh 发表于 2012-5-1 20:06
哈哈,老兄您言重了.
个人所学,也仅止于能跟各位切磋而已.

epoh老师，请问在R里面怎么判断使用ESTAR还是LSTAR呢？急，望回复！

您好！我根据您的程序跑了一下。
phi1=modmodel.specificcoefficients[1:3]
phi2=modmodel.specificcoefficients[4:6]
gamma=modmodel.specificcoefficients[7]
th=modmodel.specificcoefficients[8]
这部分程序中，会报告错误。
错误: 找不到对象'modmodel.specificcoefficients'
请问是怎么回事呀？

phi1=mod$coefficients[1:3]