Skip to contents

Title Find Z-score Hit

Usage

find_zscore_hit(table_treate_vs_control, number_barcode = 6, metric = "median")

Arguments

table_treate_vs_control

table computed with the function compute_data_table

number_barcode

Number of barcode that as to be differentially expressed (DE)in order to consider the gene associated DE. Example a gene is associated with 10 shRNA we consider a gene DE if it has at least number_barcode = 5 shRNA DE.

metric

A string containing the metric to use. The value allowed are "median" or "mean".

Value

Return a tibble containing the hit for the Z-score

Examples

object <- get0("object", envir = asNamespace("ScreenR"))
table <- compute_metrics(object,
    control = "TRT", treatment = "Time3",
    day = "Time3"
)

# For the the median
result <- find_zscore_hit(table, number_barcode = 6)
head(result)
#> # A tibble: 6 × 2
#>   Gene     numberOfBarcode
#>   <chr>              <int>
#> 1 Gene_102               8
#> 2 Gene_115               7
#> 3 Gene_117               7
#> 4 Gene_124               7
#> 5 Gene_129               7
#> 6 Gene_13                7

# For the mean
result <- find_zscore_hit(table, number_barcode = 6, metric = "mean")
#> Warning: Returning more (or less) than 1 row per `summarise()` group was deprecated in
#> dplyr 1.1.0.
#>  Please use `reframe()` instead.
#>  When switching from `summarise()` to `reframe()`, remember that `reframe()`
#>   always returns an ungrouped data frame and adjust accordingly.
#>  The deprecated feature was likely used in the ScreenR package.
#>   Please report the issue at <https://github.com/EmanuelSoda/ScreenR/issues>.
head(result)
#> # A tibble: 6 × 2
#>   Gene     numberOfBarcode
#>   <chr>              <int>
#> 1 Gene_102               8
#> 2 Gene_102               8
#> 3 Gene_102               8
#> 4 Gene_102               8
#> 5 Gene_102               8
#> 6 Gene_102               8