Description Usage Arguments Value Author(s) Examples
Lesion to symptom mapping performed on a prepared matrix. Brunner-Munzel tests are performed using each column of the matrix to split the behavioral scores in two groups. This function relies on a compiled version for fast processing.
1 2 3 4 |
lesmat |
binary matrix (0/1) of voxels (columns) and subjects (rows). |
behavior |
vector of behavioral scores. |
permuteNthreshold |
(default=9) Voxels lesioned in fewer than 9 subjects may yield incorrect p-values with Brunner-Munzel tests, so they need to identify pvalues through individualized permutations. This parameter sets the threshold to find which voxels need permutations. See Medina et al (2010). |
alternative |
(default="greater") It is assumed that healthy voxels (0) have greater behavioral scores. If your data follow an inverted relationship choose "less" or "two.sided". |
statOnly |
logical (default=FALSE), skips some computations, mostly for internal use to speed up some things. |
nperm |
(default=1000) Number of permutations to perform on entire volumes when needed for multiple comparisons corrections (i.e., in FWERperm). |
npermBM |
(default=20000) Number of permutations to perform
at every single voxel below |
FWERperm |
logical (default=FALSE) whether to perform permutation based FWER thresholding. |
v |
(default=1) which voxel to record at each permutation with FWERperm. All software use the peak voxel (v=1), but you can choose a voxel further down the list to relax the threshold (i.e., v=10 for 10 highest voxel) (see Mirman (2017)). |
pThreshold |
(default=0.05) what threshold to use for FWER |
permuteAllVoxelsBM |
(default=FALSE) whether to force the
permutation-based p-value calulation for all voxels, instead of
applying only to voxels below |
showInfo |
display info messagges when running the function. |
... |
other arguments received from |
List of objects returned:
statistic
- vector of statistical values
pvalue
- vector of pvalues
zscore
- vector of zscores
perm.vector
- (optional) vector of permuted statistics
perm.FWERthresh
- (optional) permutation threshold established
from the distribution of perm.vector
Dorian Pustina
Note on zscores qnorm gives same values as MRIcron and relies on the normal distribution. however, we are computing t-scores, and should have relied on that distribution, which is the t-score itself.
1 2 3 4 5 6 7 |
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.