inst/doc/study-wtccc-mice.R

In mvMAPIT: Multivariate Genome Wide Marginal Epistasis Test

## ----setup, include = FALSE---------------------------------------------------
knitr::opts_chunk$set(
  collapse = TRUE,
  comment = "#",
  fig.width=7,
  fig.height=5
)
library(knitr)

## ----load_dependencies, message=FALSE-----------------------------------------
library(mvMAPIT)
library(ggplot2)
library(dplyr)

## ----run_mvmapit, eval = FALSE------------------------------------------------
#  mvmapit_TRAIT <- mvmapit(
#    t(TRAIT$genotype),
#    t(TRAIT$phenotype),
#    test = "hybrid"
#  )

## ----load_data, eval = FALSE--------------------------------------------------
#  mice_SI_paper <- readRDS("mice_SI_paper.rds")
#  mice_HCTHGB_MCVMCH <- readRDS("mice_HCTHGB_MCVMCH.rds")

## ----all_traits, eval = FALSE-------------------------------------------------
#  for_ticks_chr <- aggregate(position ~ chr, mice_data$fisher, function(x) c(first = min(x), last = max(x))) %>%
#    mutate(tick = floor((position[,"first"] + position[,"last"]) / 2)) %>%
#    mutate(chr2 = case_when(chr %% 5 == 0 ~ as.character(chr),
#                            chr == 1 ~ as.character(chr),
#                            TRUE ~ ""))
#  for_facetgrid_row <- as_labeller(c(`1` = "Trait #1", `2` = "Trait #2", `3` = "Covariance", `4` = "Combined"))
#  gg <- mice_SI_paper$fisher %>% ggplot(aes(
#        x = position,
#        y = -log10(pplot),
#        color = factor(color)
#      )) +
#        geom_point_rast(
#          size = 0.7) +
#        scale_color_manual(
#          values = c("#8b8b8b", "#bfbfbf", "#1b9e77")
#        ) +
#        scale_y_continuous(breaks = c(0, 5, 10),
#                           labels = c("0", "5", ">10")) +
#        geom_hline(
#          aes(
#            yintercept = -log10(5.179737e-06),
#            linetype = "Bonferroni"
#          ),
#          color = "#d95f02",
#          size = 0.3
#        ) +
#        theme_bw() +
#        facet_grid(x ~ y) +
#        theme(
#          panel.grid.major.x = element_blank(),
#          legend.position = "bottom",
#          text = element_text(family = "Times"),
#        ) +
#        labs(
#          y = bquote(-log[10](p)),
#          color = "") +
#        scale_x_continuous("Chromosome",
#                           breaks = for_ticks_chr$tick,
#                           labels = for_ticks_chr$chr2) +
#        scale_linetype_manual(name = "", values = c('dashed'))
#  show(gg)

## ----mice_data, eval = FALSE--------------------------------------------------
#  gg <- mice_HCTHGB_MCVMCH$fisher %>% ggplot(aes(
#        x = position,
#        y = -log10(pplot),
#        color = factor(color)
#      )) +
#        geom_point_rast(
#          size = 0.7) +
#        scale_color_manual(
#          values = c("#8b8b8b", "#bfbfbf", "#1b9e77")
#        ) +
#        scale_y_continuous(breaks = c(0, 5, 10),
#                           labels = c("0", "5", ">10")) +
#        geom_hline(
#          aes(
#            yintercept = -log10(5.179737e-06),
#            linetype = "Bonferroni"
#          ),
#          color = "#d95f02",
#          size = 0.3
#        ) +
#        theme_bw() +
#        facet_grid(row ~ case, labeller = labeller(row = for_facetgrid_row)) +
#        theme(
#          panel.grid.major.x = element_blank(),
#          legend.position = "bottom",
#          text = element_text(family = "Times"),
#        ) +
#        labs(
#          y = bquote(-log[10](p)),
#          color = "") +
#        scale_x_continuous("Chromosome",
#                           breaks = for_ticks_chr$tick,
#                           labels = for_ticks_chr$chr2) +
#        scale_linetype_manual(name = "", values = c('dashed'))
#  show(gg)