tests/testthat/test-cb_2014_nga.R

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

test_that("test cb_2014_nga vs NGA-West2 spreadsheet", {
  expect_equal(cb_2014_nga(M = 5, T = 1000, Rrup = 90, Rjb = 85, Rx = 85, W = 10, Zbot = 15, dip = 80,
                           lambda = 75, Fhw = 1, Vs30 = 350, region = 1, Zhyp = 8, Z25 = 999)$med,
               c(0.004095141, 0.004062823, 0.004186984, 0.004789998, 0.006236259, 0.007934571,
                 0.010375185, 0.011039463, 0.011491364, 0.011473954, 0.009312538, 0.007578534,
                 0.004783762, 0.003035084, 0.001418271, 0.000756553, 0.000317573, 0.000184919,
                 0.000107988, 3.61424E-05, 1.80603E-05, 0.004054394, 0.256388118))

  expect_equal(round(cb_2014_nga(M = 5, T = 1000, Rrup = 90, Rjb = 85, Rx = 85, W = 10, Zbot = 15, dip = 80,
                                 lambda = 75, Fhw = 1, Vs30 = 350, region = 1, Zhyp = 8, Z25 = 999)$sigma,
                     digits = 4),
               c(0.7112, 0.7181, 0.7370, 0.7820, 0.8005, 0.7841, 0.7574, 0.7362, 0.7212, 0.7207, 0.7115,
                 0.7127, 0.7255, 0.7303, 0.7265, 0.7169, 0.7180, 0.7146, 0.7108, 0.6945, 0.6681, 0.7117,
                 0.6512))

})