Nothing
#' Annotate a "nacho" object for outliers
#'
#' Add or update `"is_outlier"` column in the `"nacho"` field of an object from
#' a call to [`load_rcc()`] or [`normalise()`] (`nacho_object$nacho`),
#' using the current quality-control thresholds.
#'
#' @inheritParams normalise
#'
#' @export
#'
#' @return A [[list]] object of class `"nacho"`.
#'
#' @examples
#'
#' data(GSE74821)
#' nacho_object <- check_outliers(GSE74821)
#' head(nacho_object$nacho)
#'
check_outliers <- function(nacho_object) {
if (missing(nacho_object)) {
stop(
'[NACHO] "nacho_object" is missing, results from "load_rcc()" and/or "normalise()" is mandatory!'
)
}
if (!attr(nacho_object, "RCC_type") %in% c("n1", "n8")) {
stop('[NACHO] RCC type must be either "n1" or "n8"!')
}
ot <- nacho_object[["outliers_thresholds"]]
is_house_factor <- "House_factor" %in% colnames(nacho_object[["nacho"]])
if (!is_house_factor) {
nacho_object[["nacho"]][, "House_factor"] <- sum(ot[["House_factor"]]) / 2
}
if (attr(nacho_object, "RCC_type") == "n1") {
nacho_object[["nacho"]][, "is_outlier"] <- {
nacho_object[["nacho"]][, "BD"] < min(ot[["BD"]]) | nacho_object[["nacho"]][, "BD"] > max(ot[["BD"]]) |
nacho_object[["nacho"]][, "FoV"] < ot[["FoV"]] |
nacho_object[["nacho"]][, "PCL"] < ot[["PCL"]] |
nacho_object[["nacho"]][, "LoD"] < ot[["LoD"]] |
nacho_object[["nacho"]][, "Positive_factor"] < min(ot[["Positive_factor"]]) |
nacho_object[["nacho"]][, "Positive_factor"] > max(ot[["Positive_factor"]]) |
nacho_object[["nacho"]][, "House_factor"] < min(ot[["House_factor"]]) |
nacho_object[["nacho"]][, "House_factor"] > max(ot[["House_factor"]])
}
} else {
nacho_object[["nacho"]][, "is_outlier"] <- {
nacho_object[["nacho"]][, "BD"] < min(ot[["BD"]]) | nacho_object[["nacho"]][, "BD"] > max(ot[["BD"]]) |
nacho_object[["nacho"]][, "FoV"] < ot[["FoV"]] |
nacho_object[["nacho"]][, "Positive_factor"] < min(ot[["Positive_factor"]]) |
nacho_object[["nacho"]][, "Positive_factor"] > max(ot[["Positive_factor"]]) |
nacho_object[["nacho"]][, "House_factor"] < min(ot[["House_factor"]]) |
nacho_object[["nacho"]][, "House_factor"] > max(ot[["House_factor"]])
}
}
if (!is_house_factor) {
nacho_object[["nacho"]][, "House_factor"] <- NULL
}
nacho_object
}
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.