Analysis of correlated circular and extremal data with a flexible cylindrical distribution