Mario, here is one way:
(* generation of similar data: *)
dd = .1;
data = Table[{2 + x^2 + RandomReal[{-dd, dd}], x + RandomReal[{-dd, dd}]}, {x, -2, 2, .1}];
(* working on x- and y-component seperately: *)
{datax, datay} = Transpose[data];
(* we then have: *)
ListPlot[{datax, datay}]

(* this then can be fitted in a regular way: *)
modelx = x2 t^2 + x1 t + x0;
modely = y1 t + y0;
paramsx = FindFit[datax, modelx, {x0, x1, x2}, t];
paramsy = FindFit[datay, modely, {y0, y1}, t];
func[t_] = ({modelx, modely} /. paramsx) /. paramsy;
ParametricPlot[func[t], {t, 1, Length[data]}, Epilog -> {Red, Point[data]}]

Does that help? Regards -- Henrik