matie computes $A$ by estimating a generalized $R^2$,
which is computed from the ratio of the likelihood
of an alternative model (allowing dependence between variables)
over the likelihood of a null model (that forces the variables to be independent).
See http://arxiv.org/abs/1303.1828 for details.
The variables should be continuous. $A$ is calculated from the ranks of the values, and ties are currently broken randomly. Ordinal data can be handled, but the properties of $A$ have not been thoroughly investigated for non-continuous data with many ties. $A$ is not currently implemented for category-valued data.
Version 1.2 provides underflow protection to fix a bug that appeared at high dimensions in earlier versions.
| Package: |
| matie |
| Type: |
| Package |
| Version: |
| 1.2 |
| Date: |
| 2013-10-11 |
| License: |
| GPL-3 |
d <- sbd(NULL,n=400,Rsq=0.9)
aScore <- ma(d)
pValue <- ma.test(d,aScore)
Run the code above in your browser using DataLab