[Rd] segfault isoreg with NAs

[Tobias Verbeke]

I was able to reproduce the problem under Windows (R 2.6.1).
When running ir4 <- isoreg(1:10, y4 <- c(5, 9, 1:2, 5:8, NA, 8)),
the following message appears:

"R for Windows GUI front-end has encountered a problem and
needs to close. We are sorry for the inconvenience"

after which R closes.

HTH,
Tobias

 > sessionInfo()
R version 2.6.1 (2007-11-26)
i386-pc-mingw32

locale:
LC_COLLATE=English_United States.1252;LC_CTYPE=English_United 
States.1252;LC_MONETARY=English_United 
States.1252;LC_NUMERIC=C;LC_TIME=English_United States.1252

attached base packages:
[1] stats     graphics  grDevices utils     datasets  methods   base

> As can be seen below, adding a NA to the y values
> in a call to isoreg results in a segfault.
> 
> ir4 <- isoreg(1:10, y4 <- c(5, 9, 1:2, 5:8, NA, 8))
> 
> Adding missing values to the x values, on the contrary,
> gives an error, but maybe the error message could be
> tailored to this particular situation.
> 
> y <- c(5, 9, 1:2, 5:8, 3, 8)
> x <- c(1:9, NA)
> isoreg(x, y)
> ## error message: Error in if (!isOrd) { : missing value where 
> TRUE/FALSE needed
> 
> Please find below a (temporary) patch (against Revision 43692)
> for both the R source and the help file.
> 
> Kind regards,
> Tobias
> 
> ### patch isoreg.R ###
> 
> --- isoreg.R	2007-12-14 19:07:47.000000000 +0100
> +++ isoreg2.R	2007-12-14 19:11:20.000000000 +0100
> @@ -18,6 +18,9 @@
>   ##
>   isoreg <- function(x, y=NULL)
>   {
> +    if (any(is.na(x))) stop("x may not contain NA values")
> +    if (any(is.na(y))) stop("y may not contain NA values")
> +
>       xy <- xy.coords(x,y)
>       x <- xy$x
>       isOrd <- (!is.null(xy$xlab) && xy$xlab == "Index") || !is.unsorted(x)
> 
> ### patch isoreg.Rd ###
> 
> --- isoreg.Rd	2007-12-14 19:08:12.000000000 +0100
> +++ isoreg2.Rd	2007-12-14 19:15:00.000000000 +0100
> @@ -20,6 +20,7 @@
>     \item{x, y}{%in \code{isoreg},
>       coordinate vectors of the regression points.  Alternatively a single
>       plotting structure can be specified: see \code{\link{xy.coords}}.
> +    The coordinate vectors may not contain missing values.
>     }
>   }
>   \details{
> 
> 
> ### sessionInfo() information segfault ###
> 
>> sessionInfo()
> R version 2.6.0 (2007-10-03)
> i486-pc-linux-gnu
> 
> locale:
> LC_CTYPE=en_US.ISO-8859-15;LC_NUMERIC=C;LC_TIME=en_US.ISO-8859-15;LC_COLLATE=en_US.ISO-8859-15;LC_MONETARY=en_US.ISO-8859-15\
> ;LC_MESSAGES=en_US.ISO-8859-15;LC_PAPER=en_US.ISO-8859-15;LC_NAME=C;LC_ADDRESS=C;LC_TELEPHONE=C;LC_MEASUREMENT=en_US.ISO-885\
> 9-15;LC_IDENTIFICATION=C
> 
> attached base packages:
> [1] stats     graphics  grDevices utils     datasets  methods   base
>> ir4 <- isoreg(1:10, y4 <- c(5, 9, 1:2, 5:8, NA, 8))
> 
>   *** caught segfault ***
> address 0x24, cause 'memory not mapped'
> 
> Process R segmentation fault (core dumped) at Fri Dec 14 17:48:22 2007
> 
> ______________________________________________
> R-devel at r-project.org mailing list
> https://stat.ethz.ch/mailman/listinfo/r-devel
> 
>