lx <- split( x, x$hr) #with the next step being lapply( lx, function( x) ...)).

HR <- unique( x$hr)
lx <- select( lx, -hr)
lx <- split( x, HR)