nr_B_in_A: Compute nr of B per A

In wolski/prolfqua: Proteomics Label Free Quantification

Compute nr of B per A

Description

Compute nr of B per A

Usage

nr_B_in_A(pdata, config, merge = TRUE)

Arguments

pdata	data.frame
config	AnalysisConfiguration

Examples

bb <- prolfqua::prolfqua_data('data_ionstar')$filtered()
bb$config <- old2new(bb$config)
stopifnot(nrow(bb$data) == 25780)
config <- bb$config$clone(deep=TRUE)
data <- bb$data |> dplyr::select(-all_of("nr_peptide_Id_IN_protein_Id"))
hierarchy <- config$table$hierarchy_keys()
res <- nr_B_in_A(data, config)

res$data |>
  dplyr::select_at(c(config$table$hierarchy_keys_depth(),  res$name)) |>
  dplyr::distinct() |>
  dplyr::pull() |> table()


bb <- prolfqua::prolfqua_data('data_skylineSRM_HL_A')
config <- old2new(bb$config_f())
data <- bb$data
data$Area[data$Area == 0] <- NA
analysis <- setup_analysis(data, config)

resDataStart <- bb$analysis(bb$data, config)


nr_B_in_A(resDataStart, config)
nr_B_in_A(resDataStart, config, merge = FALSE)
config$table$hierarchyDepth <- 2
nr_B_in_A(resDataStart, config, merge = FALSE)

bb <- prolfqua::prolfqua_data('data_IonstarProtein_subsetNorm')
bb$config <- old2new(bb$config$clone(deep=TRUE))
nr_B_in_A(bb$data, bb$config)
#undebug(nr_B_in_A)

wolski/prolfqua documentation built on May 12, 2024, 10:16 p.m.