--- a
+++ b/tests/testthat/test-fromDataAccrual.R
@@ -0,0 +1,122 @@
+context("test-from-data-accrual")
+
+test_that("Poisson",{
+  expect_error(Generate.PoissonAccrual("sd",5))
+  expect_error(Generate.PoissonAccrual("2015-10-10",0))
+  expect_error(Generate.PoissonAccrual("2015-10-10",-1))
+  ac.p <- Generate.PoissonAccrual("2015-01-01",1)
+  expect_equal("Poisson process",ac.p@model)
+  expect_equal("a poisson process with rate=1 and start date=2015-01-01.",ac.p@text)
+  res <- ac.p@f(5000)
+  expect_equal(5000,length(res))
+  expect_equal("Date",class(res))
+  expect_true(all(res >= as.Date("2015-01-01")))
+  set.seed(10)
+  r <- ac.p@f(10)
+  set.seed(20)
+  r2 <- ac.p@f(10)
+  set.seed(10)
+  r3 <- ac.p@f(10)
+  expect_equal(r,r3)
+  expect_true(any(r!=r2))
+})
+
+
+
+test_that("Deterministic",{
+  expect_error(Generate.Accrual(start.date="2015-01-01",end.date="2015-01-01",k=2))
+  expect_error(Generate.Accrual(start.date="2015-10-01",end.date="2015-01-01",k=2))
+  expect_error(Generate.Accrual(start.date="2015-01-01",end.date="2015-10-01",k=0))
+  expect_error(Generate.Accrual(start.date="2015-01-01",end.date="2015-10-01",k=-2))
+  
+  ac.d <- Generate.Accrual(start.date="2015-01-01",end.date="2015-12-31",k=1,deterministic=TRUE)
+  expect_equal("Power law allocation",ac.d@model)
+  expect_equal("a non-stochastic allocation following G(t)=t^k/B^k, where k=1 and B=364 days is the recruitment period [2015-01-01, 2015-12-31].",
+               ac.d@text)
+  
+  
+  res <- ac.d@f(365)
+  expect_equal(as.character(as.Date("2015-01-01")+0:364), as.character(res))
+  
+  ac.d <- Generate.Accrual(start.date="2015-01-01",end.date="2015-01-30",k=2,deterministic=TRUE)
+  expect_equal(3,sum(ac.d@f(12)<= as.Date("2015-01-16")))
+  expect_equal(6,sum(ac.d@f(12)<= as.Date("2015-01-22")))
+  
+  ac.d <- Generate.Accrual(start.date="2015-01-01",end.date="2015-01-30",k=0.78,deterministic=TRUE)
+  set.seed(10)
+  r <- ac.d@f(25)
+  set.seed(20)
+  r2 <- ac.d@f(25) 
+  expect_equal(r,r2)
+  expect_equal(25,length(r))
+    
+  expect_error(Generate.Accrual(start.date="2015-01-01",end.date="2015-10-01",k=2,rec.start.date="2015-01-02"))
+  ac.d <- Generate.Accrual(start.date="2015-01-01",end.date="2015-12-31",k=1,deterministic=TRUE,rec.start.date="2014-01-01")
+  expect_equal(as.character(as.Date("2015-01-01")+0:364), as.character(res))
+  
+  expect_equal(
+  "a non-stochastic allocation following G(t)=t^k/B^k, where k=1 and B=729 days is the recruitment period [2014-01-01, 2015-12-31]. New subjects will be recruited after 2015-01-01.",
+  ac.d@text)
+  
+  ac.d <- Generate.Accrual(start.date="2015-01-01",end.date="2015-01-11",k=0.5,deterministic=TRUE,rec.start.date="2014-12-31")
+  
+  expect_equal((0.5*(sqrt(11)-1)+1)^2,1+as.numeric(ac.d@f(2)[1]-as.Date("2015-01-01")))
+  
+  
+})
+
+test_that("stochastic_k",{
+  ac.s <- Generate.Accrual(start.date="2014-01-01",end.date="2014-12-31",k=1)
+  
+  out <- "a stochastic allocation following G(t)=t^k/B^k, where k=1 and B=364 days is the recruitment period [2014-01-01, 2014-12-31]."
+  expect_equal(out,ac.s@text)
+  
+  
+  set.seed(10)
+  r <- ac.s@f(10)
+  set.seed(20)
+  r2 <- ac.s@f(10)
+  set.seed(10)
+  r3 <- ac.s@f(10)
+  
+  expect_equal(10,length(r))
+  expect_equal(r,r3)
+  expect_true(any(r!=r2))
+  
+  set.seed(10)
+  ac.s <- Generate.Accrual(start.date="2015-01-01",end.date="2015-12-31",k=2,rec.start.date=as.Date("2014-01-01"))
+  res <- ac.s@f(10000)
+  expect_true(all(res < as.Date("2016-01-01")) && all(res>=as.Date("2015-01-01")) ) 
+  
+})
+
+
+test_that("new_accrual",{
+  
+  NewAccrual <- function(dates){
+  
+  myf <- function(N){
+    sample(x=dates,size=N,replace=TRUE)
+  }
+  
+  new("AccrualGenerator",f=myf,text="new",model="test")
+  
+  }
+  
+  
+  na <- NewAccrual(as.Date(c("2015-01-01","2014-01-01")))
+  d <- as.character(na@f(1000))
+  expect_true(all(d %in% c("2015-01-01","2014-01-01")))
+  
+})
+
+test_that("estimate_k",{
+  expect_error(MLestimateK(-8,c(1,2,3)))
+  expect_error(MLestimateK(5,c(-0.4,1,1)))
+  
+  expect_equal(1/4,MLestimateK(exp(4),c(1,1,1)))
+  a <- 1/(log(40)-(log(5)+log(9)+log(6)+log(10))/4)
+  expect_equal(a,MLestimateK(40,c(5,9,6,10)))
+  
+  
+})
\ No newline at end of file