post.simul.betadraw: Postprocessing of MCMC simulation

In cudaBayesreg: CUDA Parallel Implementation of a Bayesian Multilevel Model for fMRI Data Analysis

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

Details

Post-process analysis

See Also

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

Examples

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

## End(Not run)

cudaBayesreg documentation built on May 29, 2017, 6:19 p.m.