tests/testthat/test_stats.R In zFactor: Calculate the Compressibility Factor 'z' for Hydrocarbon Gases

```library(testthat)

context("z.stats")

test_that("z.stats matches dimension and mean of MAPE. default interval = coarse", {
sh_corrs <- (zFactor:::z_correlations\$short)
sapply(sh_corrs, function(corr)
expect_equal(dim(z.stats(corr)), c(112,11)))
s.mean <- sapply(sh_corrs, function(corr) mean(z.stats(corr)\$MAPE))
expected <- c(13.8862628,  0.8186578,  0.7697211,  0.8123786,  4.1530402,  0.1635157,  2.0987294)
names(expected) <- sh_corrs
expect_equal(s.mean, expected)
})

test_that("z.stats matches dimension and mean of MAPE, interval = fine", {
sh_corrs <- (zFactor:::z_correlations\$short)
sapply(sh_corrs, function(corr)
expect_equal(dim(z.stats(corr)), c(112,11)))
s.mean <- sapply(sh_corrs, function(corr) mean(z.stats(corr, interval = "fine")\$MAPE))
expected <- c(15.8031751, 0.7581761, 0.7500547, 0.7958098, 4.1555843, 0.1543264, 1.8543540)
names(expected) <- sh_corrs
expect_equal(s.mean, expected)
})

test_that("z.stats matches dimension and mean of MAPE, range = hp", {
sh_corrs <- (zFactor:::z_correlations\$short)
sapply(sh_corrs, function(corr)
expect_equal(dim(z.stats(corr)), c(112,11)))
s.mean <- sapply(sh_corrs, function(corr) mean(z.stats(corr, pprRange = "hp")\$MAPE))
expected <- c(15.8345501,  1.0276664,  0.9324033,  0.9698105,  4.2518234,  0.1668054,  2.5638032 )
names(expected) <- sh_corrs
expect_equal(s.mean, expected)
})

context("z_correlations")

test_that("z_correlations matches the short names of the correlations", {
sh_corrs <- (zFactor:::z_correlations\$long)
for (corr in sh_corrs) {
expect_true(corr %in% c("Beggs-Brill", "Hall-Yarborough", "Dranchuk-AbuKassem",
"Dranchuk-Purvis-Robinson", "Shell", "Ann10", "Papp"))
}
})

test_that("z_correlations matches the short names of the correlations", {
sh_corrs <- (zFactor:::z_correlations\$short)
for (corr in sh_corrs) {
expect_true(corr %in% c("BB", "HY", "DAK", "DPR", "SH", "N10", "PP"))
}
})

context("z.plot.range")

test_that("z.plot.range work for Hall-Yarborough", {
expect_silent(z.plot.range("HY"))
expect_silent(z.plot.range("HY", pprRange = "hp"))
expect_silent(z.plot.range("HY", pprRange = "lp"))

})

test_that("z.plot.range work for Hall-Yarborough", {
expected <- "You have to provide a z-factor correlation:  BB HY DAK DPR SH N10 PP"
expect_error(z.plot.range(), expected)

})
```

Try the zFactor package in your browser

Any scripts or data that you put into this service are public.

zFactor documentation built on Aug. 1, 2019, 5:04 p.m.