mincPlotPeak: Plotting of peaks

In Mouse-Imaging-Centre/RMINC: Statistical Tools for Medical Imaging NetCDF (MINC) Files

Plotting of peaks

Description

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

Usage

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)

## End(Not run)

Mouse-Imaging-Centre/RMINC documentation built on Nov. 12, 2022, 1:50 p.m.