R/tidyMS_aggregation.R
medpolish_protein_estimates.Rdmedian polish from normalized peptide intensities
medpolish_protein_estimates(data, config)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()
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)))