Take the product of two probability density function histograms generated from data

Question

David Gillcrist 2023-9-21

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/2024127-take-the-product-of-two-probability-density-function-histograms-generated-from-data

编辑： William Rose 2023-9-21

I have two different sets of data $D_1$ and $D_2$ both of which have data between two values $a$ and $b$. I know I can get the histograms of these two data sets displayed as probability density functions, but what I'd like to do is take the product of these two probability density functions. What would be the most efficient way of achieving this?

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

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

Answer 1

William Rose 2023-9-21

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/2024127-take-the-product-of-two-probability-density-function-histograms-generated-from-data#answer_1315347

编辑：William Rose 2023-9-21

在 MATLAB Online 中打开

@David Gillcrist,

Since you did not provide data, let's create some.

x1=rand(10000,1)+rand(10000,1);   % pdf(x1) = triangle from 0 to 2 with peak at 1
x2=0.8+0.4*randn(10000,1);       % normal with mu=0.8, sigma=0.4.
x2=x2(x2>=0 & x2<=2);           % discard x2 values <0 and >2

Compute the histograms.

bw=0.1;

h1=histogram(x1,'Normalization','pdf','BinWidth',bw);

figure

h2=histogram(x2,'Normalization','pdf','BinWidth',bw);

Multiply the histograms.

joint=h1.Values.*h2.Values; % bin-by-bin product

Plot result

bar(0.05:bw:1.95,joint,1)

How is that? Note that joint is not a PDF since the area under the curve is not 1. If you want joint to be a pdf, multiply it by the appropriate normmalization factor.