swh:1:snp:33a53053e50f7abe7d281cc0c803be827debf4a3
Tip revision: 298da2447c47ae82148ed053286b951c796a4294 authored by Edzer J. Pebesma on 12 March 2007, 14:19:41 UTC
version 0.9-36
version 0.9-36
Tip revision: 298da24
gstat.formula.R
# $Id: gstat.formula.q,v 1.7 2006-02-10 19:01:07 edzer Exp $
"gstat.formula" <-
function (formula, data)
{
m = model.frame(terms(formula), as(data, "data.frame"))
Y = model.extract(m, response)
if (length(Y) == 0)
stop("no response variable present in formula")
Terms = attr(m, "terms")
X = model.matrix(Terms, m)
has.intercept = attr(Terms, "intercept")
if (gridded(data))
grid = gridparameters(data)
else
grid = numeric(0)
list(y = Y, locations = coordinates(data), X = X, call = call,
has.intercept = has.intercept, grid = as.double(unlist(grid)))
}