In mcanouil/MiSTr: Mixed effects Score Test for continuous outcomes

knitr::opts_chunk$set(
  collapse = TRUE,
  comment = "#>"
)
options(width = 140)

Mixed effects Score Test


Test for association between a set of SNPS/genes and continuous or binary outcomes by including variant characteristic information and using (weighted) score statistics.

Note:

• From: https://cran.r-project.org/src/contrib/MiST_1.0.tar.gz
• Reference: https://doi.org/10.1002/gepi.21717

Installation

#| eval: false
# Install MiSTr from CRAN:
install.packages("MiSTr")

# Or the the development version from GitHub:
# install.packages("remotes")
remotes::install_github("mcanouil/MiSTr")

MiSTr in Action

library(MiSTr)
data(mist_data)
attach(mist_data)

Continuous Outcome

res <- mist(
  y = phenotypes[, "y_taupi"],
  X = phenotypes[, paste0("x_cov", 0:2)],
  G = genotypes,
  Z = variants_info[, 1, drop = FALSE]
)
str(res)
print(res)

Binary Outcome

res <- mist(
  y = phenotypes[, "y_binary"],
  X = phenotypes[, paste0("x_cov", 0:2)],
  G = genotypes,
  Z = variants_info[, 1, drop = FALSE]
)
str(res)
print(res)

mcanouil/MiSTr documentation built on April 7, 2023, 10:48 a.m.