Предполагая (надеясь), что у вас есть версия 7 Mma, чтобы вы могли использовать LinearModelFit
, что позволяет извлекать многие свойства встроенной модели, включая "PredictedResponse"
. Например,
ClearAll[x, y, data, pred];
data = Flatten[Table[{x, y, Sin[x + y]}, {x, 5}, {y, 5}], 1];
lm = LinearModelFit[data, {x, y}, {x, y}];
pred = lm["PredictedResponse"];
Transpose@Flatten[{Transpose@data, {pred}}, 1] // MatrixForm
дает

Если вам нужно работать с Fit
, то вы можете сделать следующее:
lmf = Fit[data, {1, x, y}, {x, y}];
fitfunc[x_, y_] := lmf[[1]] + lmf[[2, 1]] x + lmf[[3, 1]] y;
predz = fitfunc @@@ data[[All, 1 ;; 2]];
dataAndPredictions =
Transpose@Flatten[{Transpose@data, {predz}}, 1] // MatrixForm