Dear @ LINALGBS,
unfortunately it is very difficult to understand where the problem is.
you can you maybe try this,
region_sf <- sf::st_transform(x = NIC, crs = 4326)
But indeed we need to have a reprex in order to understand how to help you.
Have a look at FAQ: How to do a minimal reproducible example ( reprex ) for beginners
Best