Recreating autoplot.cca() from fortify.cca() data frame about ggvegan HOT 2 CLOSED

library('vegan')
library('ggvegan')
data(varespec, varechem)

ord <- cca(varespec ~ Al + P*(K + Baresoil), data=varechem)

autoplot(ord)

ford <- fortify(ord, axes = 1:2)  # fortify the ordination
take <- c('CCA1', 'CCA2')  # which columns contain the scores we want
arrows <- subset(ford, Score == 'biplot')  # take only biplot arrow scores
## multiplier for arrows to scale them to the plot range
mul <- ggvegan:::arrowMul(arrows[, take],
                          subset(ford, select = take, Score == 'sites'))
arrows[, take] <- arrows[, take] * mul  # scale biplot arrows

ggplot() +
  geom_point(data = subset(ford, Score == 'sites'),
             mapping = aes(x = CCA1, y = CCA2)) + 
  geom_segment(data = arrows,
             mapping = aes(x = 0, y = 0, xend = CCA1, yend = CCA2),
             arrow = arrow(length = unit(0.01, "npc"))) +
  geom_text(data = arrows, # crudely push labels away arrow heads
             mapping = aes(label = Label, x = CCA1 * 1.1, y = CCA2 * 1.1)) +
  coord_fixed()

mul <- ggvegan:::arrowMul(arrows[, take],
                          subset(ford, select = take, Score %in% c('sites')))

ggplot() +
  geom_point(data = subset(ford, Score == 'sites'),
             mapping = aes(x = CCA1, y = CCA2)) + 
  geom_segment(data = arrows,
               mapping = aes(x = 0, y = 0, xend = CCA1, yend = CCA2),
               arrow = arrow(length = unit(0.01, "npc"))) +
  geom_text(data = arrows,
            mapping = aes(label = Label, x = CCA1, y = CCA2, hjust = 0.5*(1 - sign(CCA1)), vjust = 0.5*(1-sign(CCA2)))) +
  coord_fixed() +
  scale_x_continuous(expand = c(.1, .1)) +
  scale_y_continuous(expand = c(.1, .1))