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