1
$\begingroup$

Suppose you have N correlated data points and M parameters to fit a model. Then, you can construct the $\chi^2$ function as

$\chi^2 = (\vec{y}_{data} - \vec{y}_{theo}) \cdot C^{-1} \cdot (\vec{y}_{data} - \vec{y}_{theo})$ ,

where $\vec{y}_{data}$ is the vector with the data points, $\vec{y}_{theo}$ is the theoretical vector with the M parameters that follows the model under test and $C^{-1}$ is the covariance matrix of data points. Thus, to best fit the model, you need to minimize $\chi^2$ with respect to the M parameters, i.e.,

$\frac{\partial \chi^2}{\partial a_i} = 0, \quad i=0,1,...,M-1$.

What is the best function on Mathematica to do this job? In addition, a correlation matrix must be a sub-product of this procedure. NonlinearModelFit didn't work in this situation, and I'm not finding any useful build-in function.

Thanks in advance, Gabriel.

$\endgroup$

closed as off-topic by MarcoB, m_goldberg, Henrik Schumacher, Coolwater, halirutan May 29 '18 at 21:14

This question appears to be off-topic. The users who voted to close gave this specific reason:

  • "This question cannot be answered without additional information. Questions on problems in code must describe the specific problem and include valid code to reproduce it. Any data used for programming examples should be embedded in the question or code to generate the (fake) data must be included." – MarcoB, m_goldberg, Henrik Schumacher, Coolwater, halirutan
If this question can be reworded to fit the rules in the help center, please edit the question.

  • 1
    $\begingroup$ I am quite confident that it isLinearModelFit or NonlinearModelFit what should be employed here. What happend when you tried NonlinearModelFit? What was the code? Why do you think it did not work? How the the model look? Many. many questions that are open... $\endgroup$ – Henrik Schumacher Mar 10 '18 at 16:55
  • $\begingroup$ You should use NMinimize $\endgroup$ – OkkesDulgerci Mar 10 '18 at 19:42
  • 1
    $\begingroup$ If you have real data and the model with parameters, please share to get quick response. $\endgroup$ – OkkesDulgerci Mar 10 '18 at 19:43
  • $\begingroup$ You didn't even mention if the model is linear or not. $\endgroup$ – J. M. will be back soon Mar 10 '18 at 22:14
  • $\begingroup$ The model isn't linear. I found an alternative to this type of analysis: one can minimize the function with Nminimize and then define a $\Delta \chi^2$ function such that $\Delta \chi^2 = \chi^2 - \chi_{min}^2$, where $\chi_{min}^2$ is the value of the function in the minimum. Then, you can consult tables of confidence intervals to plot the regions of interest using the ConditionalExpression in the Plot. Anyway, thanks all for the time and the help. Hope this helps, too. Best regards, Gabriel $\endgroup$ – Gabriel Das Neves Mar 11 '18 at 14:43