## Test Messages ----
test_that("removal() messages",{
  ## wrong type
  expect_error(removal(c(346,184,49),method="Derek"),
               "should be one of")
  ## wrong type for CS.se
  expect_error(removal(c(346,184,49),method="CarleStrub",CS.se="Derek"),
               "should be one of")
  ## alpha and beta are not positive
  expect_error(removal(c(346,184,49),method="CarleStrub",alpha=1,beta=0),
               "must be positive")
  expect_error(removal(c(346,184,49),method="CarleStrub",alpha=-1,beta=1),
               "must be positive")
  expect_error(removal(c(346,184,49),method="CarleStrub",alpha=-1,beta=0),
               "must be positive")
  ## Catch not in a vector
  expect_error(removal(matrix(c(346,184,49,12),nrow=2)),
               "must be a vector")
  expect_error(removal(data.frame(c(346,184),c(49,12))),
               "must be a vector")
  ## Check characters in catch=
  expect_error(removal(c(346,184,"derek")),
               "must be a vector of numeric values")
  expect_error(removal(c("jason","powell","derek")),
               "must be a vector of numeric values")
  ## only one catch
  expect_error(removal(346),"with one catch value")
  ## Try using Moran or Schnute method with not three catches
  expect_error(removal(c(346,184),method="Moran"),
               "at least three samples")
  expect_error(removal(c(346,184),method="Schnute"),
               "at least three samples")
  ## Try using 3-pass method with not three catches
  expect_error(removal(c(346,184),method="Seber3"),
               "with three samples")
  expect_error(removal(c(346,184,49,12),method="Seber3"),
               "with three samples")
  ## Try using 2-pass method with not >2 catches
  expect_error(removal(c(346,184,49),method="Seber2"),
               "with two samples")
  expect_error(removal(c(346,184,49),method="RobsonRegier2"),
               "with two samples")
  ## Schnute warns if last of three passes is 0
  expect_warning(removal(c(4,2,0),method="Schnute"),
                 "The Schnute method will fail when the catch in")
  expect_warning(removal(c(400,200,0),method="Schnute"),
                 "The Schnute method will fail when the catch in")
  ## Burnham warns if only one fish caught, all fish caught on first pass, can't estimate
  expect_warning(removal(c(1,0,0),method="Burnham"),
                 "Total catch of one fish")
  expect_warning(removal(c(0,1,0),method="Burnham"),
                 "Total catch of one fish")
  expect_warning(removal(c(0,0,1),method="Burnham"),
                 "Total catch of one fish")
  expect_warning(removal(c(38,0,0),method="Burnham"),
                 "All fish captured on first pass")
  expect_warning(removal(c(38,38,38),method="Burnham"),
                 "failed to find")
  
  ## Errors in 2- and 3-pass methods if last catch is greater than first catch
  expect_warning(removal(c(184,346),method="Seber2"),
                 "results in model failure")
  expect_warning(removal(c(184,346),method="RobsonRegier2"),
                 "results in model failure")
  expect_warning(removal(c(49,184,346),method="Seber3"),
                 "results in model failure")
  ## Warnings if all catches are zeroes (except Carle-Strub)
  expect_warning(removal(c(0,0,0),method="Zippin"),
                 "model failure")
  expect_warning(removal(c(0,0,0),method="Schnute"),
                 "will fail")
  expect_warning(removal(c(0,0,0),method="Moran"),
                 "will fail")
  expect_warning(removal(c(0,0,0),method="Seber3"),
                 "model failure")
  expect_warning(removal(c(0,0,0),method="Burnham"),
                 "model failure")
  expect_warning(removal(c(0,0),method="Seber2"),
                 "model failure")
  expect_warning(removal(c(0,0),method="RobsonRegier2"),
                 "model failure")
  
  ## wrong parm in summary and confint
  tmp <- removal(c(346,184,49))
  expect_error(summary(tmp,parm="Derek"),
               "should be one of")
  expect_error(confint(tmp,parm="Derek"),
               "should be one of")
  expect_error(removal(c(346,184,49),conf.level=0),
               "must be between 0 and 1")
  expect_error(removal(c(346,184,49),conf.level=1),
               "must be between 0 and 1")
  expect_error(removal(c(346,184,49),conf.level="R"),
               "must be numeric")
  
  ## Check whole number,
  expect_warning(removal(c(346,184,49.1)),
                 "'catch' contains non-whole numbers.")
  expect_warning(removal(c(346.1,184.1,49.1)),
                 "'catch' contains non-whole numbers.")
  
  ## Bad data leads to failure of Zippin (from Carle-Strub (1978) example 2)
  expect_warning(removal(c(5,7,8),method="Zippin"),
                 "Zippin model failure")
  
  ## Chose "p1" summary for other than Schnute method
  tmp <- removal(c(45,11,18,8),method="Zippin")
  expect_error(summary(tmp,parm="p1"),
               "method does not use 'p1' parameter")
  expect_warning(summary(tmp,parm=c("p","p1")),
                 "method does not use 'p1' parameter")
  
  ## Chose only "p" CI for Moran or Schnute method
  tmp <- removal(c(45,11,18,8),method="Schnute")
  expect_error(confint(tmp,parm="p"),
               "Confidence intervals for 'p' can not be computed")
  
  ## Chose bad value for Tmult
  expect_error(removal(c(45,11,18,8),method="Moran",Tmult=0.9),
               "greater than 1")  
  expect_warning(removal(c(45,11,18,8),method="Moran",Tmult=1.2),
                 "try increasing")  
  ## NAs in catch vector
  expect_warning(removal(c(45,11,NA,8)),"'NA's removed from")
  
  ## Formula issues
  expect_error(removal(ct~pass,
                       data=data.frame(ct=c(12,7,2),pass=c("a","b","c"))),
               "'removal' formula must have only one variable")
  expect_error(removal(~pass,
                       data=data.frame(ct=c(12,7,2),pass=c("a","b","c"))),
               "must be numeric")
})  

test_that("removal() verbose= messages",{
  expect_message(summary(removal(c(38,26,12)),verbose=TRUE),
                 "Carle & Strub")
  expect_message(summary(removal(c(38,26,12),method="Moran"),verbose=TRUE),
                 "Moran")
  expect_message(summary(removal(c(38,26,12),method="Zippin"),verbose=TRUE),
                 "Zippin")
  expect_message(summary(removal(c(38,26,12),method="Schnute"),verbose=TRUE),
                 "Schnute")
  expect_message(summary(removal(c(38,26,12),method="Seber3"),verbose=TRUE),
                 "Seber")
  expect_message(summary(removal(c(38,26,12),method="Burnham"),verbose=TRUE),
                 "Burnham")
  expect_message(summary(removal(c(38,26),method="Seber2"),verbose=TRUE),
                 "Seber")
  expect_message(summary(removal(c(38,26),method="RobsonRegier2"),verbose=TRUE),
                 "Robson & Regier")
})   


## Test Output Types ----
test_that("removal.formula() same as removal.default()",{
  ct <- c(38,26,12)
  d <- data.frame(ct=ct)
  expect_equal(removal(~ct,data=d),removal(ct))
  expect_equal(removal(~ct,data=d,method="Moran"),removal(ct,method="Moran"))
  expect_equal(removal(~ct,data=d,method="Schnute"),removal(ct,method="Schnute"))
  # formula form matches default form ... only check default form after this
})

test_that("removal() return types",{
  expect_equal(class(removal(c(38,26,12))),"removal")
  expect_equal(class(removal(c(38,26,12),method="Moran")),"removal")
  expect_equal(class(removal(c(38,26,12),method="Zippin")),"removal")
  expect_equal(class(removal(c(38,26,12),method="Schnute")),"removal")
  expect_equal(class(removal(c(38,26,12),method="Seber3")),"removal")
  expect_equal(class(removal(c(38,26),method="Seber2")),"removal")
  expect_equal(class(removal(c(38,26),method="RobsonRegier2")),"removal")
  expect_equal(class(removal(c(38,26,12),method="Burnham")),"removal")
  # do one-dimensional data.frames and matrices work?
  expect_equal(class(removal(data.frame(c(38,26,12)))),"removal")
  expect_equal(class(removal(matrix(c(38,26,12),nrow=1))),"removal")
  expect_equal(class(removal(matrix(c(38,26,12),ncol=1))),"removal")

  # coef() results
  tmp <- removal(c(38,26,12))
  tmp2 <- coef(tmp)
  expect_true(is.vector(tmp2))
  expect_equal(class(tmp2),"numeric")
  expect_equal(length(tmp2),2)
  expect_equal(names(tmp2),c("No","p"))
  tmp2 <- coef(tmp,as.df=TRUE)
  expect_true(is.data.frame(tmp2))
  expect_equal(ncol(tmp2),2)
  expect_equal(nrow(tmp2),1)
  expect_equal(names(tmp2),c("No","p"))
  tmp2 <- coef(tmp,parm="p")
  expect_true(is.vector(tmp2))
  expect_equal(class(tmp2),"numeric")
  expect_equal(length(tmp2),1)
  expect_equal(names(tmp2),"p")
  tmp <- removal(c(38,26,12),method="Schnute")
  tmp2 <- coef(tmp)
  expect_true(is.vector(tmp2))
  expect_equal(class(tmp2),"numeric")
  expect_equal(length(tmp2),3)
  expect_equal(names(tmp2),c("No","p","p1"))
  tmp2 <- coef(tmp,as.df=TRUE)
  expect_true(is.data.frame(tmp2))
  expect_equal(ncol(tmp2),3)
  expect_equal(nrow(tmp2),1)
  expect_equal(names(tmp2),c("No","p","p1"))
  
  # summary() results
  tmp <- removal(c(38,26,12))
  tmp2 <- summary(tmp)
  expect_equal(class(tmp2),c("matrix","array"))
  expect_equal(nrow(tmp2),2)
  expect_equal(ncol(tmp2),2)
  expect_equal(rownames(tmp2),c("No","p"))
  expect_equal(colnames(tmp2),c("Estimate","Std. Error"))
  tmp2 <- summary(tmp,parm="p")
  expect_equal(class(tmp2),c("matrix","array"))
  expect_equal(nrow(tmp2),1)
  expect_equal(ncol(tmp2),2)
  expect_equal(rownames(tmp2),c("p"))
  expect_equal(colnames(tmp2),c("Estimate","Std. Error"))
  tmp2 <- summary(tmp,parm="No")
  expect_equal(class(tmp2),c("matrix","array"))
  expect_equal(nrow(tmp2),1)
  expect_equal(ncol(tmp2),2)
  expect_equal(rownames(tmp2),c("No"))
  expect_equal(colnames(tmp2),c("Estimate","Std. Error"))
  tmp <- removal(c(38,26,12),method="Schnute")
  tmp2 <- summary(tmp)
  expect_equal(class(tmp2),c("matrix","array"))
  expect_equal(nrow(tmp2),3)
  expect_equal(ncol(tmp2),2)
  expect_equal(rownames(tmp2),c("No","p","p1"))
  expect_equal(colnames(tmp2),c("Estimate","Std. Error"))
  tmp2 <- summary(tmp,parm="p1")
  expect_equal(class(tmp2),c("matrix","array"))
  expect_equal(nrow(tmp2),1)
  expect_equal(ncol(tmp2),2)
  expect_equal(rownames(tmp2),c("p1"))
  expect_equal(colnames(tmp2),c("Estimate","Std. Error"))
  tmp2 <- summary(tmp,parm=c("No","p"))
  expect_equal(class(tmp2),c("matrix","array"))
  expect_equal(nrow(tmp2),2)
  expect_equal(ncol(tmp2),2)
  expect_equal(rownames(tmp2),c("No","p"))
  expect_equal(colnames(tmp2),c("Estimate","Std. Error"))
  tmp <- removal(c(38,26),method="Seber2")
  tmp2 <- summary(tmp)
  expect_equal(class(tmp2),c("matrix","array"))
  expect_equal(nrow(tmp2),2)
  expect_equal(ncol(tmp2),2)
  expect_equal(rownames(tmp2),c("No","p"))
  expect_equal(colnames(tmp2),c("Estimate","Std. Error"))

  # summary() results ... as.df=TRUE
  tmp <- removal(c(38,26,12))
  tmp2 <- summary(tmp,as.df=TRUE)
  expect_true(is.data.frame(tmp2))
  expect_equal(nrow(tmp2),1)
  expect_equal(ncol(tmp2),4)
  expect_equal(names(tmp2),c("No","No.se","p","p.se"))
  tmp2 <- summary(tmp,parm="p",as.df=TRUE)
  expect_true(is.data.frame(tmp2))
  expect_equal(nrow(tmp2),1)
  expect_equal(ncol(tmp2),2)
  expect_equal(names(tmp2),c("p","p.se"))
  tmp2 <- summary(tmp,parm="No",as.df=TRUE)
  expect_true(is.data.frame(tmp2))
  expect_equal(nrow(tmp2),1)
  expect_equal(ncol(tmp2),2)
  expect_equal(names(tmp2),c("No","No.se"))
  tmp <- removal(c(38,26,12),method="Schnute")
  tmp2 <- summary(tmp,as.df=TRUE)
  expect_true(is.data.frame(tmp2))
  expect_equal(nrow(tmp2),1)
  expect_equal(ncol(tmp2),6)
  expect_equal(names(tmp2),c("No","No.se","p","p.se","p1","p1.se"))
  tmp2 <- summary(tmp,parm="p1",as.df=TRUE)
  expect_true(is.data.frame(tmp2))
  expect_equal(nrow(tmp2),1)
  expect_equal(ncol(tmp2),2)
  expect_equal(names(tmp2),c("p1","p1.se"))
  tmp2 <- summary(tmp,parm=c("p","p1"),as.df=TRUE)
  expect_true(is.data.frame(tmp2))
  expect_equal(nrow(tmp2),1)
  expect_equal(ncol(tmp2),4)
  expect_equal(names(tmp2),c("p","p.se","p1","p1.se"))

  # confint() results
  tmp <- removal(c(38,26,12))
  tmp2 <- confint(tmp)
  expect_equal(class(tmp2),c("matrix","array"))
  expect_equal(nrow(tmp2),2)
  expect_equal(ncol(tmp2),2)
  expect_equal(rownames(tmp2),c("No","p"))
  expect_equal(colnames(tmp2),c("95% LCI","95% UCI"))
  tmp2 <- confint(tmp,parm="p")
  expect_equal(class(tmp2),c("matrix","array"))
  expect_equal(nrow(tmp2),1)
  expect_equal(ncol(tmp2),2)
  expect_equal(rownames(tmp2),c("p"))
  expect_equal(colnames(tmp2),c("95% LCI","95% UCI"))
  tmp2 <- confint(tmp,parm="No")
  expect_equal(class(tmp2),c("matrix","array"))
  expect_equal(nrow(tmp2),1)
  expect_equal(ncol(tmp2),2)
  expect_equal(rownames(tmp2),c("No"))
  expect_equal(colnames(tmp2),c("95% LCI","95% UCI"))
  tmp <- removal(c(38,26,12),method="Schnute")
  tmp2 <- suppressMessages(confint(tmp))
  expect_equal(class(tmp2),c("matrix","array"))
  expect_equal(nrow(tmp2),3)
  expect_equal(ncol(tmp2),2)
  expect_equal(rownames(tmp2),c("No","p","p1"))
  expect_equal(colnames(tmp2),c("95% LCI","95% UCI"))
  tmp2 <- confint(tmp,parm="No")
  expect_equal(class(tmp2),c("matrix","array"))
  expect_equal(nrow(tmp2),1)
  expect_equal(ncol(tmp2),2)
  expect_equal(rownames(tmp2),c("No"))
  expect_equal(colnames(tmp2),c("95% LCI","95% UCI"))
  tmp <- removal(c(38,26,12),method="Schnute")
  tmp2 <- suppressMessages(confint(tmp,incl.est=TRUE))
  expect_equal(class(tmp2),c("matrix","array"))
  expect_equal(nrow(tmp2),3)
  expect_equal(ncol(tmp2),3)
  expect_equal(rownames(tmp2),c("No","p","p1"))
  expect_equal(colnames(tmp2),c("Est","95% LCI","95% UCI"))
  tmp2 <- confint(tmp,parm="No",incl.est=TRUE)
  expect_equal(class(tmp2),c("matrix","array"))
  expect_equal(nrow(tmp2),1)
  expect_equal(ncol(tmp2),3)
  expect_equal(rownames(tmp2),"No")
  expect_equal(colnames(tmp2),c("Est","95% LCI","95% UCI"))
  tmp <- removal(c(38,26),method="Seber2")
  tmp2 <- confint(tmp)
  expect_equal(class(tmp2),c("matrix","array"))
  expect_equal(nrow(tmp2),2)
  expect_equal(ncol(tmp2),2)
  expect_equal(rownames(tmp2),c("No","p"))
  expect_equal(colnames(tmp2),c("95% LCI","95% UCI"))
  tmp2 <- confint(tmp,parm="p")
  expect_equal(class(tmp2),c("matrix","array"))
  expect_equal(nrow(tmp2),1)
  expect_equal(ncol(tmp2),2)
  expect_equal(rownames(tmp2),c("p"))
  expect_equal(colnames(tmp2),c("95% LCI","95% UCI"))
  tmp2 <- confint(tmp,parm="No")
  expect_equal(class(tmp2),c("matrix","array"))
  expect_equal(nrow(tmp2),1)
  expect_equal(ncol(tmp2),2)
  expect_equal(rownames(tmp2),c("No"))
  expect_equal(colnames(tmp2),c("95% LCI","95% UCI"))
  ## Confint with as.df=TRUE
  tmp <- removal(c(38,26),method="Seber2")
  tmp2 <- confint(tmp,incl.est=TRUE,as.df=TRUE)
  expect_true(is.data.frame(tmp2))
  expect_equal(nrow(tmp2),1)
  expect_equal(ncol(tmp2),6)
  expect_equal(names(tmp2),c("No","No.LCI","No.UCI","p","p.LCI","p.UCI"))
  tmp2 <- confint(tmp,parm="No",incl.est=TRUE,as.df=TRUE)
  expect_true(is.data.frame(tmp2))
  expect_equal(nrow(tmp2),1)
  expect_equal(ncol(tmp2),3)
  expect_equal(names(tmp2),c("No","No.LCI","No.UCI"))
  
  ## What if catches are all zeroes
  suppressWarnings(
    expect_true(all(is.na(removal(c(0,0,0),method="Zippin")$est))))
  suppressWarnings(
    expect_true(all(is.na(removal(c(0,0,0),method="Schnute")$est))))
  suppressWarnings(
    expect_true(all(is.na(removal(c(0,0,0),method="Moran")$est))))
  suppressWarnings(
    expect_true(all(is.na(removal(c(0,0,0),method="Seber3")$est))))
  suppressWarnings(
    expect_true(all(is.na(removal(c(0,0),method="Seber2")$est))))
  suppressWarnings(
    expect_true(all(is.na(removal(c(0,0),method="RobsonRegier2")$est))))
  suppressWarnings(
    expect_true(all(is.na(removal(c(0,0,0),method="Burnham")$est))))
})


## Validate Results ----
test_that("removal with 'CarleStrub' matches Carle-Strub (1978) examples",{
  tmp <- summary(removal(c(38,26,12)))
  expect_equal(round(tmp["No","Estimate"],0),91)
  expect_equal(round(tmp["No","Std. Error"],1),9.7)
  expect_equal(round(tmp["p","Estimate"],3),0.444)
  
  tmp <- summary(removal(c(5,7,8)))
  expect_equal(round(tmp["No","Estimate"],0),44)
  expect_equal(round(tmp["p","Estimate"],3),0.174)
})

test_that("removal with 'CarleStrub' matches Cowx (1983) page 77",{
  tmp <- summary(removal(c(72,56,46,30,24)))
  expect_equal(round(tmp["No","Estimate"],0),298)
  expect_equal(round(tmp["p","Estimate"],3),0.250)
  # SE does not match
  #expect_equal(round(tmp["No","Std. Error"],1),23.62)
  
  tmp <- summary(removal(c(8,23,17,8,6)))
  expect_equal(round(tmp["No","Estimate"],0),95)
  expect_equal(round(tmp["p","Estimate"],3),0.187)
})

test_that("removal with 'Seber3' matches Cowx (1983) page 75",{
  tmp <- summary(removal(c(72,56,46),method="Seber3"))
  expect_equal(round(tmp["No","Estimate"],0),353)
})

test_that("removal with 'Seber2' matches Cowx (1983) page 75",{
  tmp <- summary(removal(c(72,56),method="Seber2"))
  expect_equal(round(tmp["No","Estimate"],0),324)
  expect_equal(round(tmp["No","Std. Error"],2),178.19)
  expect_equal(round(tmp["p","Estimate"],2),0.22)
})

test_that("removal with 'Seber2' matches Seber(2012) example 7.4",{
  tmp <- summary(removal(c(79,28),method="Seber2"))
  expect_equal(round(tmp["No","Estimate"],0),122)
  expect_equal(round(tmp["No","Std. Error"],1),8.8)
  expect_equal(round(tmp["p","Estimate"],2),0.65)
})

test_that("removal with 'RobsonRegier2' matches Cowx (1983) page 75",{
  tmp <- summary(removal(c(72,56),method="RobsonRegier2"))
  # used ceiling because of weird round issue
  expect_equal(ceiling(tmp["No","Estimate"]),321)
  expect_equal(round(tmp["No","Std. Error"],2),178.19)
})


test_that("removal with 'Moran' matches Schnute (1983)",{
  data(BrookTroutNEWP1,package="FSAdata")
  Ns <- ps <- LHs <- NLCI <- NUCI <- numeric(nrow(BrookTroutNEWP1))
  for (i in seq_len(nrow(BrookTroutNEWP1))) {
    tmp <- removal(as.numeric(BrookTroutNEWP1[i,c("first","second",
                                                  "third","fourth")]),
                   method="Moran")
    Ns[i] <- round(tmp$est[["No"]],1)
    ps[i] <- round(tmp$est[["p"]],2)
    LHs[i] <- round(tmp$min.nlogLH,2)
    suppressMessages(tmp <- confint(tmp))
    NLCI[i] <- tmp["No","95% LCI"]
    NUCI[i] <- tmp["No","95% UCI"]
  }
  ## check point estimates
  tmp <- cbind(sample=seq_len(nrow(BrookTroutNEWP1)),Ns,ps,LHs,
               BrookTroutNEWP1[,c("Moran.N","Moran.p","Moran.LH")])
  ## perfect matches
  expect_equal(tmp[,"Ns"],BrookTroutNEWP1$Moran.N[])
  expect_equal(tmp[,"ps"],BrookTroutNEWP1$Moran.p[])
  expect_equal(tmp[,"LHs"],BrookTroutNEWP1$Moran.LH[])
  ## Check CIs (off by no more than 0.1 in a small handful of the UCIs)
  tmp <- cbind(sample=seq_len(nrow(BrookTroutNEWP1)),NLCI,NUCI,
               BrookTroutNEWP1[,c("Moran.NLCI","Moran.NUCI")])
  expect_true(all(abs(tmp[,2:3]-tmp[,4:5])<=0.1001))
})

test_that("removal with 'Schnute' matches Schnute (1983)",{
  data(BrookTroutNEWP1,package="FSAdata")
  Ns <- p1s <- ps <- LHs <- NLCI <- NUCI <- numeric(nrow(BrookTroutNEWP1))
  for (i in seq_len(nrow(BrookTroutNEWP1))) {
    tmp <- removal(as.numeric(BrookTroutNEWP1[i,c("first","second",
                                                  "third","fourth")]),
                   method="Schnute")
    Ns[i] <- round(tmp$est[["No"]],1)
    p1s[i] <- round(tmp$est[["p1"]],2)
    ps[i] <- round(tmp$est[["p"]],2)
    LHs[i] <- round(tmp$min.nlogLH,2)
    suppressMessages(tmp <- confint(tmp))
    NLCI[i] <- tmp["No","95% LCI"]
    NUCI[i] <- tmp["No","95% UCI"]
  }
  ## check point estimates
  tmp <- cbind(sample=seq_len(nrow(BrookTroutNEWP1)),Ns,p1s,ps,LHs,
               BrookTroutNEWP1[,c("Schnute.N","Schnute.p1",
                                  "Schnute.p","Schnute.LH")])
  ## perfect matches except sample 5 N is off by 0.1
  expect_equal(tmp[-5,"Ns"],BrookTroutNEWP1$Schnute.N[-5])
  expect_equal(tmp[,"p1s"],BrookTroutNEWP1$Schnute.p1[])
  expect_equal(tmp[,"ps"],BrookTroutNEWP1$Schnute.p[])
  expect_equal(tmp[,"LHs"],BrookTroutNEWP1$Schnute.LH[])
  ## Check CIs (off by no more than 0.1)
  tmp <- cbind(sample=seq_len(nrow(BrookTroutNEWP1)),NLCI,NUCI,
               BrookTroutNEWP1[,c("Schnute.NLCI","Schnute.NUCI")])
  expect_true(all(abs(tmp[,2:3]-tmp[,4:5])<=0.1001,na.rm=TRUE))
})

test_that("removal with 'Burnham' match results from (Van Deventer 1989) page 13",{
  tmp <- removal(c(124,61,35,14),method="Burnham",CIMicroFish=TRUE)
  ## check point estimates
  tmp2 <- summary(tmp)
  expect_equal(round(tmp2["No","Estimate"],0),249)
  expect_equal(round(tmp2["No","Std. Error"],3),6.164)
  expect_equal(round(tmp2["p","Estimate"],3),0.501)
  expect_equal(round(tmp2["p","Std. Error"],3),0.035)
  ## check CIs
  tmp2 <- confint(tmp)
  expect_equal(round(as.numeric(tmp2["No",]),3),c(237,261))
  expect_equal(round(as.numeric(tmp2["p",]),3),c(0.432,0.570))
  
})