tests/testthat/test-Dstat.R

In twopexp: The Two Parameter Exponential Distribution

test_that("test-Dstat works", {
  data <- c(12,21,26,27,29,29,48,57,59,70,74,153,326,386,502)
  theta_hat <- min(data)           #estimate theta
  beta_hat <- mean(data)-min(data) #estimate beta
  expect_equal( Dstat(data,theta_hat,beta_hat),Dstat(x=data,theta=theta_hat,beta=beta_hat) )

})