Ошибка при попытке поднять дискретную случайную величину до отрицательной степени

Question

library(discreteRV)

X <- RV(c(0, 9), c(1/2, 1/2))
X ** (-2)

всякий раз, когда я пытаюсь запустить это, я получаю следующую ошибку:

Error in RV(as.numeric(outcomes(X))^Y, probs(X), fractions = attr(X, "fractions"),  : 
  More probabilities/odds than outcomes provided

Пакет от GitHub: discreteRV

Answer 1

Проблема в том, что вы пытаетесь поднять ноль до отрицательной степени, что по сути составляет вычисление 1/0 (в данном случае вы пытаетесь найти (1/0) ^ 2) - неопределенную форму. При замене 0 на 1 код работает нормально

> X <- RV(c(1, 9), c(1/2, 1/2))
> X ** (-2)
Random variable with 2 outcomes

Outcomes 0.0123456790123457                  1
Probs                   1/2                1/2