tests/testthat/test-cy_2014_nga.R

In wltcwpf/GMPE: A pacakge of commonly used Ground Motion Prediction Equations (GMPE)

test_that("test cy_2014_nga vs NGA-West2 spreadsheet", {
  expect_equal(cy_2014_nga(M = 5, T = 1000, Rrup = 90, Rjb = 85, Rx = 85, Ztor = 999, dip = 80, lambda = 75,
                           Vs30 = 350, Vs30_code = 0, region = 1, Fhw = 1)$med[-c(6, 10, 12)],
               c(0.282410345, 0.004998711, 0.004998711, 0.005013322, 0.00522232, 0.006095921, 0.007685947,
                 0.009198078, 0.011448513, 0.012634327, 0.012522836, 0.011703983, 0.009560269, 0.007632546,
                 0.004385625, 0.002758741, 0.001336385, 0.000757718, 0.00031836, 0.000155392, 8.3825E-05,
                 2.61643E-05, 1.11451E-05))

  expect_equal(round(cy_2014_nga(M = 5, T = 1000, Rrup = 90, Rjb = 85, Rx = 85, Ztor = 999, dip = 80, lambda = 75,
                                 Vs30 = 350, Vs30_code = 0, region = 1, Fhw = 1)$sigma[-c(6, 10, 12)],
                     digits = 4),
               c(0.7346, 0.7530, 0.7530, 0.7536, 0.7648, 0.7791, 0.7903, 0.7985, 0.8093, 0.8163, 0.8211, 0.8247,
                 0.8301, 0.8331, 0.8208, 0.8020, 0.7696, 0.7601, 0.7550, 0.7513, 0.7482, 0.7426, 0.7387))

})