Während Henri in seiner Lösung geschickt Koordinaten als Bezugspunkt setzt, kann man sich auch auf die Bildabmessungen des entstandenen Plots beziehen. Direkt nach Ende des group-Plots:
\path (current bounding box.north west) west)
-- node[yshift=3ex,xshift=-2ex,rotate=90] {$y$}
(current bounding box.south west);
Etwas Shift noch, damit es perfekt an die Mitte des eigentlichen Plots rückt und in Abstand zur Achse (wenn man kein phantom-y setzt).setzt). Dann natürlich das `ylabel` weglassen.