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 <[EMAIL PROTECTED]>
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]]
_______________________________________________
R-sig-Geo mailing list
R-sig-Geo@stat.math.ethz.ch
https://stat.ethz.ch/mailman/listinfo/r-sig-geo
