Skip to contents

median polish from normalized peptide intensities

Source: R/tidyMS_aggregation.R
medpolish_protein_estimates.Rd

median polish from normalized peptide intensities

Usage

medpolish_protein_estimates(data, config)

See also

Other aggregation: INTERNAL_FUNCTIONS_BY_FAMILY, aggregate_intensity_topN(), estimate_intensity(), intensity_summary_by_hkeys(), medpolish_estimate(), medpolish_estimate_df(), medpolish_estimate_dfconfig(), plot_estimate(), plot_hierarchies_add_quantline(), plot_hierarchies_line(), plot_hierarchies_line_df(), rlm_estimate(), rlm_estimate_dfconfig()

Other deprecated: INTERNAL_FUNCTIONS_BY_FAMILY, intensity_summary_by_hkeys()

Examples


istar <- prolfqua_data("data_ionstar")$normalized()
istar$config <- old2new(istar$config)
istar_data <- istar$data
res <- medpolish_protein_estimates(
  istar_data,
  istar$config
)
#> starting aggregation
#> Warning: medpolish() did not converge in 10 iterations
#> Warning: medpolish() did not converge in 10 iterations

dr <- res("unnest")$data
stopifnot(nrow(dr) ==
  length(unique(istar$data$protein_Id)) * length(unique(istar$data$raw.file)))