aovSar.rcbd: Using a SAR model to handle spatial dependence in a...

In lrcastro/spANOVA: Spatial Analysis of Field Trials Experiments using Geostatistics and Spatial Autoregressive Model

Description

Fit a randomized complete block design when the experimental units have some degree of spatial dependence using a Spatial Lag Model (SAR).

Usage

aovSar.rcbd(resp, treat, block, coord, seq.radius)

Arguments

resp	Numeric or complex vector containing the values of response variable.
treat	Numeric or complex vector containing the treatment applied to each experimental unit.
block	Numeric or complex vector specifying the blocks.
coord	Matrix of point coordinates or a SpatialPoints Object.
seq.radius	Complex vector containing a radii sequence used to set the neighborhood pattern. The default sequence has ten numbers from 0 to half of the maximum distance between the samples.

Details

Three assumptions are made about the error in the analysis of variance (ANOVA):

1. the errors are normally distributed and, on average, zero;

2. the errors all have the same variance (they are homoscedastic), and

3. the errors are unrelated to each other (they are independent across observations).

When these assumptions are not satisfied, data transformations in the response variable are often used to circumvent this problem. For example, in absence of normality, the Box-Cox transformation can be used.

However, in many experiments, especially field trials, there is a type of correlation generated by the sample locations known as spatial correlation, and this condition violates the independence assumption. In this setting, this function provides an alternative for using ANOVA when the errors are spatially correlated, by using a data transformation discussed in Long (1996)

Y_{adj} = Y - (\hat{ρ}WY - \hat{ρ}β_0),

where \hat{ρ} denotes the autoregressive spatial parameter of the SAR model estimated by lagsarlm, β_0 is the overall mean and W is a spatial neighborhood matrix which neighbors are defined as the samples located within a radius, this radius is specified as a sequence in seq.radius. For each radius in seq.radius the model is computed as well its AIC, then the radius chosen is the one that minimizes AIC.

The aim of this transformation is converting autocorrelated observations into non-correlated observations in order to apply the analysis of variance and obtain suitable inferences.

Value

aovSar.rcbd returns an object of class "SARanova". The functions summary and anova are used to obtain and print a summary and analysis of variance table of the results. An object of class "SARanova" is a list containing the following components:

DF	degrees of freedom of rho, treatments, residual and total.
SS	sum of squares of rho, treatments and residual.
Fc	F statistic calculated for treatment.
residuals	residuals of the adjusted model.
p.value	p-value associated to F statistic for treatment.
rho	the autoregressive parameter.
Par	data.frame with the radii tested and its AIC.
y_orig	vector of response.
y_ajus	vector of adjusted response.
treat	vector of treatment applied to each experimental unit.
modelAdj	model of class aov using the adjusted response.
modelstd	data frame containing the ANOVA table using non-adjusted response.
namey	response variable name.
namex	treatment variable name.

References

Long, D.S., 1996. Spatial statistics for analysis of variance of agronomic field trials. In: Arlinghaus, S.L. (Ed.), Practical Handbook of Spatial Statistics. CRC Press, Boca Raton, FL, pp. 251–278.

Rossoni, D. F.; Lima, R. R. . Autoregressive analysis of variance for experiments with spatial dependence between plots: a simulation study. Revista Brasileira de Biometria, 2019

Scolforo, Henrique Ferraço, et al. "Autoregressive spatial analysis and individual tree modeling as strategies for the management of Eremanthus erythropappus." Journal of forestry research 27.3 (2016): 595-603.

Examples

data("rcbd_simulated")

# Fitting the model
model <- aovSar.rcbd(rcbd_simulated$y, rcbd_simulated$trat, rcbd_simulated$block,
                     cbind(rcbd_simulated$coordX, rcbd_simulated$coordY))

# Summary for class SARanova
summary(model)

# Anova for class SARanova
anova(model)

lrcastro/spANOVA documentation built on Aug. 21, 2019, 11:45 a.m.