Comments (1)
The error seems to happen here:
file: classMcmcSampler.R
function: getSample.mcmcSampler
# Sample size
if(thin == 1 && !is.null(numSamples)){
# This should already be the thinning step?
sel <- seq(1,dim(temp)[1], len = (floor(numSamples/length(sampler$chain))))
temp <- temp[sel,]
}
#############
if (!is.null(whichParameters)) temp = temp[,whichParameters]
out[[i]] = coda::mcmc(temp, start = start, end = end, thin = thin)
I don't really understand what this is supposed to do. Should the thinning be applied such that numSamples determines how many samples collected per chain or from all chains together? (seems like the latter)
The problem that happens here - to my undestanding - is that the chain is already thinned out after "if(thin == 1 && !is.null(numSamples)){ .... }". In this case the coda::mcmc() function with thinning=1 is only there to create the mcmc object from the temp matrix. However, if that is the intention we need to adjust the end parameter of coda::mcmc, whenever the temp matrix is already thinned.
from bayesiantools.
Related Issues (20)
- Pass on extra arguments to the likelihood HOT 4
- Unwritable Outputs For Some Random Seeds HOT 7
- CRAN problem (HTML5)
- Combine DE and DEzs code
- CRAN requires HTML5 for documentation HOT 1
- Add red card dataset to DHARMa?
- Metropolis sampler startup error: Matrix seems negative semi-definite HOT 1
- No longer on CRAN HOT 2
- Incoming CRAN / win builder checks - detritus in temp directory HOT 3
- bridgesample function neccessary? HOT 2
- Questions about BT
- Simulation-based calibration HOT 1
- How to tell if the parallel computation is actually turned on HOT 7
- Number of calls to likelihood function HOT 1
- CRAN DHARMa.Rd issue HOT 3
- Vignette - should there be an example for convergence check?? HOT 1
- Examples in Vignette with eval = F intentional? HOT 1
- Spatially constrained priors? HOT 3
- AR1 function sd estimate wrongly scaled HOT 3
- Writing a hierarchical model - unsure where to assign population-level parameters
Recommend Projects
-
React
A declarative, efficient, and flexible JavaScript library for building user interfaces.
-
Vue.js
🖖 Vue.js is a progressive, incrementally-adoptable JavaScript framework for building UI on the web.
-
Typescript
TypeScript is a superset of JavaScript that compiles to clean JavaScript output.
-
TensorFlow
An Open Source Machine Learning Framework for Everyone
-
Django
The Web framework for perfectionists with deadlines.
-
Laravel
A PHP framework for web artisans
-
D3
Bring data to life with SVG, Canvas and HTML. 📊📈🎉
-
Recommend Topics
-
javascript
JavaScript (JS) is a lightweight interpreted programming language with first-class functions.
-
web
Some thing interesting about web. New door for the world.
-
server
A server is a program made to process requests and deliver data to clients.
-
Machine learning
Machine learning is a way of modeling and interpreting data that allows a piece of software to respond intelligently.
-
Visualization
Some thing interesting about visualization, use data art
-
Game
Some thing interesting about game, make everyone happy.
Recommend Org
-
Facebook
We are working to build community through open source technology. NB: members must have two-factor auth.
-
Microsoft
Open source projects and samples from Microsoft.
-
Google
Google ❤️ Open Source for everyone.
-
Alibaba
Alibaba Open Source for everyone
-
D3
Data-Driven Documents codes.
-
Tencent
China tencent open source team.
from bayesiantools.