inst/classic-bugs/vol2/beetles/bench-test1.R

In nimble: MCMC, Particle Filtering, and Programmable Hierarchical Modeling

"benchstats" <-
structure(c(-60.7831272900001, 34.3079583900001, 3.56745360999999, 
9.94903490800002, 22.4989411700000, 33.89292992, 50.0425105499999, 
53.2141230900001, 59.1443245600001, 58.6841452499999, 13.21342134, 
5.15217918889568, 2.89643983703556, 0.95543584548527, 1.68101350129359, 
2.0897844158153, 1.75364152690944, 1.65285717332056, 1.10690742807124, 
0.737680372356471, 0.427056870264935, 1.99912021164428, 0.0515217918889568, 
0.0289643983703556, 0.0095543584548527, 0.0168101350129359, 0.020897844158153, 
0.0175364152690944, 0.0165285717332056, 0.0110690742807124, 0.00737680372356471, 
0.00427056870264935, 0.0199912021164428, 0.0749820368184437, 
0.0424238958019385, 0.0120982114296435, 0.0201087388192881, 0.0239164039838735, 
0.0224611804560939, 0.0243408587514327, 0.0170645630747179, 0.0114256340298251, 
0.00657670318302355, 0.0301557929820285), .Dim = as.integer(c(11, 
4)), .Dimnames = list(c("alpha", "beta", "r.hat[1]", "r.hat[2]", 
"r.hat[3]", "r.hat[4]", "r.hat[5]", "r.hat[6]", "r.hat[7]", "r.hat[8]", 
"D"), c("Mean", "SD", "Naive SE", "Time-series SE")))