set.seed(123) # Makes example reproducible
# Create example data
data <- create_synthetic_data(
n_proteins = 100,
frac_change = 0.05,
n_replicates = 3,
n_conditions = 2,
method = "effect_random"
)
# Calculate median intensities
qc_median_intensities(
data = data,
sample = sample,
grouping = peptide,
intensity = peptide_intensity_missing,
plot = FALSE
)
# Plot median intensities
qc_median_intensities(
data = data,
sample = sample,
grouping = peptide,
intensity = peptide_intensity_missing,
plot = TRUE
)
Run the code above in your browser using DataLab