For a single sample of symmetric matrices with fixed trace, test eigenvalues of the population mean.
For multiple samples of symmetric matrices with fixed trace, test for equality of the eigenvalues of the population means.
The test statistic is calculated by stat_fixedtrace().
test_fixedtrace(x, evals = NULL, B, maxit = 25)stat_fixedtrace(x, evals = NULL)
A TFORGE object (see boot_calib() or chisq_calib()) with the eigenvalues of the null hypothesis in the null_evals attribute for t0.
A single sample of symmetric matrices or multiple samples of symmetric matrices. See as_flat().
When x is a single sample, the null hypothesis is that the (extrinsic) mean of the population has eigenvalues equal to evals. For multiple samples evals must be omitted.
Number of bootstrap samples. If B = 'chisq' then a chi-squared calibration is used instead.
The maximum number of Newton steps allowed in empirical likelihood optimisation owen:2013TFORGE.
This function uses a form of weighted bootstrapping called b-boostrapping hall1999inTFORGE. An empirical distribution is defined by sampling weights for each observation in the original sample.
The sampling weights must be such that the (extrinsic) mean of the empirical distribution is $$c\hat{Q} \Lambda \hat{Q}^\top,$$
where \(\hat{Q}\) are the eigenvectors of the sample mean, \(\Lambda\) is a diagonal matrix of eigenvalues specified by either the null hypothesis (for single sample tests) or estimated as the common eigenvalues of multiple populations (for k-sample tests).
In some situations \(c\) is a free scalar to enable projection of the Euclidean mean to the extrinsic mean, otherwise \(c=1\).
If no such sampling weights exist (i.e. the convex hull of the data does not contain \(c\hat{Q} \Lambda \hat{Q}^\top\)), then the test rejects with pval=0 and a warning.
The sampling weights are also optimised to maximise empirical likelihood owen:2013TFORGE.
For a single sample the null hypothesis is that the population (extrinsic) mean has eigenvalues of evals; the alternative hypothesis is that the eigenvalues are not equal to evals.
For multiple samples, evals must be omitted and the null hypothesis is that the population (extrinsic) means have the same eigenvalues.
Test hypotheses described below. The fixed trace constraint forces the vector of eigenvalues to lie in a plane. The test statistic accounts for this constraint by using an orthonormal basis in the plane. Weighted bootstrap calibration is used (see 'Weighted Bootstrapping' below).
Eigenvalues must be distinct.
The test statistic is calculated by stat_fixedtrace().