HRM
R package providing statistical tests for high-dimensional repeated measures or split-plot designs.
Science Score: 23.0%
This score indicates how likely this project is to be science-related based on various indicators:
-
○CITATION.cff file
-
○codemeta.json file
-
○.zenodo.json file
-
✓DOI references
Found 6 DOI reference(s) in README -
○Academic publication links
-
✓Committers with academic emails
1 of 3 committers (33.3%) from academic institutions -
○Institutional organization owner
-
○JOSS paper metadata
-
○Scientific vocabulary similarity
Low similarity (16.2%) to scientific vocabulary
Keywords
Repository
R package providing statistical tests for high-dimensional repeated measures or split-plot designs.
Statistics
- Stars: 0
- Watchers: 1
- Forks: 0
- Open Issues: 0
- Releases: 0
Topics
Metadata Files
README.md
HRM 1.2.0
R package for analysing high-dimensional repeated measures for factorial designs. A description of this package can be found in [1], theoretical derivations of the test statistics are in [2] and [3].
To install the current development version:
``` r
install devtools package
if (!requireNamespace("devtools", quietly = TRUE)) { install.packages("devtools") }
install package
devtools::install_github("happma/HRM", ref = "test", dependencies = TRUE) library(HRM) ```
With this package it is possible to test for main and interaction effects of up to three whole- or subplot-factors. In total, a maximum of four factors can be used. There are two different S3 methods available. The first method requires a list of matrices in the wide table format. The second method requires a data.frame in the long table format.
``` r
hrm_test with a list of matrices
number patients per group
n = c(10,10)
number of groups
a=2
number of variables
d=40
defining the list consisting of the samples from each group
mu1 = mu2 = rep(0,d)
autoregressive covariance matrix
sigma1 = diag(d) for(k in 1:d) for(l in 1:d) sigma1[k,l] = 1/(1-0.5^2)0.5^(abs(k-l)) sigma_2 = 1.5sigma1 X = list(mvrnorm(n[1],mu1, sigma1), mvrnorm(n[2],mu2, sigma_2)) X=lapply(X, as.matrix)
hrm_test(data=X, alpha=0.05)
hrm.test with a data.frame using a 'formula' object
using the EEG dataset
hrm_test(value ~ groupregionvariable, subject = "subject", data = EEG) ```
To get confidence intervals for each factor combination you can use the generic function 'confint' for an object of class 'HRM'. This function calculates simultaneous confidence intervals which maintains the family wise error rate (FWER). See the following code:
``` r
using the EEG dataset
z <- hrm_test(value ~ groupregionvariable, subject = "subject", data = EEG)
calculate 99% confidence intervals
confint(z, level = 0.99)
```
In the data there are 4 variables with each 10 regions. We can use a multivariate approach as the variables are on different scales. For that, we can use the function 'hrm_test' with the argument 'variable' set to the column name which contains the factor variable for the variables.
``` r
using the EEG dataset
hrm_test(value ~ group*region, subject = subject, variable = variable, data = EEG) ```
Additionally, the package can be used with a GUI.
r
hrm_GUI()
References
[1] Happ, M., Harrar, S. W., and Bathke, A. C. (2018). HRM: An R Package for Analysing High-dimensional Multi-factor Repeated Measures. The R Journal 10(1), 534--548. https://journal.r-project.org/archive/2018/RJ-2018-032/index.html
[2] Happ, M., Harrar S. W. and Bathke, A. C. (2017). High-dimensional Repeated Measures. Journal of Statistical Theory and Practice. 11(3), 468-477. URL: doi:10.1080/15598608.2017.1307792.
[3] Happ, M., Harrar, S. W., & Bathke, A. C. (2016). Inference for low‐and high‐dimensional multigroup repeated measures designs with unequal covariance matrices. Biometrical Journal, 58(4), 810-830. doi:10.1002/bimj.201500064
Owner
- Name: happma
- Login: happma
- Kind: user
- Repositories: 6
- Profile: https://github.com/happma
PHD in statistics | data scientist | actuary
GitHub Events
Total
Last Year
Committers
Last synced: over 2 years ago
Top Committers
| Name | Commits | |
|---|---|---|
| happma | m****p@a****t | 30 |
| happma | b****1@i****t | 28 |
| Martin Happ | m****p@a****t | 8 |
Committer Domains (Top 20 + Academic)
Issues and Pull Requests
Last synced: 12 months ago
All Time
- Total issues: 0
- Total pull requests: 0
- Average time to close issues: N/A
- Average time to close pull requests: N/A
- Total issue authors: 0
- Total pull request authors: 0
- Average comments per issue: 0
- Average comments per pull request: 0
- Merged pull requests: 0
- Bot issues: 0
- Bot pull requests: 0
Past Year
- Issues: 0
- Pull requests: 0
- Average time to close issues: N/A
- Average time to close pull requests: N/A
- Issue authors: 0
- Pull request authors: 0
- Average comments per issue: 0
- Average comments per pull request: 0
- Merged pull requests: 0
- Bot issues: 0
- Bot pull requests: 0
Top Authors
Issue Authors
Pull Request Authors
Top Labels
Issue Labels
Pull Request Labels
Packages
- Total packages: 1
-
Total downloads:
- cran 300 last-month
- Total dependent packages: 0
- Total dependent repositories: 0
- Total versions: 12
- Total maintainers: 1
cran.r-project.org: HRM
High-Dimensional Repeated Measures
- Homepage: http://github.com/happma/HRM
- Documentation: http://cran.r-project.org/web/packages/HRM/HRM.pdf
- License: GPL-2 | GPL-3
- Status: removed
-
Latest release: 1.2.1
published about 6 years ago
Rankings
Maintainers (1)
Dependencies
- MASS * depends
- R >= 3.4.0 depends
- ggplot2 * depends
- matrixcalc * depends
- plyr * depends
- Rcpp >= 0.12.16 imports
- data.table * imports
- doBy * imports
- mvtnorm * imports
- pseudorank >= 0.3.7 imports
- reshape2 * imports
- tcltk * imports
- xtable * imports
- RGtk2 >= 2.8.0 suggests
- cairoDevice * suggests
- testthat * suggests