如何從不引發異常的庫函數中捕獲錯誤（`ars`）？

Question

我正在使用ars（自適應拒絕採樣）庫來對凹面密度分佈中的給定變量進行採樣。

事情是，ars功能打印「ifault代碼」，當出現錯誤，但顯然它不會拋出任何異常。

library(ars) 

f<-function(x,mu=0,sigma=1){-1/(2*sigma^2)*(x-mu)^2} 
fprima<-function(x,mu=0,sigma=1){-1/sigma^2*(x-mu)} 
# mysample<-ars(1,f,fprima,mu=2,sigma=3) # it gives no errors 
mysample<-ars(1,f,fprima,mu=2,sigma=3, x=c(2,4)) # bad choice of x 
hist(mysample) 

的問題是，我稱之爲從蟒蛇的ars功能使用rpy2，我需要我的程序知道ars是否已經能夠生成樣本。也就是說，我想ars拋出一個真正的異常或返回一些特殊的價值。

我該如何發現錯誤？

Answer 1

重寫代碼，使結論部分做你想要的。這可能是函數的最後部分的一種可能的變體：

ars2 <- function (n = 1, f, fprima, x = c(-4, 1, 4), ns = 100, m = 3, 
    emax = 64, lb = FALSE, ub = FALSE, xlb = 0, xub = 0, ...) 
{ 
# leave top portion the same 
#...... 
      else { 
       return(list(msg= "Error in sobroutine sample_...", ifault= sample$ifault)) 
      } 
     } 
    } 
    else { 
     return(list(msg="\nError in sobroutine initial_...", ifault= initial$ifault)) 
    } 
    return(mysample) 
} 
# Then finish up by setting the environment for ars2 
environment(ars2) <- environment(ars)