r - 检查函数 Gammad 和 Truncate 的返回值（来自包 distr 和 truncdist）

Question

搜索论坛后，我没有找到类似的问题。如果你找到了，请告诉我。我真的很感激。

在 R 中，我需要检查函数 Gammad 和 Truncate（来自 lib distr 和 truncdist）的返回值。

这意味着如果他们无法生成 Gammad 和 Truncate pdf，则可以返回失败值或异常，以便我可以处理它。

   G0 <- Gammad(scale = s, shape = sh) 
   # what if Gammad() fails ? 
   TG <- Truncate(G0, lower = lowerbound, upper = upperbound) 
   # what if Truncate() fails ? 

谢谢 ！

Answer 1

从rgamma帮助页面：“无效的参数将导致返回值 NaN，并带有警告。”

如果这是您所看到的，您可以使用

ow <- options("warn")
options(warn=2)
G0 <- try(Gammad(scale = s, shape = sh), silent=TRUE)
if(inherits(G0, "try-error")) # handle invalid arguments
options(warn=ow)