Package: ggmix 0.0.2

ggmix: Variable Selection in Linear Mixed Models for SNP Data

Fit penalized multivariable linear mixed models with a single random effect to control for population structure in genetic association studies. The goal is to simultaneously fit many genetic variants at the same time, in order to select markers that are independently associated with the response. Can also handle prior annotation information, for example, rare variants, in the form of variable weights. For more information, see the website below and the accompanying paper: Bhatnagar et al., "Simultaneous SNP selection and adjustment for population structure in high dimensional prediction models", 2020, <doi:10.1371/journal.pgen.1008766>.

Authors:Sahir Bhatnagar [aut, cre], Karim Oualkacha [aut], Yi Yang [aut], Celia Greenwood [aut]

ggmix_0.0.2.tar.gz
ggmix_0.0.2.zip(r-4.7)ggmix_0.0.2.zip(r-4.6)ggmix_0.0.2.zip(r-4.5)
ggmix_0.0.2.tgz(r-4.6-any)ggmix_0.0.2.tgz(r-4.5-any)
ggmix_0.0.2.tar.gz(r-4.7-any)ggmix_0.0.2.tar.gz(r-4.6-any)
ggmix_0.0.2.tgz(r-4.6-emscripten)
manual.pdf |manual.html
DESCRIPTION |NEWS
card.svg |card.png
ggmix/json (API)

# Install 'ggmix' in R:
install.packages('ggmix', repos = c('https://sahirbhatnagar.r-universe.dev', 'https://cloud.r-project.org'))

Bug tracker:https://github.com/sahirbhatnagar/ggmix/issues

Datasets:
  • admixed - Simulated Dataset with 1D Geography
  • karim - Karim's Simulated Data

On CRAN:

Conda:

5.64 score 12 stars 24 scripts 277 downloads 1 mentions 5 exports 11 dependencies

Last updated from:eb1c8a71dd. Checks:9 OK. Indexed: yes.

TargetResultTimeFilesSyslog
linux-devel-x86_64OK137
source / vignettesOK197
linux-release-x86_64OK135
macos-release-arm64OK151
macos-oldrel-arm64OK144
windows-develOK126
windows-releaseOK95
windows-oldrelOK87
wasm-releaseOK108

Exports:gen_structured_modelggmixgicrandom.effectsranef

Dependencies:codetoolsforeachglmnetiteratorslatticeMASSMatrixRcppRcppEigenshapesurvival

Introduction to the ggmix package
Model | Installation | Quick Start | Fit the linear mixed model with Lasso Penalty | Find the Optimal Value of the Tuning Parameter | Get Coefficients Corresponding to Optimal Model | Diagnostic Plots | Observed vs. Predicted Response | QQ-plots for Residuals and Random Effects | Tukey-Anscombe Plot

Last update: 2020-03-09
Started: 2018-06-28

Alternative Inputs for Population Structure
The kinship argument | The U and D arguments | The K argument

Last update: 2020-03-09
Started: 2018-06-28

Incorporating Prior Annotation Weights
Setting the Penalty Factor Argument

Last update: 2020-03-09
Started: 2018-06-28

Readme and manuals

Help Manual

Help pageTopics
Simulated Dataset with 1D Geographyadmixed
Simulation Scenario from Bhatnagar et al. (2018+) ggmix papergen_structured_model
Fit Linear Mixed Model with Lasso or Group Lasso Regularizationggmix
Constructor functions for the different ggmix objectsggmix_data_object new_fullrank_K new_fullrank_kinship new_fullrank_UD new_lowrank_K new_lowrank_kinship new_lowrank_UD
Generalised Information Criteriongic gic.default gic.ggmix_fit
Functions related to eta parameter used in optim and kkt checksfn_eta_lasso_fullrank gr_eta_lasso_fullrank
Karim's Simulated Datakarim
Check of KKT Conditions for Linear Mixed Modelgrr_beta0 grr_sigma2 kkt_check
Estimation of Lambda Sequence for Linear Mixed Model with Lasso Penaltylambdalasso lambdalasso.default lambdalasso.fullrank
Estimation of Linear Mixed Model with Lasso Penaltylmmlasso lmmlasso.default lmmlasso.fullrank
Estimation of Log-likelihood for Linear Mixed Model with Lasso Penaltylogliklasso logliklasso.default logliklasso.fullrank
Plot Method for 'ggmix_fit' objectplot.ggmix_fit plotCoef
Plot the Generalised Information Criteria curve produced by 'gic'plot.ggmix_gic plotGIC
Make predictions from a 'ggmix_fit' objectcoef.ggmix_fit predict.ggmix_fit
Make predictions from a 'ggmix_gic' objectcoef.ggmix_gic predict.ggmix_gic
Print Method for Objects of Class 'ggmix_fit'print.ggmix_fit print.ggmix_gic
Extract Random Effectsrandom.effects random.effects.default ranef ranef.default ranef.ggmix_gic
Estimation of Sigma2 for Linear Mixed Model with Lasso Penaltysigma2lasso sigma2lasso.default sigma2lasso.fullrank