An alternative approach to get standardization parameters would be to calculate them (again): tmp <- monthly_retail_tbl %>% group_by(Industry) %>% arrange(Month) %>% mutate(Turnover = log1p(x = Turnover)) %>% group_map(~ c(mean = mean(.x$Turnover, na.rm = TRUE ), sd = sd(.x$Turnover, na.rm = TRUE ))) %>% bind_rows() std_mean <- tmp$mean std_sd <- tmp$sd rm( 'tmp' )