Palmer penguins dataset (bill length, bill depth, flipper length; NA rows dropped). The model is fit by gradient ascent on the log-likelihood. Drag the 3D view to rotate.