Scale data using a subset of the data

this should reduce the overall variance.

Usage

scale_with_subset(
  data,
  subset,
  config,
  preserveMean = FALSE,
  get_scales = TRUE
)

Arguments

data

the whole dataset

subset

a subset of the dataset

config

configuration

preserveMean

default FALSE - sets mean to zero

get_scales

return a list of transformed data and the scaling parameters

See also

Other preprocessing: INTERNAL_FUNCTIONS_BY_FAMILY, apply_to_response_matrix(), filter_proteins_by_peptide_count(), get_robscales(), normalize_log2_robscale(), robust_scale(), scale_with_subset_by_factors()

Examples

bb <-sim_lfq_data_peptide_config(Nprot = 100)
#> creating sampleName from fileName column
#> completing cases
#> completing cases done
#> setup done
conf <- bb$config$clone(deep=TRUE)
sample_analysis <- bb$data

res <- transform_work_intensity(sample_analysis, conf, log2)
#> Column added : log2_abundance
s1 <- get_robscales(res, conf)
res <- scale_with_subset(res, res, conf)
#> Joining with `by = join_by(sampleName, isotopeLabel, protein_Id, peptide_Id)`
s2 <- get_robscales(res$data, conf)
stopifnot(abs(mean(s1$mads) - mean(s2$mads)) < 1e-6)