R Under development (unstable) (2024-09-06 r87103 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.

> # check reproducibility over a range of randomly chosen seeds
> if (file.exists("_options.R")) source("_options.R")
> library(panelPomp,quietly=TRUE)
> 
> ppo <- panelRandomWalk(U=3,N=4)
> 
> seeds <- c(21125715L,81902071L,86680005L,29971689L,73376915L)
> 
> for (s in seq_along(seeds)) {
+   set.seed(seeds[s])
+   pf <- pfilter(ppo,Np=10)
+   print(unitLogLik(pf))
+   print(logLik(pf))
+ }
      rw1       rw2       rw3 
-9.693242 -7.958418 -5.702933 
[1] -23.35459
      rw1       rw2       rw3 
-9.939198 -9.169043 -5.679824 
[1] -24.78806
      rw1       rw2       rw3 
-8.860156 -7.318119 -6.049635 
[1] -22.22791
      rw1       rw2       rw3 
-7.809090 -8.144534 -5.255597 
[1] -21.20922
      rw1       rw2       rw3 
-8.512685 -6.536669 -6.036917 
[1] -21.08627
> 
> for (s in seq_along(seeds)) {
+   set.seed(seeds[s])
+   mf <- mif2(ppo,Np=10,rw.sd=rw_sd(sigmaX=0.05,X.0=0.5),
+              cooling.fraction.50=0.5,cooling.type="geometric")
+   print(mf@shared)
+   print(mf@specific)
+   print(unitLogLik(mf))
+   print(logLik(mf))
+   print(traces(mf))
+ }
  sigmaX   sigmaY 
1.078101 1.000000 
        unit
variable       rw1       rw2      rw3
     X.0 0.9995142 0.6630938 1.546127
       rw1        rw2        rw3 
-10.308758  -7.274029  -5.310702 
[1] -22.89349
         variable
iteration    loglik   sigmaX sigmaY unitLoglik[rw1]  X.0[rw1] unitLoglik[rw2]
        0 -22.89349 1.000000      1       -10.30876 1.0000000       -7.274029
        1        NA 1.078101      1              NA 0.9995142              NA
         variable
iteration  X.0[rw2] unitLoglik[rw3] X.0[rw3]
        0 1.0000000       -5.310702 1.000000
        1 0.6630938              NA 1.546127
   sigmaX    sigmaY 
0.8366933 1.0000000 
        unit
variable       rw1     rw2      rw3
     X.0 0.1609146 1.56302 1.238118
      rw1       rw2       rw3 
-8.853153 -8.303707 -5.563768 
[1] -22.72063
         variable
iteration    loglik    sigmaX sigmaY unitLoglik[rw1]  X.0[rw1] unitLoglik[rw2]
        0 -22.72063 1.0000000      1       -8.853153 1.0000000       -8.303707
        1        NA 0.8366933      1              NA 0.1609146              NA
         variable
iteration X.0[rw2] unitLoglik[rw3] X.0[rw3]
        0  1.00000       -5.563768 1.000000
        1  1.56302              NA 1.238118
   sigmaX    sigmaY 
0.7177845 1.0000000 
        unit
variable       rw1      rw2      rw3
     X.0 0.3572861 1.608191 1.214683
       rw1        rw2        rw3 
-10.032511  -8.120632  -5.503718 
[1] -23.65686
         variable
iteration    loglik    sigmaX sigmaY unitLoglik[rw1]  X.0[rw1] unitLoglik[rw2]
        0 -23.65686 1.0000000      1       -10.03251 1.0000000       -8.120632
        1        NA 0.7177845      1              NA 0.3572861              NA
         variable
iteration X.0[rw2] unitLoglik[rw3] X.0[rw3]
        0 1.000000       -5.503718 1.000000
        1 1.608191              NA 1.214683
   sigmaX    sigmaY 
0.9576311 1.0000000 
        unit
variable       rw1       rw2      rw3
     X.0 0.2445654 0.6539408 1.133491
      rw1       rw2       rw3 
-9.367004 -8.050974 -5.949910 
[1] -23.36789
         variable
iteration    loglik    sigmaX sigmaY unitLoglik[rw1]  X.0[rw1] unitLoglik[rw2]
        0 -23.36789 1.0000000      1       -9.367004 1.0000000       -8.050974
        1        NA 0.9576311      1              NA 0.2445654              NA
         variable
iteration  X.0[rw2] unitLoglik[rw3] X.0[rw3]
        0 1.0000000        -5.94991 1.000000
        1 0.6539408              NA 1.133491
   sigmaX    sigmaY 
0.8674928 1.0000000 
        unit
variable      rw1      rw2      rw3
     X.0 1.352745 1.253761 1.459401
      rw1       rw2       rw3 
-9.720848 -6.949583 -5.832371 
[1] -22.5028
         variable
iteration   loglik    sigmaX sigmaY unitLoglik[rw1] X.0[rw1] unitLoglik[rw2]
        0 -22.5028 1.0000000      1       -9.720848 1.000000       -6.949583
        1       NA 0.8674928      1              NA 1.352745              NA
         variable
iteration X.0[rw2] unitLoglik[rw3] X.0[rw3]
        0 1.000000       -5.832371 1.000000
        1 1.253761              NA 1.459401
> 
> proc.time()
   user  system elapsed 
   0.76    0.07    2.50