Coordinate Center

想到該如何計算多個座標點中心的問題

找了一下好像沒有適合的

就寫了一個小程式來處理這個問題

Y=matrix(c(
5,10,
12,18,
20,13),ncol=2,byrow=TRUE)

fn=function(x,Y) {
x01=x[1]
x02=x[2]
sum(sqrt((Y[,1]-x01)^2+(Y[,2]-x02)^2))
}

optim(par=apply(Y,2,mean), fn=fn,Y=Y,control=list(factr=1e-10,reltol=1e-10))

## $par
## [1] 12.33 15.94
## 
## $value
## [1] 19.74
## 
## $counts
## function gradient 
##       59       NA 
## 
## $convergence
## [1] 0
## 
## $message
## NULL