function (obj, test, test_type = "wt", which_model = "full",
sig_level = 0.1, point_alpha = 0.2, sig_color = "red", highlight = NULL)
{
stopifnot(is(obj, "sleuth"))
if (test_type == "lrt") {
stop("Currently only works for the Wald test.", " Eventually we will do something for the likelihood ratio test.",
" Suggestions? Email us.")
}
res <- sleuth_results(obj, test, test_type, which_model,
rename_cols = FALSE, show_all = FALSE)
res <- dplyr::mutate(res, significant = qval < sig_level)
p <- ggplot(res, aes(b, -log10(qval)))
p <- p + geom_point(aes(colour = significant), alpha = point_alpha)
p <- p + scale_colour_manual(values = c("black", sig_color))
p <- p + xlab("beta_value")
p <- p + ylab("-log10(qval)")
p <- p + geom_vline(xintercept = 0, colour = "black", linetype = "longdash")
p
}
res <- sleuth_results(obj, test, test_type, which_model,
rename_cols = FALSE, show_all = FALSE)
it's false, it should be : pval_aggregate = FALSE
function (obj, test, test_type = "wt", which_model = "full", sig_level = 0.1, point_alpha = 0.2, sig_color = "red", highlight = NULL) { stopifnot(is(obj, "sleuth")) if (test_type == "lrt") { stop("Currently only works for the Wald test.", " Eventually we will do something for the likelihood ratio test.", " Suggestions? Email us.") } res <- sleuth_results(obj, test, test_type, which_model, rename_cols = FALSE, show_all = FALSE) res <- dplyr::mutate(res, significant = qval < sig_level) p <- ggplot(res, aes(b, -log10(qval))) p <- p + geom_point(aes(colour = significant), alpha = point_alpha) p <- p + scale_colour_manual(values = c("black", sig_color)) p <- p + xlab("beta_value") p <- p + ylab("-log10(qval)") p <- p + geom_vline(xintercept = 0, colour = "black", linetype = "longdash") p }
res <- sleuth_results(obj, test, test_type, which_model, rename_cols = FALSE, show_all = FALSE) it's false, it should be : pval_aggregate = FALSE