Estimation of GLMs with unspecified link function: a P-spline approach