R/partSsd.R

In sib-swiss/dsSwissKnife: DataSHIELD Tools and Utilities - server side

# sum of squared distances (for a saner calculation of global variance):
partSsd <- function (xvect, xpoint ) {
  xect <- xvect[!is.na(xvect)]
  # if no xpoint we use the local mean
  xpoint <- .decode.arg(xpoint)
  if(is.null(xpoint) || length(xpoint) == 0){
    xpoint <- mean(xvect, na.rm = TRUE)
  }


  check <- .dsBase_isValidDSS(xvect)
  if (check) {
    res <- sum((xvect - xpoint)^2)
    len <- length(xvect)
  }
  else {
    res <- NA
    len <- NA
  }
  return(list(ssd = res, len = len))
}