请问R能做ESTAR模型吗?

epoh 发表于 2013-1-6 19:06
package tsDyn的架构是class nlar, subclass lstar
所以如果你是要在此架构下修改为estar,
自然需要再 ...

那这些缺少的函数是不是不需要改代码，只要source一下就行了？

bertf 发表于 2013-1-7 02:50
那这些缺少的函数是不是不需要改代码，只要source一下就行了？

嗯,缺少的函数不需要改代码.
package tsDyn的架构既是class nlar
所以配套函数大都放在 "nlar-methods.R"
你只要source("estar.R")
      source("nlar-methods.R")
此外就剩几个函数了,
你可以跟estar.R,也可以跟nlar-methods.R,放在同一文件

PS:不需要加载package tsDyn

epoh 发表于 2013-1-7 14:00
嗯,缺少的函数不需要改代码.
package tsDyn的架构既是class nlar
所以配套函数大都放在 "nlar-methods. ...

高人，请问你对stata熟悉吗？我看到stata里面有个nl命令，专门回归非线性模型的，如果我把estar模型的方程直接写上去能不能用这个nl命令在stata上直接得到结果？

bertf 发表于 2013-1-8 19:15
高人，请问你对stata熟悉吗？我看到stata里面有个nl命令，专门回归非线性模型的，如果我把estar模型的方程 ...

抱歉,stata我不熟.

epoh 发表于 2013-1-9 16:24
抱歉,stata我不熟.

高人，不好意思又要麻烦你了，还是这个模型的问题。我按照你说的方法复制了R文件，然后运行显示少了一个extend函数，然后我查到raster包里有一个extend函数，然后用相同的方法把extend.R复制下来，但是source
之后显示：
> source("extend.R")
in method for ‘extend’ with signature ‘x="Extent"’: no definition for class “Extent”
in method for ‘extend’ with signature ‘x="Raster"’: no definition for class “Raster”
运行lstar.R(我暂时没改函数名)之后显示
> mod<-lstar(Frn,m=1)
错误于function (classes, fdef, mtable)  :
  unable to find an inherited method for function ‘extend’ for signature ‘"ts"’
请问这个怎么解决？

bertf 发表于 2013-1-15 16:06
高人，不好意思又要麻烦你了，还是这个模型的问题。我按照你说的方法复制了R文件，然后运行显示少了一个e ...

function extend()
function extend.list()
function sigmoid()
在misc.R

epoh 发表于 2013-1-15 16:18
function extend()
function extend.list()
function sigmoid()

回复的真快，谢谢！那buildConstants这个函数呢？我怎么能找到哪个函数在哪个R文件里呢？

bertf 发表于 2013-1-15 16:34
回复的真快，谢谢！那buildConstants这个函数呢？我怎么能找到哪个函数在哪个R文件里呢？

miscSETAR.R

epoh 发表于 2013-1-6 19:06
package tsDyn的架构是class nlar, subclass lstar
所以如果你是要在此架构下修改为estar,
自然需要再 ...

epoh老师，您好！
为什么我运行您的代码会出现以下错误？
> source("estar_2013.R")
错误于file(filename, "r", encoding = encoding) : 无法打开链结
此外: 警告信息：
In file(filename, "r", encoding = encoding) :
  无法打开文件'estar_2013.R': No such file or directory
>  svpdx=read.table("svpdx.dat")
错误于file(file, "rt") : 无法打开链结
此外: 警告信息：
In file(file, "rt") : 无法打开文件'svpdx.dat': No such file or directory
>  mod.estar <- estar(svpdx[,1], m=2, d=1, control=list(maxit=3000))
错误: 没有"estar"这个函数
>  summary(mod.estar)
错误于summary(mod.estar) :
  在为'summary'函数选择方法时评估'object'参数出了错: 错误: 找不到对象'mod.estar'

>

本文来自: 人大经济论坛 S-Plus&R专版 版，详细出处参考： https://bbs.pinggu.org/forum.php? ... 1&from^^uid=11232

zhangtao 发表于 2013-1-15 18:12
epoh老师，您好！
为什么我运行您的代码会出现以下错误？
> source("estar_2013.R")

哈哈!zhangtao兄,
代码我还没上传.
estar_2013.rar(need pw)
estar_2013.rar (18.32 KB)

2013-1-15 20:53:31 上传

###########
source("estar_2013.R")
source("nlar-methods.R")
svpdx=read.table("svpdx.dat")
mod.estar <- estar(svpdx[,1], m=2, d=1, control=list(maxit=3000))
summary(mod.estar)