tests/testthat/test_vgSA.R

In ppernot/rgumlib: Uncertainty Propagation (GUM & GUM-Supp1) and Sensitivity Analysis

fExpr = expression(x1+x2)
x.mu  = c(1,1)          
x.u   = c(0.1,1.0)      
x.pdf = c('norm','norm')
x.df  = c(Inf, Inf)
names(x.u) = names(x.mu) = names(x.pdf) = names(x.df) = c('x1','x2')

test_that('vgSA results fit Combination of Variances for linear model',
          {
            nrun = 1000
            Sref = gumCV(fExpr,x.mu,x.u,silent=TRUE)
            S    = gumS1(fExpr,x.mu,x.u,x.pdf,x.df,nrun=nrun,silent=TRUE)
            SA   = vgSA(fExpr,x.mu,x.u,S$X,S$Y,silent=TRUE)
            expect_equal(SA$vg[['x1']],Sref$anova[['x1']], 
                         tolerance = 3/sqrt(nrun))
          })