cudaBayesreg (version 0.3-9)

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, ymaskdata, vreg=2)

Arguments

out

output of MCMC simulation

slicedata

list(slice=slice, niislicets=niislicets, mask=mask, dsgn=dsgn); input slicedata used in simulation as returned by read.fmrislice

ymaskdata

list(yn = yn, kin = kin, nreg = nreg); masked and standardised slice data as returned by premask

vreg

number of the active variable to visualize; default(vreg=2).

Details

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

References

Adelino Ferreira da Silva, A Bayesian Multilevel Model for fMRI Data Analysis, Computer Methods and Programs in Biomedicine, to be published.

See Also

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

Examples

## read data and load MCMC simulation 
slicedata <- read.fmrislice(fbase="fmri", slice=3, swap=TRUE)
ymaskdata <- premask(slicedata)
fsave <- "/tmp/simultest1.sav"
load(fsave)
post.tseries(out=out, slicedata=slicedata, ymaskdata=ymaskdata, vreg=2)