To plot the bei dataset in Google Earth use:

library(spatstat)
library(maptools)
library(rgdal)

data(bei)

lon0lat0 <- data.frame(lat=9.15125, lon=-79.8553)
coordinates(lon0lat0) <- ~lon+lat
proj4string(lon0lat0) <- CRS("+proj=longlat +datum=WGS84")
x0y0 <- data.frame(spTransform(lon0lat0, CRS("+init=epsg:32617")))

bei.pnt <- data.frame(x=x0y0[1,1]+bei$x, y=x0y0[1,2]+bei$y, no=rep(1, length(bei$x)))
coordinates(bei.pnt) <- ~x+y
proj4string(bei.pnt) <- CRS("+init=epsg:32617")
bei.ll <- spTransform(bei.pnt, CRS("+proj=longlat +datum=WGS84"))
writeOGR(bei.ll, "bei.kml", "bei", "KML")