Re: [R-sig-Geo] bad position of png & klm file in google earth

2008-10-24 Thread Marta Rufino
Hi, Thank you all very much :-) I finally managed to put it working. It was a problem in the code, and not the coordinates. >> A second doubt I have is if we can set apriori, in R the markers projected >> in google earth... how do I do if instead of the pins, I want a ball? Do I >> have to chan

Re: [R-sig-Geo] bad position of png & klm file in google earth

2008-10-24 Thread Marta Rufino
; >> write(" ", filename, append = TRUE) >> write(paste("\t\t\t", getwd(), "/legend.png", sep = ""), >> filename, append = TRUE) >> write(" ", filename, append = TRUE) >> write(" > yunits=\"fraction\

[R-sig-Geo] bad position of png & klm file in google earth

2008-10-24 Thread Marta Rufino
Dear list, I would like to start by acknowledging the wonderful work of this support list :-) I have a question, which should be prety dummy,... but I was looking around in the net, in reports, etc. and could not sort it out. I am trying to export a image plot to google earth. I manage to do

Re: [R-sig-Geo] cross-validation

2008-02-08 Thread Marta Rufino
ternatives? Best wishes and thank you very much once again, Marta Edzer Pebesma wrote: > Marta Rufino wrote: >> Great! > Hi Marta, > > I now have in the example section of gstat.cv help: > > # multivariable; thanks to M. Rufino: > meuse.g <- gstat(id = "zn"

Re: [R-sig-Geo] cross-validation

2008-02-07 Thread Marta Rufino
k if you have duplicate observations, duplicate >> observations lead to a singular matrix. Use the function zerodist() >> to check where the observations are and remove.duplicates() to remove >> them. >> >> cheers, >> Paul >> >> Marta Rufino schreef:

Re: [R-sig-Geo] cross-validation

2008-02-06 Thread Marta Rufino
Hello, yes, I know it is suppose to do it, but I could not find how, because it gives me an error... for example: require(gstat); require(lattice) data(meuse) coordinates(meuse) = ~x + y data(meuse.grid) gridded(meuse.grid) = ~x + y meuse.g <- gstat(id = "zn", formula = log(zinc) ~ 1, data = me

[R-sig-Geo] cross-validation

2008-02-06 Thread Marta Rufino
Dear list members, Is it possible to do cross-validation on multivariate kriging? i.e. apply krige.cv to multivariate kriging in gstat? thank you very much, Best wishes, Marta ___ R-sig-Geo mailing list R-sig-Geo@stat.math.ethz.ch https://stat.ethz.ch/

[R-sig-Geo] review on autocorrelation methods

2008-01-30 Thread Marta Rufino
Hello, There is recent article about 'methods to account for spatial auto-correlation in the analysis of species distribution data", which I think it might clarify on some of the points being adressed and it is interesting to bring up to this discussion list. The authors compare several spatia

[R-sig-Geo] azymuth error, polygons and Lisbon projection

2007-07-26 Thread Marta Rufino
Hello, I have four doubts: 1) why gzAzymuth is giving me Northern direction in this case, for example? require(sp); require(rgdal) kk=data.frame(x=c(575552.8,576295.9), y=c(4103900,4103611)) coordinates(kk)=~x+y plot(kk, typ="o", col=4) gzAzimuth(coordinates(kk)[1,], coordinates(kk)[2,]) 2) h

[R-sig-Geo] graphical doubts: lines colour, order objects and print in only one graph

2007-05-15 Thread Marta Rufino
Dear list members, I have three doubts: :-) 1. How to colour the lines from a shape file (which has been inserted into a spplot, through sp.layout): #For example nc <- readShapeLines("costa.shp", proj4string=CRS("+proj=longlat +datum=WGS84")) spplot(grid[1], sp.layout=list(nc))