how_to_continue_gibbs_sampling_from_an_accidental_stop
Differences
This shows you the differences between two versions of the page.
Both sides previous revisionPrevious revisionNext revision | Previous revisionLast revisionBoth sides next revision | ||
how_to_continue_gibbs_sampling_from_an_accidental_stop [2018/05/19 21:57] – yutaka | how_to_continue_gibbs_sampling_from_an_accidental_stop [2018/05/19 22:05] – yutaka | ||
---|---|---|---|
Line 1: | Line 1: | ||
====== Recovering from an accidental stop in Gibbs sampling ====== | ====== Recovering from an accidental stop in Gibbs sampling ====== | ||
- | ===== Quick answer | + | ===== Overview ===== |
+ | |||
+ | ==== Quick answer ==== | ||
Our Gibbs sampling programs can continue the sampling only if you have prepared something for the accidental stop before running the program. | Our Gibbs sampling programs can continue the sampling only if you have prepared something for the accidental stop before running the program. | ||
Line 11: | Line 13: | ||
* Otherwise, use a shell script shown in below. | * Otherwise, use a shell script shown in below. | ||
- | ===== Details | + | ==== Details ==== |
The idea is to save the current samples of " | The idea is to save the current samples of " | ||
Line 23: | Line 25: | ||
===== Methods ===== | ===== Methods ===== | ||
+ | |||
+ | In both methods, the continued sampling may give you slightly different results from the non-breaking run. | ||
+ | This is because of numerical error. | ||
+ | The saved samples of variance components have limited precision and it creates errors in subsequent sampling. | ||
+ | In general, this error is trivial and you can ignore it. | ||
==== Option for THRGIBBS1F90 ==== | ==== Option for THRGIBBS1F90 ==== | ||
Line 56: | Line 63: | ||
# Put your favorite value here. | # Put your favorite value here. | ||
- | total_samples=50 | + | total_samples=100000 |
# Save the intermediate samples every this rounds. | # Save the intermediate samples every this rounds. | ||
- | save_rounds=10 | + | save_rounds=10000 |
# set your parameter file | # set your parameter file | ||
- | parameter_file=exmr99s | + | parameter_file=renf90.par |
# thinning: samples saved in every n rounds | # thinning: samples saved in every n rounds |
how_to_continue_gibbs_sampling_from_an_accidental_stop.txt · Last modified: 2024/03/25 18:22 by 127.0.0.1