R/cr_functions.R

In reyzaguirre/st4gi: Statistical tools for genetic improvement

# Create fieldplan array

fp <- function(nr, nc, serpentine) {
  
  plan.id <- t(array(1:(nr*nc), dim = c(nc, nr)))
  
  if (serpentine == 'yes' & nr > 1)
    for (i in seq(2, nr, 2))
      plan.id[i, ] <- sort(plan.id[i, ], decreasing = TRUE)
    
  plan.id 
  
}

# Get most rectangular field shape

gnc <- function(nplot) {

  foo <- function(x, nplot) {
    if (nplot %% x == 0)
      x
  }
  
  f1 <- unlist(sapply(1:nplot, foo, nplot))
  
  d <- data.frame(f1, f2 = nplot / f1)
  d$s <- d$f1 + d$f2
  d <- d[d$s == min(d$s), ]
  d <- d[1, ]
  
  d$f2
  
}