急：R软件加载程序包tsDyn问题 [推广有奖]

epoh老师，我想问一下，前面您在7楼发上了个estar的程序包，那个好像是S-PLUS软件中用的，这个是在R中用的吗？谢谢您了！

你用91楼estar_gscale.R

老师麻烦问一下，这个要怎么安装到R中呢？为什么我加载的过程中失败了？

> utils:::menuInstallLocal()
错误于read.dcf(file.path(pkgname, "DESCRIPTION"), c("Package", "Type")) :
  无法打开链结
此外: 警告信息：
In read.dcf(file.path(pkgname, "DESCRIPTION"), c("Package", "Type")) :
  无法打开压缩文件'estar_gscale/DESCRIPTION'，可能是因为'No such file or directory'
> load("C:\\Users\\nanyan\\Desktop\\estar_gscale.R")
错误: 復原文件幻数出错(文件可能有损坏)-- 没有载入任何数据
此外: 警告信息：
‘estar_gscale.R’文件的魔号是est
   在R版本2号之前的保存功能不再有用
谢谢您了！

不需要载入package
只需要
source("estar_gscale.R")
#########
source("estar_gscale.R")
ndx <- read.table("ndx.rvol.txt");
ndx=ndx[,1]
mod.estar <- estar(log(ndx), m=2,  control=list(maxit=3000))
mod.estar
ESTAR model
Coefficients:
Low regime:
      phi1.0      phi1.1      phi1.2
-0.02512237  0.82408380  0.19486576
High regime:
      phi2.0      phi2.1      phi2.2
-2.14572870 -0.74500270  0.05390954
Smoothing parameter: gamma = 1.013
Threshold
Variable: Z(t) = + (1) X(t) + (0) X(t-1)
Value: -2.775

哦，好的，我试了一下，可以估计！但是我想得到参数的t值和P值的统计量，在LSTAR模型中估计完以后利用“summary(mod)”就可以得到，那为什么这个估计ESTAR之后得到的我用“summary(mod)”得到的结果是这个呢？
> source("estar_gscale.R")
> new=read.table("zero.txt",header=TRUE)
> y=new$dy
> x=new$dx
> mod=estar(y,m=2,thVar=x)
Using maximum autoregressive order for low regime: mL = 2
Using maximum autoregressive order for high regime: mH = 2
Using only first 119 elements of thVar
Performing grid search for starting values...
Starting values fixed: gamma =  0.5 , th =  6.946 ; SSE =  56.08782
Optimization algorithm converged
Optimized values fixed for regime 2  : gamma =  0.6174869 , th =  3.660855
> summary(mod)
               Length Class       Mode   
str              8    nlar.struct list   
coefficients     8    -none-      numeric
fitted.values  119    -none-      numeric
residuals      119    -none-      numeric
k                1    -none-      numeric
model            7    data.frame  list   
model.specific  14    -none-      list   
epoh老师老麻烦您，我实在不好意思，确实是有点菜！谢谢您了！

哈哈,你的水平很高的,因为你说对了
2010.10,当时,package tsDyn,function summary()有小bug
就如同29楼所讨论
#error dims [product 224] do not match the length of object [0]
所以estar_gscale.R,当时我并没有写入function summary()
而是写了小程序来计算系数的t-value
自然你用summary(mod)就跑不出结果了

那老师您方便把那个程序贴给我吗？我做个模型，很想要那个，看看模型能用不？谢谢您了！

那老师您方便把那个程序贴给我吗？我做个模型，很想要那个，看看模型能用不？谢谢您了！

呵呵,我说的小程序指的就是170楼的程序

Transition function要改为exponential transition functions

请注意短信息

我看到了短消息了，可以做，非常感谢您！