R/fordiff.R

#' Forward Difference Contrast Maker
#'
#' Creates a contrast matrix for forward difference coding used in regression
#' @param k Number of levels in the factor
#' @export

fordiff <- function(k) {
  cmat = matrix(nrow = k, ncol = k -1)  # Create an empty matrix

  for (i in 1:ncol(cmat)) {
    cmat[1:i,i] = k-i
    cmat[(i+1):nrow(cmat),i] = i * -1
  }

  cmat = cmat / k
  return(cmat)
}
ndiquattro/quatts documentation built on May 23, 2019, 1:28 p.m.