R/normW.R

In hughesevoanth/moosefun: Functions by Moose

#' A function to provide a Shapiro-Wilk w-statistic
#'
#' This function performs a Shapiro-Wilk test of normality and returns just the W-statistic
#' @param y a vector of numerical values
#' @keywords shapiro, normality
#' @return a vector of lenght 1 providing the W-statistic
#' @export
#' @examples
#' normW()
normW = function( y ){
  ## verify that x is a vector.
  if(is.vector(y) == FALSE){ stop("id_outliers() parameter y is expected to be a vector of data.") }

  ## verify that there is some variability
  if (length(unique(y)) == 1)
    stop("trait is monomorphic")
  if (length(unique(y)) == 2)
    stop("trait is binary")

  ## if there are more than 5000 observations
  ## sample down to n=5000
  if(length(y) <= 5000){
    W = shapiro.test( y )$stat; names(W) = "W"
  } else {
    W = shapiro.test( sample(y,5000) )$stat; names(W) = "W"
  }

  return(W)
}