tests/testthat/test-i_2014_nga.R

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

test_that("test i_2014_nga vs NGA-West2 spreadsheet", {
  expect_equal(i_2014_nga(M = 5.5, Vs30 = 550, T = 1000, Rrup = 90, Fault_Type = 1)$med,
               c(0.008279822, 0.008506422, 0.009335473, 0.010635557, 0.01316374, 0.015581321, 0.018375866,
                 0.018416791, 0.017333836, 0.016135616, 0.013115548, 0.010408766, 0.006875145, 0.004682499,
                 0.002461634, 0.001487336, 0.000705068, 0.000424701, 0.000284073, 0.000138556, 7.75439E-05,
                 0.008279822))

  expect_equal(round(i_2014_nga(M = 5.5, Vs30 = 550, T = 1000, Rrup = 90, Fault_Type = 1)$sigma,
                     digits = 4),
               c(0.7451, 0.7451, 0.7451, 0.7451, 0.7593, 0.7694, 0.7836, 0.7937, 0.8015, 0.8079, 0.8179,
                 0.8257, 0.8399, 0.8500, 0.8642, 0.8743, 0.8885, 0.8885, 0.8885, 0.8885, 0.8885, 0.7451))

})