post.simul.betadraw: Postprocessing of MCMC simulation

Description Usage Arguments Details See Also Examples

View source: R/post.simul.R

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

1
	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

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
## 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.