[Rd] R crashes (PR#7037)

From: <cornulier_at_cebc.cnrs.fr>
Date: Thu 01 Jul 2004 - 06:06:29 EST


Full_Name: thomas cornulier
Version: 1.9.0 and 1.9.1
OS: Win XP
Submission from: (NULL) (194.254.155.62)

the following function produces R crashes under windows XP platform i386-pc-mingw32

arch     i386           
os       mingw32        
system   i386, mingw32  
status                  
major    1              
minor    9.1            
year     2004           
month    06             
day      21             

language R

start code-----
ts.cor<- function(Variable, Xcoord, Ycoord, Year){   ts<- tapply(Variable, list(Year, paste(Xcoord, Ycoord)), mean)   for(j in 1:500) {
    samp<- sample(seq(ncol(ts)), replace= T)     cor.ts.boot<- cor(ts[, samp], use= "p", method= "spearman")     cor.ts.boot<- cor.ts.boot[lower.tri(cor.ts.boot)]   }
}

V<- rpois(180, 0.05)
X<- rep.int(runif(20), 9)
Y<- rep.int(runif(20), 9)

Yr<- rep(1:9, each= 20)
b<- ts.cor(V, X, Y, Yr)

end code------

use= "c" or use= "a" in cor() cause systematic crash, whereas use= "p" produces crashes or inconsistent error messages:
Error in cor.ts.boot[lower.tri(cor.ts.boot)] :

        object is not subsettable
In addition: There were 50 or more warnings (use warnings() to see the first 50)

Error in as.vector(x, mode) : cannot coerce to vector In addition: There were 22 warnings (use warnings() to see them)

Error: invalid type/length (4/399) in vector allocation In addition: There were 50 or more warnings (use warnings() to see the first 50)

Error in cor.ts.boot[lower.tri(cor.ts.boot)] :

        Unimplemented feature in type2str In addition: There were 50 or more warnings (use warnings() to see the first 50)

Error in cor.ts.boot[lower.tri(cor.ts.boot)] :

        Unimplemented feature in type2str In addition: There were 22 warnings (use warnings() to see them)

Behaviour similar whatever the method specification in cor(), but V has to contain many ties for the errors/crashes being produced. (problems fixed by jittering V)
V<- jitter(rpois(180, 0.05), a=0.001)

Crash occured at line 6 of the function when using debug()

Hope I didn't miss something obvious,
Thomas



R-devel@stat.math.ethz.ch mailing list
https://www.stat.math.ethz.ch/mailman/listinfo/r-devel Received on Thu Jul 01 06:10:03 2004

This archive was generated by hypermail 2.1.8 : Wed 03 Nov 2004 - 22:45:00 EST