Learn R Programming

cudaBayesreg (version 0.3-6)

plot.hcoef.post: Plot Method for Hierarchical Model Coefficients

Description

plot.hcoef.post plots arrays of hierarchical coefficients.

Usage

plot.hcoef.post(x,spmname,spm, burnin=trunc(.1*R), ...)

Arguments

x
betadraw object generated by the MCMC simulation
spmname
name associated with the thresholded voxels, e.g. "activated", "non-activated"
spm
threshold active, or non-active voxel, coefficients
burnin
n. of draws to burnin, def: .1*R
...
standard graphics parameters

concept

  • MCMC
  • plot
  • hierarchical model

Details

See description of similar function in bayesm::plot.bayesm.hcoef.

See Also

cudaMultireg.slice, pmeans.hcoef, regpostsim, post.simul.hist, post.simul.betadraw

Examples

Run this code
## load simulation
fsave <- "/tmp/simultest1.sav"
load(file=fsave)
cat("loaded",fsave,"\n")
vreg <- 2
pmeans <- pmeans.hcoef(out$betadraw) 
px <- regpostsim(pmeans, vreg=vreg)
spma <- px$spma # active voxels
spmn <- px$spmn # non-active voxels
plot(out$betadraw,spmname="activated",spm=spma)
plot(out$betadraw,spmname="non-activated",spm=spmn)

Run the code above in your browser using DataLab