[R] setar function error message
Pascal Oettli
kridox at ymail.com
Mon Jul 23 11:28:04 CEST 2012
Hello,
It works for me (with a warning message), by adding this line before the
setar procedure:
> a <- as.ts(a)
As I don't know the time step, thus it is only a pseudo time-series.
Regards,
Pascal
Le 23/07/2012 18:18, Ario Ario a écrit :
> Hi,
> I apologize for not attaching a complete dataset. The dataset contains only
> 1000 observations. Please find it in the attachment. Thanks
>
> (Ario)
>
>
>
> On Mon, Jul 23, 2012 at 3:46 PM, Pascal Oettli <kridox at ymail.com> wrote:
>
>> Hello,
>>
>> Your problem is nto reproducible without the complete "a" dataset.
>>
>> Regards,
>> Pascal
>>
>>
>> Le 23/07/2012 09:49, Ario Ario a écrit :
>>
>>> Hi all,
>>>
>>> I have problem to estimate a SETAR model. I always get an error message.
>>> Here is the code:
>>>
>>> ## there are 4175 observation in the series (a).
>>>
>>> a[1:10,1] [1] 1.496498 1.496602 1.496636 1.496515 1.496515 1.496463
>>>> 1.496429 1.496549 1.496480
>>>>
>>> [10] 1.496498
>>>
>>> library("tsDyn")
>>>>
>>>
>>> selectSETAR(a, m=2)Using maximum autoregressive order for low regime: mL
>>>> = 2
>>>>
>>> Using maximum autoregressive order for high regime: mH = 2
>>> Searching on 1552 possible threshold values within regimes with
>>> sufficient ( 15% ) number of observations
>>> Searching on 6208 combinations of thresholds ( 1552 ), thDelay ( 1
>>> ), mL ( 2 ) and MM ( 2 )
>>> Results of the grid search for 1 threshold
>>> thDelay mL mH th pooled-AIC
>>> 1 0 2 2 1.674402 -43590.12
>>> 2 0 2 2 1.674218 -43588.81
>>> 3 0 2 2 1.673758 -43588.24
>>> 4 0 2 2 1.674011 -43588.19
>>> 5 0 2 2 1.674480 -43587.85
>>> 6 0 2 2 1.673988 -43587.62
>>> 7 0 2 2 1.673666 -43587.56
>>> 8 0 2 2 1.674471 -43587.16
>>> 9 0 2 2 1.673620 -43586.51
>>> 10 0 2 2 1.673574 -43585.83
>>>
>>> mod.star <- setar(a, m=2, steps=5, thDelay=0, th=1.674402)Error in
>>>> rep(NA, length(x) - length(reg)) : invalid 'times' argument
>>>>
>>>
>>>
>>> Would someone tell me the meaning of this error message? and how to solve
>>> the problem?
>>>
>>>
>>> Kind Regards,
>>> Ario
>>>
>>> [[alternative HTML version deleted]]
>>>
>>> ______________________________**________________
>>> R-help at r-project.org mailing list
>>> https://stat.ethz.ch/mailman/**listinfo/r-help<https://stat.ethz.ch/mailman/listinfo/r-help>
>>> PLEASE do read the posting guide http://www.R-project.org/**
>>> posting-guide.html <http://www.R-project.org/posting-guide.html>
>>> and provide commented, minimal, self-contained, reproducible code.
>>>
>>>
>>
>
More information about the R-help
mailing list