R/misc.R

Defines functions pooledSD scaleLambdasBySide calcSlopeLambdas

## Daniel Gerlanc and Kris Kirby (2010-2012)
## Miscellaneous helper functions

calcSlopeLambdas <- function(slope.levels) {
  ## Converts unit measures into contrast for calculating slope
  ##
  ## Args:
  ##   slope.levels: a named numeric vector of slope levels
  ##
  ## Returns:
  ##  a named vector of scaled lambdas

  if (!is.numeric(slope.levels))
    stop("'slope.levels' must be numeric.")

  if (is.null(names(slope.levels)))
    stop("'slope.levels' must have names.")
  
  avg = mean(slope.levels)
  lmbds = (slope.levels - avg) / sum((slope.levels - avg)^2)
}

scaleLambdasBySide <- function(lambdas) {
  ## Scales negative and positive lambdas to sum to +/- 1, respectively
  ##
  ## Args:
  ##   lambdas: A numeric vector of lambdas
  ##
  ## Returns:
  ##  a scaled numeric vector of lambdas

  signs = sign(lambdas)
  sum.by.side = c(tapply(lambdas, signs, sum))
  neg.idx = signs == -1 
  pos.idx = signs == +1
  
  lambdas[neg.idx]  = lambdas[neg.idx] / abs(sum.by.side["-1"])
  lambdas[pos.idx]  = lambdas[pos.idx] / sum.by.side["1"]
  return(lambdas)
}

pooledSD <- function(vals, grps, pop.sd=FALSE) {
  ## Calculate the pooled standard deviation of a vector.
  ## Args:
  ##  @param x a numeric vector
  ##  @param grps a numeric, character, or factor vector specifying groups
  ##  @param pop.sd calculate the biased population standard deviation instead
  ##         of the unbiased sample standard deviation
  ##
  ##  @return
  ##    An unbiased point estimate of the pooled standard deviation

  grps.ls = split(vals, grps, drop=TRUE)
  ss      = sum(sapply(grps.ls, function(x) sum((x - mean(x))^2)))
  denom   = if (pop.sd) length(vals) else length(vals) - length(grps.ls)
  res     = sqrt(ss / denom)
  return(res)
}

Try the bootES package in your browser

Any scripts or data that you put into this service are public.

bootES documentation built on April 11, 2023, 5:57 p.m.