## ----knitr-opts, include = FALSE----------------------------------------------
knitr::opts_chunk$set(
  collapse = TRUE, comment = "#>",
  fig.width = 7.2, fig.height = 4.4, out.width = "100%",
  dpi = 110, fig.align = "center"
)

## ----setup--------------------------------------------------------------------
library(BorderEffect)
data(field_trial)
head(field_trial)

## -----------------------------------------------------------------------------
be <- border_effect(response ~ trt + blk, data = field_trial, coords = field_trial)
be

## -----------------------------------------------------------------------------
bt <- border_test(be, pattern = "outer1", nsim = 300, seed = 42)
bt
plot(bt)

## -----------------------------------------------------------------------------
moran_border(be)

## -----------------------------------------------------------------------------
csv <- system.file("extdata", "field_trial.csv", package = "BorderEffect")
d <- read.csv(csv)
names(d)[names(d) == "block"] <- "blk"   # the CSV names the block column `block`
be2 <- border_effect(response ~ trt + blk, data = d, coords = d)
border_test(be2, pattern = "outer1", nsim = 300, seed = 1)

## ----eval=FALSE---------------------------------------------------------------
# fl <- field_layout(nx = 12, ny = 8, ntrt = 3, nblk = 2,
#                    arrangement = "triangular", width = 550, height = 210,
#                    seed = 3000)
# fl$response <- my_recorded_yields          # length-96 numeric vector, plot order
# be <- border_effect(response ~ trt + blk, data = fl, coords = fl)
# border_test(be, pattern = "outer1")