Take the 2-minute tour ×  I am very new at using GIS data and
only modestly experienced with R. I've been reading about how to
analyze spatial data using the spatial-analyst.net PDF book, so
I'm not completely lost, but I thought I could describe my
problem and people might suggest ideas.  I have a dataset with
about 2000 measurements at different lat/long coordinates,
although I will probably subdivide this dataset as the data was
collected over 3 years and conditions changed over time. Let's
call the variable being measured "IP."  I want to create a map
of IP in the full area under question using Kriging or some
other interpolation method on the sample data. Then I want to
create a histogram measuring the amount of land in various IP
buckets. I will also need to create a histogram which shows the
number of samples in each bucket (note a sample could have a
higher or lower actual IP than what kriging predicts for its
land).  I follow how to load the data into a
SpatialPointsDataFrame and run a kriging analys