[R-sig-Geo] Problem with points.geodata and plot.geodata

[Rubén Roa-Ureta]

Hi ComRades,

The functions points.geoadata and plot.geodata have an underisable 
behavior with the limits of the X axis of the plot, as shown by the toy 
example below. The range covered by the X axis is too wide. Setting 
specific limits with xlim in, for example points.geoadata, doesn't solve 
the problem.
Does anybody knows of a solution or a way around this feature.
Thanks
Rubén

x <- runif(100,100,200)
y <- runif(length(x),1000,2000)
z <- rnorm(length(x),50,15)
a <- sample(c(1,2),length(x),replace=TRUE)
Q <- data.frame(cbind(x,y,z,a))
z1 <- subset(Q,a==1,select=c(x,y,z))
library(geoR)
#Loading required package: sp
#
#-------------------------------------------------------------
#Analysis of geostatistical data
#For an Introduction to geoR go to http://www.leg.ufpr.br/geoR
#geoR version 1.6-22 (built on 2008-10-08) is now loaded
#-------------------------------------------------------------
z1 <- as.geodata(z1)
plot(z1)
points(z1,pt.divide=c("data.proportional"),lambda=1,col="gray",yl="",xl="")