gates_inf <-
function(dx, dy, marks, par=list(a=1, b=4, smark=1)) {
# Gates et al influence functions [(s/b)^a - R^a]^(1/a)
with(as.list(par),
pmax(0, marks[[smark]]^a - (b * sqrt(dx^2 + dy^2))^a)^(1/a)
)
}
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.