R/spatcen.R

In musto101/wilcox_R: What the package does (short line)

spatcen <-
function(x){
#
# compute spatial median
# x is an n by p matrix
#
if(!is.matrix(x))stop("x must be a matrix")
x<-elimna(x)
START<-apply(x,2,median)
val=optim(START,spat.sub,x=x,method='BFGS')$par             
list(center=val)
}