"Calculations based on either Euclidean or Manhattan distance require projected data to accurately measure distances." therefore it does not work with longitude latitude (or using a package like [circular (Circular Statistics)][1])

With projected data

    x = c(203595.92857397342, 203711.9297397244, 203992.57772137996, 203756.83341678928, 203610.89646632838, 203393.86202718143, 203524.83108528735, 203696.96184736944, 203363.9262424715, 203887.8024748952, 204018.77153300116, 204131.03072566338, 204426.6465996739)
    y = c(89463.34428931243, 89343.60115047272, 89246.30985016546, 89152.76052294695, 89253.79379634294, 89343.60115047272, 89122.82473823702, 88991.85568013109, 89609.28123977332, 89687.86267463688, 89556.89361653096, 89901.15514069512, 89833.79962509777)
    mean_centerX <- mean(x)
    mean_centery <- mean(y)
    standard_deviationX = sd(x)
    standard_deviationY = sd(y)
    standard_distance <- sqrt(sum(((x -mean_centerX)^2(y-mean_centerY)^2))/(length(x)))
    library(plotrix)
    plot(x,y)
    points(mean_centerX,mean_centery,col="red",pch=16)
    draw.circle(mean_centerX,mean_centery,radius=standard_distance,border="red",lwd=2)

[![enter image description here][2]][2]


  [1]: https://cran.r-project.org/web/packages/circular/index.html
  [2]: https://i.sstatic.net/drFCy.png