procrustes.R

# rotate Y to align it to X
mypro <- function(X, Y) {
  XY <- crossprod(X, Y)
  s <- svd(XY)
  rotation <- s$v %*% t(s$u)
  Yrot <- Y %*% rotation
  diff <- X - Yrot
  sum_squared_error <- sum(diff^2)
  list(
    s = s,
    rotation = rotation,
    Yrot = Yrot,
    diff = diff,
    sum_squared_error = sum_squared_error
  )
}