predict.surface.se: Standard errors of predictions

Description

Evaluates the standard error of predictions on a surface.

Usage

predict.surface.se(out, grid.list=NA, extrap=FALSE, chull.mask,...)

Arguments

out

A fitted model object of a certain class

grid.list

A list with as many components as variables describing the surface. All components should have a single value except the two that give the grid points for evaluation. If the matrix or data frame has column names, these must appear in the grid list. Se

chull.mask

Whether to restrict the fitted surface to be on a convex hull, NA's are assigned to values outside the convex hull. chull.mask should be a sequence of points defining a convex hull. Default is to form the convex hull from the observations if this argument

extrap

Extrapolation beyond the range of the data. If false function will be restricted to the convex hull of the observed data or the convex hull defined from the points from the argument chull.mask.

...

Any additional arguments that will passed to the predict.se function specific to the fit object.

Value

A surface object with components A vector of standard errors for the predicted values.

Details

This function is generic and will call the appropriate function to calculate the standard errors for the object class. It operation is simple a grid is created based on the grid list or from the fit object. The prediction standard error are evaluated on the grid using predict.se. Finally the standard errors are reformed into a surface object suitable for plotting.

Examples

Run this code

fit<-Tps(ozone$x,ozone$y)                # tps fit 
out<- predict.surface.se( fit)

surface( out)

# or ... 

image.plot( out)