r - 尝试将离散随机变量提高到负幂时出错

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

标签： r

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