RUVphospho: RUV for phosphoproteomics data normalisation

Description Usage Arguments Value Examples

View source: R/RUVs.R

Description

This is a wrapper implementation of RUVIII for phosphoproteomics data normalisation. This function will call tailImpute function to impute all the missing values (if there is any) in the phosphoproteomics data for applying RUVIII. It will then return the normalised values for quantified phosphosites and remove imputed values.

Usage

1
RUVphospho(mat, M, ctl, k = NULL, m = 1.6, s = 0.6, keepImpute = FALSE, ...)

Arguments

mat

a matrix with rows correspond to phosphosites and columns correspond to samples.

M

is the design matrix as defined in RUVIII.

ctl

is the stable phosphosites (or negative controls as defined in RUVIII).

k

is the number of unwanted factors as defined in RUVIII.

m

a numeric number for controlling mean downshifting.

s

a numeric number for controlling standard deviation of downshifted sampling values.

keepImpute

a boolean to keep the missing value in the returned matrix.

...

additional parameters that may be passed to RUVIII.

Value

A normalised matrix.

Examples

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
19
20
21
data('phospho_L6_ratio')
data('SPSs')

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

# Cleaning phosphosite label
phospho.site.names = rownames(phospho.L6.ratio)
L6.sites = gsub(' ', '', sapply(strsplit(rownames(phospho.L6.ratio), '~'),
                                function(x){paste(toupper(x[2]), x[3], '',
                                                sep=';')}))
phospho.L6.ratio = t(sapply(split(data.frame(phospho.L6.ratio), L6.sites),
                            colMeans))
phospho.site.names = split(phospho.site.names, L6.sites)

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

# phosphoproteomics data normalisation using RUV
ctl = which(rownames(phospho.L6.ratio) %in% SPSs)
phospho.L6.ratio.RUV = RUVphospho(phospho.L6.ratio, M = design, k = 3,
                                ctl = ctl)

PhosR documentation built on Nov. 8, 2020, 6:54 p.m.