Hi Rubén,

Try this,

my.func <- set.coords.lims #if you want this geoR function latter.
set.coords.lims <-function(coords)apply(coords,2,range) 
plot(z1)
points(z1,pt.divide=c("data.proportional"),lambda=1,col="gray",yl="",xl="")

## to get set.coords.lims back.
set.coords.lims <- my.func
plot(z1)

 

______________________
Gledson Luiz Pichasrki
PET Estatística - UFPR
http://www.leg.ufpr.br/~gledson
http://www.leg.ufpr.br/pet




----- Mensagem original ----
De: Rubén Roa-Ureta <rroa@udec.cl>
Para: r-sig-geo@stat.math.ethz.ch
Enviadas: Domingo, 19 de Outubro de 2008 21:11:48
Assunto: [R-sig-Geo] Problem with points.geodata and plot.geodata

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="")

_______________________________________________
R-sig-Geo mailing list
R-sig-Geo@stat.math.ethz.ch
https://stat.ethz.ch/mailman/listinfo/r-sig-geo


__________________________________________________


	[[alternative HTML version deleted]]


