Nanotime formatting with underscore separator for fractional seconds

Question

Working with nanoseconds from epoch is a breeze with nanotime, e.g.:

> x  <- 1639142792837258898
> xn <- nanotime::nanotime(x)
> xn
[1