stat_img: Create Statistics Image
In muschellij2/extrantsr: Extra Functions to Build on the 'ANTsR' Package

Creates output from a list of images, performing voxel-wise operations to create a statistic image

stat_img(
  imgs,
  func = c("mean", "median", "mode", "pct", "peak", "sd", "var", "mad", "sum", "prod",
    "z", "quantile", "staple_prob", "staple_label"),
  finite = TRUE,
  masks = NULL,
  na_masks = TRUE,
  ...
)

`imgs`	Character vector, list of characters, or object of class `nifti`
`func`	Function to perform voxel-wise on list of images
`finite`	Should non-finite values be removed?
`masks`	Character vector, list of characters, or object of class `nifti` same length as images.
`na_masks`	if `masks` are given, should values of `0` be turned into `NA`?
`...`	Addictional arguments to pass to `func`

Object of class nifti

When func = "mode", the data is tabulated and then max.col is run. The user can pass ties.method to determine how to break ties. The default is "first" as compared to max.col where it is "random". When func = "peak", density is run on each voxel and the value with the maximum. If the number of unique values is only 1, that value is returned.