From btyner at gmail.com  Wed Oct 10 16:35:24 2007
From: btyner at gmail.com (btyner at gmail.com)
Date: Wed, 10 Oct 2007 16:35:24 +0200 (CEST)
Subject: [Rd] corMatrix crashes with corARMA structure (PR#9952)
Message-ID: <20071010143524.E4A3F288008C@mail.pubhealth.ku.dk>

Full_Name: Benjamin Tyner
Version: 2.6.0 RC 2007-10-01 r43043
OS: WinXP
Submission from: (NULL) (171.161.224.10)


platform       i386-pc-mingw32                       
arch           i386                                  
os             mingw32                               
system         i386, mingw32                         
status         RC                                    
major          2                                     
minor          6.0                                   
year           2007                                  
month          10                                    
day            01                                    
svn rev        43043                                 
language       R                                     
version.string R version 2.6.0 RC (2007-10-01 r43043)

I have seen this in other versions/platforms as well. Brian Ripley informs me
the segfault is in corStruct.c

Code to reproduce:

n <- 100

# example from Box and Jenkins p. 83
arcoefs <- c(0.8)
macoefs <- c(-0.6)
p <- length(arcoefs)
q <- length(macoefs)

require(nlme)
tmp <- corARMA(value=c(arcoefs,macoefs), form=~1, p=p, q=q)
Sigma <- corMatrix(tmp, covariate = 1:n) # segfault



[Rd] corMatrix crashes with corARMA structure (PR#9952)
btyner at gmail.com