We use ReIm@p
to get the real pairs then use FindFit
.
solsol1 = FindFit[ReIm@p, a x + b, {a, b}, x];
eqnR[x_]eqnR1[x_] := a x + b /. sol;sol1;
eqnR[x]eqnR1[x]
Show[ListPlot[ReIm@p, Joined -> False],
Plot[eqnR[x]Plot[eqnR1[x], {x, 0, 1.4}, PlotStyle -> Red]]