In the detailed walk-through, section Diverse cluster assignment, there is this line:
diversity_score <- sweep(E / O, 2, theta, '/') %*% as.matrix(Phi)
According to the equation (E/O)^theta, shouldn't '/' be replaced by '^', like:
diversity_score <- sweep(E / O, 2, theta, '^') %*% as.matrix(Phi)
Or am I missing something?