Here is the entire code. Trying to fit the data into a piecewise exponential fucntion and testing the maximum likelihood. Where is the mistake. Could not figure it out. Help please..

3 次查看（过去 30 天）
Marwan Abukhaled
Marwan Abukhaled 2018-1-1
评论： Marwan Abukhaled 2018-1-2
The entire code is attached. Please check it out. I spent hours trying to figure the source of the mistake but could not get around it.

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

采纳的回答

John BG
John BG 2018-1-2
Hi Marwan
this is John BG <mailto:jgb2012@sky.com jgb2012@sky.com>
when correcting the following line
[mldln,conf1] = mle(xdata,'pdf',@pdfdln,'start',[mu,sig],'upperbound',[5,5], 'lowerbound',[0.1,0.1],'options',options);
with
[mldln,conf1] = mle(xdata,'pdf',@dln,'start',[mu,sig],'upperbound',[5,5], 'lowerbound',[0.1,0.1],'options',options);
then the function produces a histogram, no errors
.
Does it fix it?
if so would you please be so kind to consider marking my answer as Accepted Answer?
To any other reader, if you find this answer useful please consider clicking on the thumbs-up vote link
thanks in advance for time and attention
John BG
<mailto:jgb2012@sky.com jgb2012@sky.com>
  2 个评论
Walter Roberson
Walter Roberson 2018-1-2
Good catch about pdfdln .
It turns out, though, that another fix is required. In the function dbn, right before the for loop you need to add
z = zeros(rs);
This is needed so that the output is the same orientation as the input; without this change, you get a crash in fminsearch that is difficult to track down.
Marwan Abukhaled
Marwan Abukhaled 2018-1-2
Thank you Walter. This was smart note that I immediately accept it. Together with the above suggestion by *John BG,* the problem is solved.
Thank you for your time and cooperation.

请先登录，再进行评论。

更多回答（1 个）

Walter Roberson
Walter Roberson 2018-1-1
You have
[mldln,conf1] = mle(xdata,'pdf',@pdfdln,'start',[mu,sig],'upperbound',[5,5], 'lowerbound',[0.1,0.1],'options',options);
but pdfdln is not a function you define, and is not a Mathworks supplied function or any function I can find anywhere on the Internet.

类别

在 Help Center 和 File Exchange 中查找有关 Descriptive Statistics 的更多信息

标签

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by