From bolker at ufl.edu  Sat Nov 24 17:45:16 2007
From: bolker at ufl.edu (Ben Bolker)
Date: Sat, 24 Nov 2007 08:45:16 -0800 (PST)
Subject: [Rd] Bug in package stats function ar() (PR#10459)
In-Reply-To: <20071124010521.AEE0B282EFF7@mail.pubhealth.ku.dk>
References: <20071124010521.AEE0B282EFF7@mail.pubhealth.ku.dk>
Message-ID: <13926484.post@talk.nabble.com>




Steven McKinney wrote:
> 
> Full_Name: Steven McKinney
> Version: 2.6.0
> OS: OS X
> Submission from: (NULL) (142.103.207.10)
> 
> 
> 
> Function ar() in package "stats" is showing
> a quirky bug.  Some calls to ar() run to
> completion, others throw an error.
> 
> The bug is reproducible by several people on different
> machines, however, the ar() function itself ends
> up throwing the error sporadically.  Several calls to
> ar() may be necessary to trip the error condition.
> 
> Code to reproduce:
> 
> x<-ts(c(-0.2052083,-0.3764986,-0.3762448,0.3740089,0.2737568,2.8235722,-1.7783313,0.2728676,-0.3273164),start=c(1978,3),frequency=4,end=c(1980,3))
> # ar function
> res.ar<-ar(x,aic=TRUE,demean=F)
> # call "ar" again and ............
> res.ar<-ar(x,aic=TRUE,demean=F)
> 
> 
> Example output:
> (Note that on this attempt the first call to ar()
>  tripped the error.)
> 
>> x<-ts(c(-0.2052083,-0.3764986,-0.3762448,0.3740089,0.2737568,2.8235722,-
> +
> 1.7783313,0.2728676,-0.3273164),start=c(1978,3),frequency=4,end=c(1980,3))
>> 
> 

 Slightly more detail: in the "eureka" function, in stats/src/eureka.f, the
last
element of "vars" in the return list is bogus and varies among calls; if it
comes out to NaN or NA it triggers the bug.

  Haven't spent any more time yet tracking this down, but hope that
will save someone a few minutes.

   Ben Bolker
-- 
View this message in context: http://www.nabble.com/Bug-in-package-stats-function-ar%28%29-%28PR-10459%29-tf4864536.html#a13926484
Sent from the R devel mailing list archive at Nabble.com.



[Rd] Bug in package stats function ar() (PR#10459)
Ben Bolker