I have a data frame containing a time series with two time stamp columns, d$day and d$time, and say, for simplicity, one measured variable d$
d$day
d$time
d$
Did you try this?
val_eod <- aggregate(d$val1, by = list(d$day), FUN = tail, n = 1)