R: Shouldn't generic methods work internally within a package without it being attached?

103 views Asked by At

I am writing a package that uses classes and functions from the spatial sp library. sp exports methods for rbind (am I correct in calling rbind a generic?).

For instance, the following code creates two SpatialPoints objects and then uses rbind.SpatialPoints to join them together:

> crdsA <- matrix(c(1,2,3,4), ncol = 2)
> crdsB <- matrix(c(7,8), ncol = 2)
> 
> sptsA <- sp::SpatialPoints(crdsA)
> sptsB <- sp::SpatialPoints(crdsB) 
> 
> sp::rbind.SpatialPoints(sptsA, sptsB)
SpatialPoints:
     coords.x1 coords.x2
[1,]         1         3
[2,]         2         4
[3,]         7         8
Coordinate Reference System (CRS) arguments: NA 

However, if I then convert the SpatialPoints to SpatialPointsDataFrame (a higher level object class within the sp library), and then use rbind.SpatialPointsDataFrame, I get an error:

> sptsdfA <- sp::SpatialPointsDataFrame(sptsA, data.frame(IDs = c(1,2)))
> sptsdfB <- sp::SpatialPointsDataFrame(sptsB, data.frame(IDs = 3))
> 
> sp::rbind.SpatialPointsDataFrame(sptsdfA, sptsdfB)
Error in rbind2(..1, r) : 
  no method for coercing this S4 class to a vector

A look at the rbind.SpatialPointsDataFrame source code reveals that it calls rbind for SpatialPoints:

rbind.SpatialPointsDataFrame <- function(...) {
    dots = list(...)
    names(dots) <- NULL # bugfix Clement Calenge 100417
    sp = do.call(rbind, lapply(dots, function(x) as(x, "SpatialPoints")))
    df = do.call(rbind, lapply(dots, function(x) x@data))
    SpatialPointsDataFrame(sp, df, coords.nrs = dots[[1]]@coords.nrs)
}

So this seems to be the problem, but I do not understand why. If I attach the sp library, then none of these problems occur, but I thought that since rbind was being called internally within rbind.SpatialPointsDataFrame, then the rest of the library did not to be attached.

Within the context of the package I am creating, even if I include import(sp) and importFrom(sp,rbind.SpatialPoints) in the NAMESPACE, the code above does not work.

I guess there is clearly something I am not understanding with regards to loading, attaching and importing packages. Could anyone explain why sp::rbind.SpatialPointsDataFrame does not work without the rest of the library being attached, and how I can get it to work within my package?

Thanks a lot!

1

There are 1 answers

1
Edzer Pebesma On

rbind is not a usual generic: since its signature has ... as its first (and only) argument, it cannot dispatch on the first argument. It has been problementic in all cases to program methods for it. Have you looked at maptools::spRbind?