fitting a model to correlated data

2 次查看（过去 30 天）

Alex 2014-1-21

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/113147-fitting-a-model-to-correlated-data

评论： dpb 2014-1-22

Trying to fit a model (bloch equations from MRI) to experimental data (12 parameters, 12 equation). Monte Carlo Multistart ('lsqcurvefit' in MATLAB) gives expected values for all parameters except for two, which are correlated:

http://imgur.com/1PRayGj (fig. 1) The fig. is a result of the Monte Carlo simulation.

Briefly about the model. It is a system of differential equations (solved numerically by Cramer's rule). The fig. 2 contains 9 diff. eq. My model is a bit more complicated and contain 12 diff. eq.

http://imgur.com/E5LiYdq (the model. fig. 2)

I am passing initial guess, lower and upper bounds for all parameters to the solver. The problem - In the optimized output, 1st parameter always takes the value of its upper limit and 2nd one the value of its lower limit (please, see the fig. 1).

Bad luck, because I am only after those two parameters. Are there any good fitting algorithm for correlated data?

3 个评论
显示 1更早的评论隐藏 1更早的评论

Alex 2014-1-22

Sorry for being unclear. I did Monte Carlo simulation in MATLAB using Multi Start with "lsqcurvefit" (500 repeats). The Fig. 1 is the results of the simulation. From the fig. 1 it is clear that 2 parameters are correlated. Did not get anything better with Monte Carlo simulated annealing. Trying Monte Carlo Random Walk now...

dpb 2014-1-22

I seriously doubt repeating the same calculation w/ various tools is going to change anything...the problem is in the basic formulation that you need to be able to recast or add another independent expression in order to remove the correlation.

Or, potentially, as suggested before, estimate the one in terms of the other and then back-substitute that into the remainder of the system.

请先登录，再进行评论。

请先登录，再回答此问题。