As a part of my work I need the Maximum Likelihood Estimate that maximizes the log-likelihood of a vector. Can I use the in-built function 'mle' for this purpose?

21 次查看（过去 30 天）
BIPIN SAMUEL
BIPIN SAMUEL 2022-8-25
评论： Torsten 2022-9-6
In my work I have a vector of particular length 'L'. I want to maximize the log-likelihood estimate of the vector using maximum likelihood estimation. I have followed the following procedure:
x=[vector(1:L)];
[phat,pci]=mle(x);
des_signal=normpdf(x,phat(1),phat(2));
figure();
plot(des_signal)
However, I donot know whether the above command maximizes the log-likelihood. Also, the output I got was not the exact output that I needed but is somewhat similar. The distribution is assumed to be Normal distribution. Is there any other way to do the problem?
  1 个评论
Torsten
Torsten 2022-9-2
编辑：Torsten 2022-9-2
However, I donot know whether the above command maximizes the log-likelihood.
phat is such that it maximizes the log-likelihood.
Also, the output I got was not the exact output that I needed but is somewhat similar.
I already asked you to explain what you need.
Is there any other way to do the problem?
Why ? What is insufficient with what you get from "mle" ?

请先登录，再进行评论。

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

回答（1 个）

Shreyas
Shreyas 2022-9-2
编辑：Shreyas 2022-9-2
The ‘mle’ command does not return maximum log-likelihood estimates. Log of the likelihood is used for convenience as it simplifies the calculation for many of the derivations.
Please go through the following links for more information: -
log likelihood can be generated by taking log of the output after applying the ‘normpdf’ function.
Example,
>> des_signal = log(normpdf(x, phat(1), phat(2)));
The following functions can be used to get the sum of negative of log likelihood: -
Example,
>> pd = fitdist(x, >Normal>);
>> nll = negloglik(pd);
Regarding the second part of the question, in the code the function ‘normpdf’ evaluates the specified array of scalar values. For generating the normal curve, the input values must be sorted.
Example,
>> x = sort(x);
>> des_signal = normpdf(x, phat(1), phat(2));
  4 个评论
显示 2更早的评论
Torsten
Torsten 2022-9-6
No. You have the parameters of the distribution (mu and sigma) using the "fitdist" command. Knowing these values, you can plot the distribution. The maximum likelihood value cannot be sufficient to fix the distribution.

请先登录，再进行评论。

类别

在 Help Center 和 File Exchange 中查找有关 Exploration and Visualization 的更多信息

标签

Community Treasure Hunt

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

Start Hunting!

Translated by