Модифицированное экспоненциальное распределение в R

Question

Я хотел бы создать свою собственную функцию плотности вероятности в R для имитации некоторых вещей из бумаги.

Это похоже на экспоненциальное распределение, но я действительно хочу переопределить экспоненциальное распределение в "модифицированное" ...

Есть ли такой способ сделать это? Спасибо.

---

Я хочу смоделировать это:

b (x) = {(µ / p) (e ^ (- (µx-q) / p), x> q (xbar) и 0 в противном случае}

xbar - это х с линией над ним, среднее значение

Answer 1

Я считаю, что вы можете делать то, что вы хотите с ?sample. Используйте вашу известную функцию распределения b(x), чтобы сгенерировать вектор вероятностей, скажем bprob, затем

sample(x,size, replace=TRUE,prob=bprob)

Существует несколько очень интересных методов генерации выборок из произвольных распределений. См., Например, Нормальные случайные числа: использование машинного анализа, чтобы выбрать лучший алгоритм W. H. PAYNE WashlngtonState University (где-то в Интернете) или Численные рецепты на C Искусство научных вычислений, второе издание Уильям Х. Пресс, Сол А. Теукольский и др., Глава 7.3

Answer 2

Знаете ли вы, что экспоненциальные распределения, усеченные снизу, остаются экспоненциальными распределениями?