R: Round datetime to nearest time of day

Question

Say I have a POSIXct vector like

timestamps = seq(as.POSIXct("2021-01-23"), as.POSIXct("2021-01-24"), length.out = 6)

I would