So, I’m curious… who are you?  What are you using?  Who knows, if a lot of people are using something, I may just improve on it.

If you know me, say hi.  Heck, even if you don’t know me, say hi also.

grid <- expand.grid(x=lon,y=lat)

grid2 <- data.frame(x=grid$x,y=grid$y,eof=eof)

I want you to add one thing to part 1: After you calculate and plot the L function, calculate and plot the K function.

Data for the second part are in this zip file.

The first one, indicatorSim.R, looks at how you could use multiple indicator kriging to simulate an entire dataset. In the last lab, we used single indicator kriging to simulate the probability that a pixel was above 100 ppb with the jura dataset. Here, we will look at a range of values, from low to high.

The second one, lab6.GLS.R, applies Generalized Least Squares to the meuse data. The lead concentrations are first regressed on distance to the river using OLS. The residuals are then extracted, and a semivariogram is fit to them, and GLS is calculated.