cudaBayesreg (version 0.2-1)

post.tseries: Show fitted time series of active voxel

Description

post.tseries plots the fitted time series of a voxel estimated as active.

Usage

post.tseries(out, slicedata, vreg=2)

Details

Use the estimated regression coefficients to vizualize the fitted time series in an active voxel.

References

Adelino R. Ferreira da Silva "A Bayesian Multilevel Model for fMRI Data Analysis", to be published.

See Also

cudaMultireg.slice post.simul.hist post.simul.betadraw post.ppm

Examples

slicedata <- read.fmrisample(slice=3)
fsave <- "/tmp/simultest.sav"
answ <- readline("Run MCMC simulation first ? ")
run <- FALSE
if (substr(answ, 1, 1) == "y") { run <- TRUE }
if(run) {
	out <- cudaMultireg.slice(slicedata, fsave=fsave )
} else {
	load(file=fsave)
	cat("loaded",fsave,"\n")
}
post.tseries(out=out, slicedata=slicedata)