xsample(): An R Function for Sampling Linear Inverse Problems
View Researcher's Other CodesDisclaimer: The provided code links for this paper are external links. Science Nest has no responsibility for the accuracy, legality or content of these links. Also, by downloading this code(s), you agree to comply with the terms of use as set out by the author(s) of the code(s).
Authors | Karel Van den Meersche, Karline Soetaert, Dick Van Oevelen |
Journal/Conference Name | JOURNAL OF STATISTICAL SOFTWARE |
Paper Category | Other |
Paper Abstract | An R function is implemented that uses Markov chain Monte Carlo (MCMC) algorithms to uniformly sample the feasible region of constrained linear problems. Two existing hit-and-run sampling algorithms are implemented, together with a new algorithm where an MCMC step reflects on the inequality constraints. The new algorithm is more robust compared to the hit-and-run methods, at a small cost of increased calculation time. |
Date of publication | 2009 |
Code Programming Language | R |
Comment |