Computation of derivatives of the rotation number for parametric families of circle diffeomorphisms