Генерация случайной выборки с загрязненным логарифмическим нормальным распределением

Question

Я хотел бы сгенерировать случайную выборку с загрязненным логнормом, распределенным в Matlab.Примером загрязненной логнормально распределенной выборки может быть бимодальное логнормальное распределение, где у вас будет 1 большой узел и 1 маленький узел.

Это то, что я использовал для создания 5 выборок, каждая из которых содержит 10 точек данных.

trueGM=1;
trueGSD=2;
n=10;
m=5;
o = lognrnd(log(trueGM), log(trueGSD),n,m);

Как я могу изменить это для получения загрязненного логнормального распределения?

Answer 1

Это загрязняет логнормальный образец 25% другого логнормального образца.

sampleSize = 10000
pureln = exp(rnorm(n=sampleSize, 3, 4))
contam = runif(n=sampleSize) < .25
other = exp(rnorm(n=sum(contam), 15, 0.01))
mixture = pureln
mixture[contam] = other
hist(log(mixture))