Learn R Programming

siland (version 1.0)

siland: Estimation of spatial influence of landscape

Description

siland is used to fit spatial influence of landscape.

Usage

siland(loc.model, land=NULL, data, initSIF = NULL, sif = "exponential", 
family = "gaussian", test = FALSE)

Arguments

loc.model

a symbolic description (see lm() or glm()) of the response variable concerning local variables. Random effects are also allowed according to the syntax in package lme4 (see lmer() function in package lme4).

land

a list containing the location of the landscape variables. Each component of the list is a dataframe with two columns "X" and "Y" indicating the locations of pixels for a given landscape variable.

data

a dataframe containing the response variable and the local variables.

initSIF

a vector indicating the starting values for the estimation of the mean distance of the spatial influence functions. The length of initSIF has to be equal to the length of list land. By default, initialisation is equal to 100 for each landscape variable.

sif

the family of the spatial influence function. sif can be "exponential", "gaussian" or "uniform".

family

the distribution of response variable. family can be "gaussian", "poisson" or "binomial" and the associated link function are identity, log and logit respectively.

test

logical value. If TRUE, a ratio likelihood test is performed for each explanatory variable (local and spatial). By default, the value is FALSE since performing tests for all the variables can be heavy computing.

Value

siland returns an object of type list.

coefficients

vector of estimated coefficients

local

an object of class formula that indicates the local model used

landcontri

a dataframe of estimated contributions of each spatial variable (in column) to each observation (in row). The number of columns is equal to the length of list land

loglik

log-likelihood for the estimated parameters

loglik0

log-likelihood for the local model

fitted

fitted values

sif

the family of the spatial influence function

resoptim

an object of class optim or optimize giving informations about the optimization procedure see optim() or optimize() for further details.

AIC

akaike information criterion

AIC0

akaike information criterion for local model (no landscape variable)

nparam

number of parameters

pval0

p.value of the test of global effect of spatial variables. Obtained from the likelihood ratio test between the complete model and the local model.

pval

if test is TRUE, vector of p.values of the test of effect of each landscape variables

family

family distribution for the model

sd.error

standard error for gaussian family, NA in other case

model.Type

type of local model: GLM for generalised model, LMM for linear mixed model or GLMM for generalised linear mixed model

rand.StdDev

standard deviation of random effects for LMM or GLMM

nparam

number of model parameters

References

Chandler R. and Hepinstall-Cymerman J. (2016) Estimating the spatial scales of landscape effects on abundance. Landscape ecology, 31: 1383-1394.

Examples

Run this code
# NOT RUN {
data(dataSiland)
data(landSiland)
resE=siland(y~locvar,land=landSiland,data=dataSiland,sif="exponential",family=gaussian)
resE
resE$AIC
# }
# NOT RUN {
# }

Run the code above in your browser using DataLab