The objective seems to be a straightforward nonlinear regression. I don't understand why you'd want to construct a chisquared statistic to do so when using NonlinearModelFit
will give you want you want directly:
(* Put in an id for each dataset and join all of the data into a single dataset *)
data = Join[Insert[#, 1, 1] & /@ data1, Insert[#, 2, 1] & /@ data2, Insert[#, 3, 1] & /@ data3];
(* Run the nonlinear regression *)
nlm = NonlinearModelFit[data, Boole[id == 1] e1[x] + Boole[id == 2] e2[x] + Boole[id == 3] e3[x], {a, v}, {id, x}];
nlm["BestFitParameters"]
(* {a -> 0.018304629031942354, v -> 0.3633047083635527} *)
This assumes that all of the regressions have the same error variance.