Hi. I need to construct income distributions for several countries. The underlying assumption is that income follows a lognormal distribution and i have values for the two parameters. Of course, each distribution comes from a different population size (say China 1371220000 individuals and so on).
I used the code:
x=[0:100:500000]
china=lognpdf(x,'mu','sigma')
but this way matlab returns the values of the PDF.
I tried creating a distribution object like
china=makedist('Lognormal','mu','sigma')
but, in this case, i do not know how to assign population sizes and make operations between objects.
Substantially, how can i create income distributions knowing the two parameters and population sizes? Once created, i need to integrate them (e.g. to see how many people are below a certain threshold) or sum them (e.g. obtain a unique distribution for Chinese + Indians) so take this into account when making suggestions.
