Description Usage Arguments Value Examples

Estimates regression coefficients and spatial autocorrelation
parameters, given spatial coordinates, a model formula, and a
stratification variable. Arguments are the same here as they are
for `slmfit()`

, with an extra argument for
`stratacol`

, the name of the stratification column. Note that
stratum can either by incorporated as a covariate in
`slmfit()`

, in which case the errors have the same
spatial covariance, or, models with differing spatial covariances
for the errors can be fit to each level of stratum, as is done here
in `stratafit()`

.

1 2 3 4 5 6 7 8 9 10 |

`formula` |
is an |

`data` |
is the data set with the response column, the covariates to be used for the block kriging, and the spatial coordinates for all of the sites. |

`xcoordcol` |
is the name of the column in the data frame with x coordinates or longitudinal coordinates |

`ycoordcol` |
is the name of the column in the data frame with y coordinates or latitudinal coordinates |

`stratacol` |
is the name of the stratification column |

`areacol` |
is the name of the column with the areas of the sites. By default, we assume that all sites have equal area, in which case a vector of 1's is used as the areas. |

`CorModel` |
is the covariance structure. By default, |

`estmethod` |
is either the default |

a list of class `slmfit`

with

the spatial covariance estimates

the regression coefficient estimates

the covariance matrix of the fixed effects

minus two times the log-likelihood of the model

the names of the predictors

the sample size

the name of the covariance model used

a vector of residuals

the design matrix

a vector of the sampled densities

a list containing

formula, the model formula

data, the data set input as the

`data`

argumentxcoordcol, the name of the x-coordinate column

ycoordcol, the name of the y-coordinate column

estmethod, either REML or ML

CorModel, the correlation model used

estimated covariance matrix of all sites

Inverted covariance matrix on the sampled sites

the vector of areas.

1 2 3 4 5 6 | ```
data(exampledataset) ## load a toy data set
exampledataset$strata <- c(rep("A", 19), rep("B", 21))
strataobj <- stratafit(formula = counts ~ pred1 + pred2,
data = exampledataset, stratacol = "strata",
xcoordcol = 'xcoords', ycoordcol = 'ycoords', areacol = 'areavar')
summary(strataobj)
``` |

Embedding an R snippet on your website

Add the following code to your website.

For more information on customizing the embed code, read Embedding Snippets.