the average of multiple Probability density functions

Question

Rima Habib 2019-12-5

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/494981-the-average-of-multiple-probability-density-functions

编辑： Ridwan Alam 2019-12-6

Temp.mat

Hallo,

I have a cell array of {1000,1} pds, they are all the same distribution, but they differ slightly in the parameter values

Does any one know how to get the average of these pds? I have got them in iterative process and would like now to know the mean or average of them

I though of taking the mean of the parameter values in the pds in for loop but I am truggling with the right way to do it because I keep getting errors.

it would be helpful if some one gives me a hint :)

I have attached the cell array (Temp) in mat file

Thanks in advanced!

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

Image Analyst 2019-12-5

在 MATLAB Online 中打开

J, it should work regardless of what toolboxes you have

 s=load('temp.mat')
 cellArray = s.Temp
s = 
  struct with fields:
    Temp: {1000×1 cell}
    
cellArray =
  1000×1 cell array
    {1×1 prob.GeneralizedParetoDistribution}
    {1×1 prob.GeneralizedParetoDistribution}
    {1×1 prob.GeneralizedParetoDistribution}
    {1×1 prob.GeneralizedParetoDistribution}
etc.

Jeff Miller 2019-12-5

What do you mean by "the average of the probability distributions"?

For example, suppose one distribution is normal(0,1) and another distribution is normal(10,2). Is the average normal(5,1.5)? Or is the average bimodal with one mode around 0 and another around 10? Or is the average distribution something else?

请先登录，再进行评论。

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

Answer 1

Ridwan Alam 2019-12-5

1
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/494981-the-average-of-multiple-probability-density-functions#answer_405038

The idea of "average" is ill-defined for pdfs.

In your case, all the 1000 distridutions are generalized paretos. Your theta is fixed in all pdfs, so you can certainly take the average of the sigma's and the k's. But what meaning/significance those "averages" carry, you need to come up with that.

If you are looking for the average expected value, you can calculate the mean (expected value) for each pareto and then average those. mean = theta + (sigma/(1-k)). ref: https://en.wikipedia.org/wiki/Generalized_Pareto_distribution

Also, if you are looking for the overall distribution parameters, you can combine the input data [available in Temp{}.InputData.data] and use gppdf().

Hope this helps.

2 个评论
显示无隐藏无

Rima Habib 2019-12-6

@Jiff Miller @Ridwan Alam

I am trying to generate artificial Time series, so after an iterative methode, taking the average of the resulted pdfs is the best way to get the gppdf that suits the data.

so the only idea I could came up with is to take the average of the sigmas and ks but I am also interessted in your other approch which is combinig the inputdata then using gppdf.

how to do that exactly with a cell array full of pdfs. I couldnt figure it out?

Thanks you for your Help!

Ridwan Alam 2019-12-6

编辑：Ridwan Alam 2019-12-6

在 MATLAB Online 中打开

Rima, as I mentioned one approach could be to stack all the inputdata from all the pdfs.

allInputData=[];
for i = 1:length(Temp)
    allInputData = [allInputData; Temp{i}.InputData.data];
end
% in your case, theta is fixed
newTheta = Temp{1}.theta;
paramhat = gpfit(allInputData-newTheta);
newK = paramhat(1);
newSigma = paramhat(2);