cudaBayesreg (version 0.3-13)

post.simul.betadraw: Postprocessing of MCMC simulation

Description

Postprocessing of MCMC simulation. Boxplots of posterior distributions for regressor coefficient beta[vreg] in two cases: estimates for 30 time series of random voxels in active cortex areas; estimates for 30 time series of random voxels in non-active cortex areas.

Usage

post.simul.betadraw(out,vreg = 2)

Arguments

out

List of output objects of MCMC simulation

vreg

regression variable to map; default vreg=2

concept

• MCMC
• plot
• hierarchical model

Details

Post-process analysis

See Also

cudaMultireg.slice, regpostsim, post.ppm, post.tseries

Examples

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")
}
post.simul.betadraw(out=out, vreg=2)
post.simul.betadraw(out=out, vreg=4)