data-raw/R/create_ref.nl5defSGAlengte.R

#' Create table of Dutch birth length references
#' @export
create_ref.nl5defSGAlengte <- function() {
  x <- data.frame(matrix(NA, 36L, 8L))
  x[, 1L] <- rep("pinkeltje", 36L)
  x[, 2L] <- c(rep("NL", 36L))
  x[, 3L] <- c(rep(c("male", "female"), 18L))
  x[, 4L] <- rep(0L, 36L)
  x[, 5L] <- rep(1, 36L)
  x[, 6L] <- c(28.49, 28.60, 30.68, 30.66, 32.80, 32.66,
               34.85, 34.58, 36.82, 36.44, 38.71, 38.22,
               40.52, 39.93, 42.26, 41.57, 43.91, 43.12,
               45.47, 44.59, 46.94, 45.97, 48.32, 47.27,
               49.60, 48.47, 50.79, 49.59, 51.88, 50.62,
               52.87, 51.56, 53.76, 52.42, 53.76, 52.42)
  x[, 7L] <- c(.1261, .0933, .1134, .0955, .1025, .0951,
               .0931, .0924, .0846, .0874, .0769, .0811,
               .0699, .0740, .0638, .0670, .0586, .0605,
               .0541, .0549, .0502, .0504, .0468, .0470,
               .0438, .0447, .0413, .0432, .0395, .0424,
               .0382, .0422, .0377, .0425, .0377, .0425)
  x[, 8L] <- rep(25:42, each = 2L)

  dimnames(x)[[2L]] <- c("pop", "sub", "sex", "x", "L", "M", "S", "ga")
  x
}
growthcharts/growthscreener documentation built on Oct. 25, 2024, 4:07 p.m.