Hi @Jason1,
Take a look at this:
library("tidyverse")
my_data <- tibble(bmi = runif(n = 100,
min = 20,
max = 40))
my_data_aug <- my_data %>%
mutate(bmigp = case_when(bmi < 25 ~ 1,
25 <= bmi & bmi < 30 ~ 2,
30 <= bmi ~ 3))
Hope it helps