Plots a slice containing a peak. Optionally plots a graph of that peak alongside.

mincPlotPeak(peak, anatomy, statistics, dim = 2, crossCol = "green",
  crossSize = 4, plotFunction = NULL, ...)

Arguments

peak

a row from mincFindPeaks

anatomy

a mincArray for the underlying anatomy

statistics

a mincArray for the stats volume

dim

the dimension (1:3)

crossCol

the colour for the cross-hair

crossSize

the size (cex) of the cross-hair

plotFunction

a function which will produce a graph

...

other details to pass on to mincPlotAnatAndStatsSlice

Examples

# NOT RUN {
peaks <- mincFindPeaks(-log10(qvs), "Neonatal:time.to.sac", "pos",
posThreshold=1.3, minDistance=1)
p <- function(peak) {
  gfTiming$voxel <- mincGetWorldVoxel(gfTiming$reljacobians02,
                                      peak["x"], peak["y"], peak["z"])
  qplot(time.to.sac, exp(voxel), data=gfTiming, colour=Neonatal,
        geom="boxplot") + theme_classic()
}
mincPlotPeak(peaks[1,], anatVol, -log10(mincArray(qvs, "Neonatal:time.to.sac")),
             anatLow=700, anatHigh=1400, low=1, high=4, col=heat.colors(244),
             crossCol = "blue", crossSize = 3, plotFunction = p)
# }