Nothing
## ---- echo = FALSE, message = FALSE-------------------------------------------
knitr::opts_chunk$set(fig.width = 10, fig.height = 5,
echo = FALSE, collapse = TRUE, comment = "#>")
## ---- eval = TRUE, echo = TRUE------------------------------------------------
library("seagull")
data("seagull_data")
## ---- eval = TRUE, echo = TRUE------------------------------------------------
fit_sgl1 <- seagull(y = phenotypes[, 1], Z = genotypes, groups = groups)
## ---- eval = TRUE, echo = TRUE------------------------------------------------
attributes(fit_sgl1)
## ---- eval = TRUE, echo = TRUE------------------------------------------------
last_solution <- fit_sgl1$loops_lambda
plot(x = seq(1, dim(genotypes)[2], 1),
y = fit_sgl1$random_effects[last_solution,],
xlab = "position", ylab = "effect estimate",
col = "gray80", pch = 16)
## ---- eval = TRUE, echo = FALSE-----------------------------------------------
print(c("The number of ZEROS in the last solution is: ", length(which(fit_sgl1$random_effects[last_solution,]==0))))
## ---- eval = TRUE, echo = 3---------------------------------------------------
last_solution <- fit_sgl1$loops_lambda
plot(x = seq(1, dim(genotypes)[2], 1),
y = fit_sgl1$random_effects[last_solution,],
xlab = "position", ylab = "effect estimate",
col = "gray80", pch = 16)
points(x = seq(1, dim(genotypes)[2], 1),
y = fit_sgl1$random_effects[20,],
pch = 16)
## ---- eval = TRUE, echo = FALSE-----------------------------------------------
print(c("The number of ZEROS in the solution in line 20 is: ", length(which(fit_sgl1$random_effects[20,]==0))))
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.