R/Gmedian.R

In Gmedian: Geometric Median, k-Medians Clustering and Robust Median PCA

Gmedian <- function (X, init = NULL, gamma = 2, alpha = 0.75, nstart=2, epsilon=1e-08) 
{
  X <- as.matrix(X)
  if (!is.null(init)) X = rbind(init,X) ### initialisation   
  med.X = Gmedianrowvec_rcpp(X, gamma=gamma, alpha=alpha, nstart=nstart, epsilon=epsilon)
  return(med.X)
}


Weiszfeld <- function (X, weights = NULL, epsilon=1e-08, nitermax = 100) 
{
  X <- as.matrix(X)
  if (is.null(weights)) poids <- rep(1,nrow(X)) 
  else poids <- weights
  return(Weiszfeld_rcpp(X, poids, epsilon=epsilon, nitermax=nitermax))
}