Index: tol/stdlib/math/stat/models/bayesian/bysMcmc/_inference_functions.tol
===================================================================
--- tol/stdlib/math/stat/models/bayesian/bysMcmc/_inference_functions.tol	(revision a4e42555295dda0b27cf1203a05b25d3f5f8ca6f)
+++ tol/stdlib/math/stat/models/bayesian/bysMcmc/_inference_functions.tol	(revision c6d5b9eb509b88828dc44c2e3e33c1982d66ea2b)
@@ -475,6 +475,5 @@
 ////////////////////////////////////////////////////////////////////////////
 {
-  Matrix aux = MatReadRows(path, burnin, length*thinning);
-  SubRow(aux,Range(1,length*thinning,thinning))
+  MatReadRows(path, burnin, length, thinning)
 };
 
@@ -484,6 +483,6 @@
 /////////////////////////////////////////////////////////////////////////////
 {
-  Matrix mcmc.full = MatReadFile(path);
-  Real max = (Rows(mcmc.full)-burnin)/thinning;
+  Set mcmc.dim = MatReadDimensions(path);
+  Real max = (mcmc.dim::Rows-burnin)/thinning;
   Real length = Min(max,length_);
   Matrix mcmc_ = BysMcmc::Inference::loadMcmc(path, burnin, length, thinning);