cameron@macbookpro:mv_extremes$ . start_from_scratch.sh rm -rf src/*.o src/*.mod fit_model cameron@macbookpro:mv_extremes$ . run_all.sh (cd src; make all) gfortran -O2 -c functions.f90 gfortran -O2 -c model.f90 gfortran -O2 -c mfu_sampler.f90 gfortran -O2 -c flib.f gfortran -O2 -c flib_blas.f gfortran -O2 -framework Accelerate fit_model.f90 -o fit_model src/*.o -Isrc usage: cmdline [OPTIONS] options: -h, --help print this message -s, --sample-file file name of the file to save samples [default=samples] -l, --lp-file file name of the file to save samples [default=none] -i, -iterations n number of iterations to run [default=20000] -r, -refresh n print after each refresh iterations [default=1] --no-sample-output do not write out any samples (useful for model selection). Sample file is: output/bayes_joint/samples/samples_01 log posterior file is: output/bayes_joint/samples/lp_01 Iterations: 5000 Refresh: 100 Initial lp: -1007.9657794324854 Sampling started. Iteration 100 of 5000, 1.928s, lp: 100.30 Iteration 200 of 5000, 1.649s, lp: 99.20 Iteration 300 of 5000, 1.341s, lp: 98.63 Iteration 400 of 5000, 1.167s, lp: 105.09 Iteration 500 of 5000, 1.194s, lp: 101.65 Iteration 600 of 5000, 1.108s, lp: 104.39 Iteration 700 of 5000, 1.094s, lp: 99.07 Iteration 800 of 5000, 1.035s, lp: 95.75 Iteration 900 of 5000, 1.097s, lp: 101.36 Iteration 1000 of 5000, 1.038s, lp: 106.96 Iteration 1100 of 5000, 0.984s, lp: 100.79 Iteration 1200 of 5000, 1.002s, lp: 100.25 Iteration 1300 of 5000, 1.005s, lp: 101.77 Iteration 1400 of 5000, 1.004s, lp: 92.49 Iteration 1500 of 5000, 0.998s, lp: 104.86 Iteration 1600 of 5000, 0.996s, lp: 103.34 Iteration 1700 of 5000, 1.056s, lp: 103.60 Iteration 1800 of 5000, 0.992s, lp: 102.04 Iteration 1900 of 5000, 1.041s, lp: 106.12 Iteration 2000 of 5000, 1.050s, lp: 98.61 Iteration 2100 of 5000, 1.014s, lp: 93.44 Iteration 2200 of 5000, 0.993s, lp: 101.62 Iteration 2300 of 5000, 1.043s, lp: 100.25 Iteration 2400 of 5000, 1.001s, lp: 87.97 Iteration 2500 of 5000, 1.017s, lp: 102.17 Iteration 2600 of 5000, 1.015s, lp: 96.54 Iteration 2700 of 5000, 1.033s, lp: 108.24 Iteration 2800 of 5000, 1.112s, lp: 96.93 Iteration 2900 of 5000, 0.953s, lp: 94.80 Iteration 3000 of 5000, 0.947s, lp: 98.46 Iteration 3100 of 5000, 0.973s, lp: 100.48 Iteration 3200 of 5000, 1.002s, lp: 103.24 Iteration 3300 of 5000, 0.948s, lp: 103.36 Iteration 3400 of 5000, 0.991s, lp: 104.17 Iteration 3500 of 5000, 1.013s, lp: 97.52 Iteration 3600 of 5000, 0.976s, lp: 95.01 Iteration 3700 of 5000, 0.941s, lp: 94.03 Iteration 3800 of 5000, 0.968s, lp: 97.18 Iteration 3900 of 5000, 0.944s, lp: 96.05 Iteration 4000 of 5000, 0.954s, lp: 93.17 Iteration 4100 of 5000, 0.948s, lp: 101.35 Iteration 4200 of 5000, 0.940s, lp: 91.19 Iteration 4300 of 5000, 0.937s, lp: 98.09 Iteration 4400 of 5000, 0.992s, lp: 101.31 Iteration 4500 of 5000, 0.916s, lp: 101.10 Iteration 4600 of 5000, 0.955s, lp: 104.90 Iteration 4700 of 5000, 0.996s, lp: 106.11 Iteration 4800 of 5000, 1.007s, lp: 97.89 Iteration 4900 of 5000, 0.946s, lp: 101.88 Iteration 5000 of 5000, 0.934s, lp: 94.13 Computation time: 52.2 seconds Sample file is: output/bayes_joint/samples/samples_02 log posterior file is: output/bayes_joint/samples/lp_02 Iterations: 5000 Refresh: 100 Initial lp: -1007.9657794324854 Sampling started. Iteration 100 of 5000, 2.134s, lp: 103.99 Iteration 200 of 5000, 1.431s, lp: 100.69 Iteration 300 of 5000, 1.288s, lp: 91.56 Iteration 400 of 5000, 1.232s, lp: 91.31 Iteration 500 of 5000, 1.215s, lp: 101.58 Iteration 600 of 5000, 1.208s, lp: 100.40 Iteration 700 of 5000, 1.158s, lp: 104.48 Iteration 800 of 5000, 1.145s, lp: 92.16 Iteration 900 of 5000, 1.118s, lp: 99.56 Iteration 1000 of 5000, 1.090s, lp: 99.07 Iteration 1100 of 5000, 1.035s, lp: 96.26 Iteration 1200 of 5000, 1.053s, lp: 102.13 Iteration 1300 of 5000, 1.054s, lp: 96.67 Iteration 1400 of 5000, 1.026s, lp: 100.87 Iteration 1500 of 5000, 1.072s, lp: 104.57 Iteration 1600 of 5000, 1.028s, lp: 96.10 Iteration 1700 of 5000, 1.038s, lp: 100.81 Iteration 1800 of 5000, 1.043s, lp: 102.06 Iteration 1900 of 5000, 1.053s, lp: 96.90 Iteration 2000 of 5000, 1.170s, lp: 101.27 Iteration 2100 of 5000, 1.176s, lp: 106.45 Iteration 2200 of 5000, 1.116s, lp: 100.19 Iteration 2300 of 5000, 1.087s, lp: 97.71 Iteration 2400 of 5000, 1.140s, lp: 104.48 Iteration 2500 of 5000, 1.111s, lp: 96.84 Iteration 2600 of 5000, 1.141s, lp: 108.51 Iteration 2700 of 5000, 1.097s, lp: 99.46 Iteration 2800 of 5000, 1.045s, lp: 100.29 Iteration 2900 of 5000, 0.992s, lp: 100.85 Iteration 3000 of 5000, 1.009s, lp: 99.45 Iteration 3100 of 5000, 0.984s, lp: 101.59 Iteration 3200 of 5000, 0.960s, lp: 96.21 Iteration 3300 of 5000, 0.972s, lp: 99.87 Iteration 3400 of 5000, 0.997s, lp: 101.27 Iteration 3500 of 5000, 0.981s, lp: 105.57 Iteration 3600 of 5000, 0.985s, lp: 91.32 Iteration 3700 of 5000, 0.974s, lp: 98.78 Iteration 3800 of 5000, 1.008s, lp: 97.35 Iteration 3900 of 5000, 0.996s, lp: 103.93 Iteration 4000 of 5000, 0.999s, lp: 91.11 Iteration 4100 of 5000, 0.969s, lp: 95.45 Iteration 4200 of 5000, 0.952s, lp: 93.55 Iteration 4300 of 5000, 0.988s, lp: 99.48 Iteration 4400 of 5000, 0.936s, lp: 96.58 Iteration 4500 of 5000, 0.941s, lp: 97.14 Iteration 4600 of 5000, 0.993s, lp: 101.12 Iteration 4700 of 5000, 1.002s, lp: 98.31 Iteration 4800 of 5000, 1.008s, lp: 96.29 Iteration 4900 of 5000, 0.986s, lp: 105.25 Iteration 5000 of 5000, 0.959s, lp: 93.66 Computation time: 54.1 seconds Sample file is: output/bayes_joint/samples/samples_03 log posterior file is: output/bayes_joint/samples/lp_03 Iterations: 5000 Refresh: 100 Initial lp: -1007.9657794324854 Sampling started. Iteration 100 of 5000, 1.041s, lp: 97.18 Iteration 200 of 5000, 0.993s, lp: 90.43 Iteration 300 of 5000, 0.962s, lp: 103.07 Iteration 400 of 5000, 0.939s, lp: 102.81 Iteration 500 of 5000, 0.933s, lp: 100.84 Iteration 600 of 5000, 0.939s, lp: 96.24 Iteration 700 of 5000, 0.950s, lp: 94.83 Iteration 800 of 5000, 0.925s, lp: 98.90 Iteration 900 of 5000, 0.934s, lp: 101.65 Iteration 1000 of 5000, 0.935s, lp: 98.46 Iteration 1100 of 5000, 0.934s, lp: 99.74 Iteration 1200 of 5000, 0.929s, lp: 101.26 Iteration 1300 of 5000, 0.931s, lp: 100.34 Iteration 1400 of 5000, 0.933s, lp: 102.67 Iteration 1500 of 5000, 0.899s, lp: 101.60 Iteration 1600 of 5000, 0.905s, lp: 91.68 Iteration 1700 of 5000, 0.897s, lp: 102.14 Iteration 1800 of 5000, 0.915s, lp: 96.96 Iteration 1900 of 5000, 0.921s, lp: 95.55 Iteration 2000 of 5000, 0.918s, lp: 97.26 Iteration 2100 of 5000, 0.918s, lp: 102.43 Iteration 2200 of 5000, 0.892s, lp: 93.50 Iteration 2300 of 5000, 0.935s, lp: 103.73 Iteration 2400 of 5000, 1.021s, lp: 101.62 Iteration 2500 of 5000, 1.031s, lp: 112.83 Iteration 2600 of 5000, 1.080s, lp: 100.91 Iteration 2700 of 5000, 1.020s, lp: 101.56 Iteration 2800 of 5000, 0.967s, lp: 101.15 Iteration 2900 of 5000, 0.952s, lp: 102.86 Iteration 3000 of 5000, 0.911s, lp: 90.74 Iteration 3100 of 5000, 0.929s, lp: 89.91 Iteration 3200 of 5000, 0.927s, lp: 102.35 Iteration 3300 of 5000, 0.956s, lp: 106.26 Iteration 3400 of 5000, 0.929s, lp: 94.79 Iteration 3500 of 5000, 0.961s, lp: 100.89 Iteration 3600 of 5000, 0.935s, lp: 93.73 Iteration 3700 of 5000, 0.957s, lp: 99.65 Iteration 3800 of 5000, 0.906s, lp: 94.87 Iteration 3900 of 5000, 0.884s, lp: 97.42 Iteration 4000 of 5000, 0.911s, lp: 104.08 Iteration 4100 of 5000, 0.903s, lp: 103.68 Iteration 4200 of 5000, 0.957s, lp: 93.49 Iteration 4300 of 5000, 0.913s, lp: 101.17 Iteration 4400 of 5000, 0.931s, lp: 94.35 Iteration 4500 of 5000, 0.932s, lp: 97.34 Iteration 4600 of 5000, 0.921s, lp: 100.80 Iteration 4700 of 5000, 0.946s, lp: 97.90 Iteration 4800 of 5000, 0.941s, lp: 100.47 Iteration 4900 of 5000, 0.942s, lp: 98.93 Iteration 5000 of 5000, 0.906s, lp: 99.23 Computation time: 47.0 seconds Read 110000 items Read 110000 items Read 110000 items The max Rhat is: 1.16 Plotting posterior histograms. Loading required package: ggmcmc Loading required package: dplyr Attaching package: ‘dplyr’ The following objects are masked from ‘package:data.table’: between, last The following objects are masked from ‘package:stats’: filter, lag The following objects are masked from ‘package:base’: intersect, setdiff, setequal, union Loading required package: tidyr Loading required package: ggplot2 a b beta_loc_flow beta_loc_snow c cor_flow_elev cor_snow_elev cor_snow_flow scale_elev scale_flow scale_snow shape_elev shape_flow shape_snow Plotting traces. a b beta_loc_flow beta_loc_snow c cor_flow_elev cor_snow_elev cor_snow_flow scale_elev scale_flow scale_snow shape_elev shape_flow shape_snow Read 5000 items Read 5000 items Read 5000 items The DIC is: -162.88 null device 1 null device 1