meanAbundance: Obtain average expression from replicates

In PYangLab/PhosR: A set of methods and tools for comprehensive analysis of phosphoproteomics data

Obtain average expression from replicates

Description

Obtain average expression from replicates

Usage

meanAbundance(mat, grps)

Arguments

mat	a matrix with rows correspond to phosphosites and columns correspond to samples.
grps	a string specifying the grouping (replciates).

Value

a matrix with mean expression from replicates

Examples

data('phospho_L6_ratio_pe')
data('SPSs')
data('PhosphoSitePlus')

grps = gsub('_.+', '', colnames(phospho.L6.ratio.pe))

# Construct a design matrix by condition
design = model.matrix(~ grps - 1)

# phosphoproteomics data normalisation using RUV
L6.sites = paste(sapply(GeneSymbol(phospho.L6.ratio.pe), function(x)paste(x)),
                 ";",
                 sapply(Residue(phospho.L6.ratio.pe), function(x)paste(x)),
                 sapply(Site(phospho.L6.ratio.pe), function(x)paste(x)),
                 ";", sep = "")
ctl = which(L6.sites %in% SPSs)
phospho.L6.ratio.pe = RUVphospho(phospho.L6.ratio.pe,
                                 M = design, k = 3,ctl = ctl)

phosphoL6 = SummarizedExperiment::assay(phospho.L6.ratio.pe, "normalised")
phosphoL6.mean <- meanAbundance(phosphoL6, grps = grps)

PYangLab/PhosR documentation built on Nov. 1, 2024, 1:47 p.m.