Edzer Pebesma

As in ```r plot.variogramModel = function(x, ..., cutoff, type = 'l') { if (missing(cutoff)) stop("parameter cutoff needs to be specified") plot(variogramLine(x, cutoff, ...), type = type) } ``` ?

I'll look into; it should be similar to plot.gstatVariogram (using lattice, setting the x axis minimum to 0); you're right about the axis labels.

Yes, if you specify the direction vector: ```r > plot(model2, cutoff = 1600, dir = c(.5 *sqrt(2),.5* sqrt(2),0)) > plot(model2, cutoff = 1600, dir = c(.5 *sqrt(2), -.5* sqrt(2),0)) ```...

It doesn't sound like trouble, but it is too easy to forget. On https://cran.r-project.org/src/contrib/Archive/gstat/ you can see the release dates; I do my best to keep important changes or bug...

I've looked at it, and managed to get the commands run without error up to the last command, but then get ``` > var

Try not naming the arguments, as in ``` krige(lst~longtitude+latitude+slope+evelation, sa, grid.grid, model=fvlst) ```

It needs to be in `grid.grid` as well.

`slope` is on the RHS of the formula, meaning it is a predictor; for prediction, you need to have all predictors at prediction locations.

@BenGraeler @henningte here you go!

This requires spacetime from https://github.com/edzer/spacetime , currently, to run.