Новое распределение вероятностей геометрическое распределение со средним 1 / р

Question

Я хотел бы объявить новое распределение вероятностей, мне нужно геометрическое распределение, я знаю, что в R этот тип распределения уже объявлен, но мне нужно геометрическое распределение вида p (1-p) ^ (k- 1), поэтому среднее значение равно 1 / p и k = {1, 2, ...}. И тогда я бы хотел использовать все функции типа rgeom и т. Д.

Большое спасибо за любую помощь.

Answer 1

Геометрическое распределение в R определяется как p(1-p)^(x), т. Е. Количество испытаний, не включающих первый успех.

Чтобы получить распределение для всех испытаний, включая первый успех, вы можете просто настроить формулы соответствующим образом,

my_rgeom <- function(n, prob) rgeom(n, prob) + 1
my_dgeom <- function(x, prob, log = FALSE) {
  p <- dgeom(x, prob) / (1 - prob)
  if (isTRUE(log)) log(p) else p
}