Nothing
library(qtl2pleio)
library(testthat)
context("Simulation of one phenotype vector")
X1 <- as.matrix(rbinom(n = 100, size = 1, prob = 1 / 2))
X2 <- as.matrix(rbinom(n = 100, size = 1, prob = 1 / 2))
gemma2::stagger_mats(X1, X2) -> X
Sigma <- calc_Sigma(Vg = diag(1, 2), Ve = diag(1, 2), kinship = diag(100))
test_that("sim1 outputs numeric vector", {
expect_true(is.numeric(sim1(X = X, B = c(1, 2), Sigma = Sigma)))
})
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.