R Under development (unstable) (2024-12-07 r87428 ucrt) -- "Unsuffered Consequences"
Copyright (C) 2024 The R Foundation for Statistical Computing
Platform: x86_64-w64-mingw32/x64

R is free software and comes with ABSOLUTELY NO WARRANTY.
You are welcome to redistribute it under certain conditions.
Type 'license()' or 'licence()' for distribution details.

R is a collaborative project with many contributors.
Type 'contributors()' for more information and
'citation()' on how to cite R or R packages in publications.

Type 'demo()' for some demos, 'help()' for on-line help, or
'help.start()' for an HTML browser interface to help.
Type 'q()' to quit R.

> if (file.exists("_options.R")) source("_options.R")
> library(panelPomp,quietly=TRUE)
> 
> TESTS_PASS <- NULL
> test <- function(expr1,expr2,all="TESTS_PASS",env=parent.frame(),...)
+   panelPomp:::test(expr1,expr2,all=all,env=env,...)
> 
> ## basic checks for logLik, matrix-method
> repMargin <- 1
> unit1.ll <- 1
> unit2.ll <- 10
> object <- ulls <- matrix(c(unit1.ll,unit1.ll,unit2.ll,unit2.ll),nr=2,nc=2,
+                          byrow=ifelse(repMargin==1,FALSE,TRUE))
> exp.ll <- unit1.ll + unit2.ll
> exp.se <- 0
> 
> chck.avl <- logLik(ulls,repMargin=repMargin,first="aver",aver="logmeanexp",se=TRUE)
> chck.avm <- logLik(ulls,repMargin=repMargin,first="aver",aver="mean",se=TRUE)
> chck.agl <- logLik(ulls,repMargin=repMargin,first="aggr",aver="logmeanexp",se=TRUE)
> chck.agm <- logLik(ulls,repMargin=repMargin,first="aggr",aver="mean",se=TRUE)
> 
> test(TRUE,all(c(
+   chck.avl==c(exp.ll,exp.se),
+   chck.avm==c(exp.ll,exp.se),
+   chck.agl==c(exp.ll,exp.se),
+   chck.agm==c(exp.ll,exp.se))))
[1] TRUE
> 
> ## test se=FALSE
> test(c("","se"),
+      names(
+        logLik(ulls,repMargin=repMargin,first="aver",aver="logmeanexp",se=TRUE)
+      ))
[1] TRUE
> test(unit1.ll+unit2.ll,
+      logLik(ulls,repMargin=repMargin,first="aver",aver="logmeanexp",se=FALSE))
[1] TRUE
> 
> ## test first='aver' & aver='logmeanexp' for logLik, matrix-method
> repMargin <- 2
> unit1.ll <- 2
> unit2.ll <- 10
> object <- ulls <- matrix(
+   c(log(0.5*exp(unit1.ll)),log(1.5*exp(unit1.ll)),unit2.ll,unit2.ll),
+   nr=2,nc=2,byrow=ifelse(repMargin==1,FALSE,TRUE)
+ )
> chck <- logLik(ulls,repMargin=repMargin,first="aver",aver="logmeanexp",se=TRUE)
> test(unname(chck[1]),c(unit1.ll + unit2.ll))
[1] TRUE
> 
> ## test aver='mean' for logLik, matrix-method
> repMargin <- 1
> unit1.ll <- 1
> unit2.ll <- 10
> object <- ulls <- matrix(
+   c(unit1.ll+0.5,unit1.ll-0.5,unit2.ll,unit2.ll),
+   nr=2,nc=2,byrow=ifelse(repMargin==1,FALSE,TRUE)
+ )
> chck <- logLik(ulls,repMargin=repMargin,first="aggr",aver="mean",se=TRUE)
> test(unname(chck[1]),c(unit1.ll + unit2.ll))
[1] TRUE
> 
> ## test panel_logmeanexp
> repMargin <- 1
> ul1 <- 1
> ul2 <- 10
> uls <- matrix(c(ul1,ul1,ul2,ul2),nr=2,byrow=ifelse(repMargin==1,FALSE,TRUE))
> expected_output <- setNames(c(ul1+ul2,0),nm=c("","se"))
> 
> test(panel_logmeanexp(uls,MARGIN=ifelse(repMargin==1,2,1),se=FALSE),
+      unname(expected_output[1]))
[1] TRUE
> test(panel_logmeanexp(uls,MARGIN=ifelse(repMargin==1,2,1),se=TRUE),
+      expected_output)
[1] TRUE
> 
> ## check whether all tests passed
> all(get(eval(formals(test))$all))
[1] TRUE
> if (!all(get(eval(formals(test))$all))) stop("Not all tests passed!")
> 
> proc.time()
   user  system elapsed 
   0.71    0.17    0.87