sdfBeta.R

# See: http://stats.stackexchange.com/questions/12148/looking-for-a-name-for-a-mean-influencing-statistic
# Somewhat personalized
 
sdfBeta <- function(numer, denom = 1){
  if(identical(1, denom)){denom <- rep(1, length(numer))}
  Bj <- sum(numer, na.rm = T) / sum(denom, na.rm = T)
  Bjni <- (sum(numer, na.rm = T) - numer) / (sum(denom, na.rm = T) - denom)
  StdError <- sd(Bjni, na.rm = T)
  Value <- (Bj - Bjni) / StdError
  names(Value) <- names(numer)
  return(Value)
  }