update icREML spats example

AparicioJohan · May 3, 2024 · 54ce7eb · 54ce7eb
1 parent b1cd250
commit 54ce7eb
Show file tree

Hide file tree

Showing 10 changed files with 190 additions and 57 deletions.
diff --git a/R/SpATS.R b/R/SpATS.R
@@ -427,6 +427,7 @@ coef_SpATS <- function(model) {
 #' @examples
 #' \donttest{
 #' library(SpATS)
+#' library(agriutilities)
 #' data(wheatdata)
 #' wheatdata$R <- as.factor(wheatdata$row)
 #' wheatdata$C <- as.factor(wheatdata$col)

diff --git a/R/icREML.R b/R/icREML.R
@@ -119,6 +119,46 @@ ic_reml_asr <- function(fm, scale = 1) {
 #' }
 #' @export
 #'
+#' @examples
+#' \donttest{
+#' library(SpATS)
+#' library(agriutilities)
+#' data(wheatdata)
+#' wheatdata$R <- as.factor(wheatdata$row)
+#' wheatdata$C <- as.factor(wheatdata$col)
+#'
+#' m1 <- SpATS(
+#'   response = "yield",
+#'   spatial = ~ PSANOVA(col, row, nseg = c(10, 20), nest.div = 2),
+#'   genotype = "geno",
+#'   genotype.as.random = TRUE,
+#'   fixed = ~ colcode + rowcode,
+#'   random = ~ R + C,
+#'   data = wheatdata,
+#'   control = list(tolerance = 1e-03, monitoring = 0)
+#' )
+#'
+#' m2 <- SpATS(
+#'   response = "yield",
+#'   spatial = ~ PSANOVA(col, row, nseg = c(10, 20), nest.div = 2),
+#'   genotype = "geno",
+#'   genotype.as.random = TRUE,
+#'   fixed = ~colcode,
+#'   random = ~ R + C,
+#'   data = wheatdata,
+#'   control = list(tolerance = 1e-03, monitoring = 0)
+#' )
+#'
+#' rbind.data.frame(
+#'   ic_reml_spt(m1, label = "colcode_rowcode"),
+#'   ic_reml_spt(m2, label = "colcode_no_rowcode")
+#' )
+#'
+#' rbind.data.frame(
+#'   h_cullis_spt(m1),
+#'   h_cullis_spt(m2)
+#' )
+#' }
 #' @author Johan Aparicio
 #' @references
 #' Verbyla, A. P. (2019). A note on model selection using information

diff --git a/docs/articles/how-to-start.html b/docs/articles/how-to-start.html