diff --git a/R/bearing_to_group_centroid.R b/R/bearing_to_group_centroid.R
new file mode 100644
index 0000000..f9a35b2
--- /dev/null
+++ b/R/bearing_to_group_centroid.R
@@ -0,0 +1,28 @@
+#' Calculate absolute bearing to group centroid
+#'
+#' @param DT expects group_mean columns generated with group_centroid
+#' @param coords character vector of column names for x, y
+bearing_to_group_centroid <- function(DT, coords = NULL) {
+  pre <- 'group_mean_'
+
+  stopifnot(length(coords) == 2)
+
+  xcol <- first(coords)
+  ycol <- last(coords)
+  group_xcol <- paste0(pre, xcol)
+  group_ycol <- paste0(pre, ycol)
+
+  stopifnot(xcol %in% colnames(DT))
+  stopifnot(ycol %in% colnames(DT))
+  stopifnot(group_xcol %in% colnames(DT))
+  stopifnot(group_ycol %in% colnames(DT))
+
+  DT[, bearing_centroid := fifelse(
+    .SD[[xcol]] == .SD[[group_xcol]] &
+      .SD[[ycol]] == .SD[[group_ycol]],
+    NaN,
+    atan2(.SD[[group_ycol]] - .SD[[ycol]],
+          (.SD[[group_xcol]] - .SD[[xcol]]))
+  )]
+  return(DT[])
+}
diff --git a/R/distance_to_group_centroid.R b/R/distance_to_group_centroid.R
new file mode 100644
index 0000000..5bd1039
--- /dev/null
+++ b/R/distance_to_group_centroid.R
@@ -0,0 +1,34 @@
+#' Calculate distance to group centroid
+#'
+#' @param DT expects group_mean columns generated with group_centroid
+#' @param coords character vector of column names for x, y
+distance_to_group_centroid <- function(DT, coords, group = 'group',
+                                     return_rank = FALSE) {
+  pre <- 'group_mean_'
+
+  stopifnot(length(coords) == 2)
+
+  xcol <- first(coords)
+  ycol <- last(coords)
+  group_xcol <- paste0(pre, xcol)
+  group_ycol <- paste0(pre, ycol)
+
+  stopifnot(xcol %in% colnames(DT))
+  stopifnot(ycol %in% colnames(DT))
+  stopifnot(group_xcol %in% colnames(DT))
+  stopifnot(group_ycol %in% colnames(DT))
+  stopifnot(group %in% colnames(DT))
+
+
+  DT[, dist_group_centroid :=
+       sqrt((.SD[[xcol]] - .SD[[group_xcol]])^2 +
+              (.SD[[ycol]] - .SD[[group_ycol]])^2)]
+
+  if (return_rank) {
+    DT[, N_by_group := .N, by = c(group)]
+    DT[, rank_dist_group_centroid :=
+         data.table::frank(dist_group_centroid),
+       by = c(group)]
+  }
+  return(DT[])
+}