similar to: merge maps from shapefile to lattice

Hi everybody: I just generate interpolation maps with differents methods, like IDW and kriging, and now i want to compare the predict values versus real, and i don't who is the commant to do it. I want for example, if I pass from console the coordinates of a place, return me the predict value. Could anybody tell me how should do?. Thanks in advance --

Hi everybody: I just generate different maps of interpolation with different methods but with the same package -(gstat, i used the krige command). And now for evaluate the best method i want to compare the predict values versus experimental, but i don't know how do it. Could anybody tell me how should do it?. Thanks in advance. -- ######################################## Antoni Vi??dez

hi all, I have a plot to which i would like to add text labels. And i cant find a way...here is the code : enaD2<-idw(D2~1, loca=dva, newdata=grd) pts = list("sp.points", dva, pch = 20, cex=1.5, col = "darkred spplot(enaD2, "var1.pred",sp.layout=pts, main = "globina 60 cm", sub="D2",

Hi! Is there a way to manually costumize the color legend in an spplot() - especially where to draw ticks and labels for the ticks? The reason I'm asking: Usually spplot() automatically divides the data into fitting slices and makes a color legend (also automatically). I want to assign the slices myself and have a fixed scale instead of an automatic/dynamic scale. I think what I want

Hi! How does one change the background color of the map-panel in spplot()? Example: library(sp) data(meuse.grid) gridded(meuse.grid) = ~x+y spplot(meuse.grid, "part.a") How would I get another background-color for the map-panel (but not for the whole plot) here? Thank you! Marcel

Hi! Im new to R and I need to interpolate a shapefile using kriging. I've been able to plot/read the shapefile using the package maptools or rgdal. I've searched the internet for sample codes but most of the kriging codes that I've found done in R is done using txtfiles or CSVs.  An example could be of great help. Thanks. [[alternative HTML version deleted]]

Dear all can you help me? i have this problem: i have a dataset in a text file in a matrix of 3 columns: x, y, z where x and y are the coordinates and Z are the mesurements. How can i do a variogram with R?

Hello everyone, When looking for information about clustering of spatial data in R I was directed towards DBSCAN. I've read some docs about it and theb new questions have arisen. DBSCAN requires some parameters, one of them is "distance". As my data are three dimensional, longitude, latitude and temperature, which "distance" should I use? which dimension is related to

Hi, I am very new to R and I am trying to write KML files. I am running R-2.13.1 and I have the following packages: - AKIMA - RGDAL - RSAGA - SP I was trying to run the simple example that I found at the following web page: http://spatial-analyst.net/wiki/index.php?title=Export_maps_to_GE The lines that I can run are: data(meuse.grid)

Hi, I would know how to plot two variograms on a same graph. I can plot one by one but I would draw both on the same one. Is it possible? Do i need any special package? Thanks! Cordialement Damien Landais

Hello! I have a dataset with the dates and positions - Lat and Long columns How do I make a scatter plot? Do I have to cbind to combine the Lat and Long columns? Any suggestions will be much appreciated! Thank you Julia [[alternative HTML version deleted]]

Hi ALL, With the pakage gstat How can one create a grid of points as in the meuse.grid example, from measuring sites as in the Meuse example Thanks in advance! Alain PONSERO Conservateur Réserve Naturelle de la Baie de Saint-Brieuc site de l'Etoile 22120 HILLION tel/fax : 02.96.32.31.40 <mailto:aponsero@cabri22.com> aponsero@cabri22.com

Hello, I'm using the gstat package within R for an automated procedure that uses ordinary kriging. I can see that there is a logical ("singular") atrtibute of some adjusted model semivariograms: .- attr(*, "singular")= logi TRUE I cannot find documentation about the exact meaning and the implications of this attribute, and I dont know anything about the inner calculations

hello, i do the following in order to get an persp-plot x<-c(2,2,2,2,2,2,3,3,3,3) y<-c(41,41,83,83,124,166,208,208,208,208) z<-c(90366,90366,92240,92240,92240,96473,100995,100995,100995,100995) x<-data$x y<-data$y z<-matrix(data$z,length(y),length(x)) persp(x,y,z, col="gray") but i always get the error message increasing 'x' and 'y' values expected,

Hi everyone As I couldn't succeed with manual installation of Rmpi I decided to start again from the beginning. I removed R and MPICH in my Ubuntu Hardy installation. Then, to avoid any dependencies problems I have installed MPICH and R from synaptic, not from sources. But now I can't install Rmpi. An error message appears when trying to install Rmpi, you can find in

I downloaded the iid.test, but I can't run it. I get the following message: Error: could not find function "iid.test" Where am I supposed to save this package in order that it works? Thanks, EZ [[alternative HTML version deleted]]

Hi, I'm using R 2.5.1 and gstat 0.9-39 library and I'm working with the jura data set provided by gstat library. I tried to plot a graph of metal concentrations (let's say Cd) with the command spplot, but I realized that the default lags are equally distributed between Min and Max. I did: library(gstat) data(jura) jura.pred.xy <- jura.pred names(jura.pred.xy)[1:2] <-

Hi R-users I'm a beginner with R and statistics, so I need some help to start my data analysis. I've been reading some docs and tutorials on R and cluster analysis. I've got a large dataset (102000 points) with values of longitude, latitude and temperature and want to see if I can find groups (clusters). Following some tutorials I can look for principal components but get an error

Hello, I am learning how to use functions, but I'm running into a roadblock. I would like my function to do two things: 1) convert an object to a dataframe, 2) and then subset the dataframe. Both of these commands work fine outside the function, but I would like to wrap them in a function so I can apply the code iteratively to many such objects. Here's what I wrote, but it doesn't

Hi everyone I am new to R and have a question about missing values. I am trying to do a cluster analysis of monthly temperatures and my data are 14 columns with spatial coordinates (lat,lon) and 12 monthly values: /lat - lon - temp1 - //temp2 - temp3 - .... - //temp12/ If I omit missing values (my missing values are 99.00) with /mydata <- na.omit(mydata)/ every row with a