*** V2.7.1 ***

- Fixes to several functions including bmd, bmdMA, bmdBoot and many more
- qplotBmd works for multiple curves
- Added tests using testthat

Author: JensBaalkilde

DESCRIPTION

@@ -1,16 +1,24 @@
 Package: bmd
 Type: Package
 Title: Benchmark dose estimation for dose-response data
-Version: 2.6.7
-Date: 2024-12-11
+Version: 2.7.1
+Date: 2025-03-24
 Author: Signe M.Jensen, Christian Ritz and Jens Riis Baalkilde
 Maintainer: Signe M. Jensen <smj@plen.ku.dk>
 Description: Benchmark dose analysis for continuous, quantal, count and ordinal dose-response data
 Imports: drc, ggplot2, dplyr
 Suggests: 
+    sandwich,
     CVXR,
     multcomp,
     gridExtra,
+    isotone,
+    reshape2,
+    dplyr,
+    car,
+    Matrix,
+    RLRsim,
+    scales,
     testthat (>= 3.0.0)
 License: GPL
 Encoding: UTF-8

NAMESPACE

@@ -1,11 +1,13 @@
 import(drc, ggplot2, dplyr)
 export(bmd, bmdBoot, bmdIso, bmdIsoBoot, PAV, bmdMA, bootDataGen, bmdMACurve, BCa, invBmd, expandBinomial, 
-       getStackingWeights, drmOrdinal, bmdOrdinal, bmdOrdinalMA, expandOrdinal, bootDataGenOrdinal,
+       getStackingWeights, drmOrdinal, bmdOrdinal, bmdOrdinalMA,
+       expandOrdinal, bootDataGenOrdinal,
        qplotDrc, qplotBmd, MACurve, monotonicityTest, trendTest, bmdHetVar, drmHetVar)
 
 ## S3 methods
 S3method(logLik, drcOrdinal)
 S3method(AIC, drcOrdinal)
+S3method(BIC, drcOrdinal)
 S3method(plot, drcOrdinal)
 S3method(print, drcOrdinal)
 S3method(print, bmdOrdinal)

R/BIC.drcOrdinal.R

@@ -0,0 +1,5 @@
+BIC.drcOrdinal <- function(object, epsilon = 10^(-16)){
+  n.parameters <- sum(sapply(object$drmList, function(mod) length(mod$coefficients)))
+  n.obs <- sum(object$data[[object$weights]])
+  n.parameters * log(n.obs) - 2 * logLik(object, epsilon) 
+}