Description Usage Arguments Details Value See Also Examples
pmeans.hcoef
processes the MCMC simulation to evaluate the posterior mean of the regression variables.
1 | pmeans.hcoef(x,burnin=trunc(.1*R))
|
x |
betadraw object generated by the MCMC simulation |
burnin |
n. of draws to burnin, def: .1*R |
Post-process MCMC simulation
pmeans |
Posterior Means of Coefficients |
cudaMultireg.slice
,
read.fmrislice
,
regpostsim
,
post.simul.betadraw
,
post.simul.hist
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 | ## Not run:
slicedata <- read.fmrislice(fbase="fmri", slice=3, swap=FALSE)
ymaskdata <- premask(slicedata)
fsave <- paste(tempdir(),"/simultest1",fileext = ".sav", sep="")
answ <- readline("Run MCMC simulation first ? ")
run <- FALSE
if (substr(answ, 1, 1) == "y") { run <- TRUE }
if(run) {
out <- cudaMultireg.slice(slicedata, ymaskdata, R=2000, keep=5, nu.e=3,
fsave=fsave, zprior=FALSE)
} else {
load(file=fsave)
cat("loaded",fsave,"\n")
}
pmeans <- pmeans.hcoef(out$betadraw)
## End(Not run)
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.