R/mkBmat.R

In smacofx: Flexible Multidimensional Scaling and 'smacof' Extensions

#' Auxfunction1
#' 
#' only used internally
#' @param x matrix
mkBmat <- function (x) {
    d <- rowSums (x)
    x <- -x
    diag (x) <- d
    if(any(is.na(x))) stop("An error from deep within the bowels: NAs in the MY matrix (the majorizing function). This is likely because of too many 0s as dissimilarities or distances after the respective transformation and weighting is applied. You can try to supply other transformation parameters or remove offending observations.")
    if(any(!is.finite(x))) x[!is.finite(x)] <- sign(x[!is.finite(x)])*(max(x[is.finite(x)]))#New as this was a bug
    return (x)
}