orderGoAnnotations: Orders annotation information

Description Usage Arguments Details Value Author(s) See Also

View source: R/goannotations.R


For a given matrix of annotation information, this function returns the information ordered according to the best fit with the data.


orderGoAnnotations(object, fcol = "GOAnnotations", k = 1:5, n = 5,
  p = 1/3, verbose = TRUE, seed)



An instance of class MSnSet.


The name of the annotations matrix. Default is GOAnnotations.


The number of clusters to test. Default is k = 1:5


The minimum number of proteins per component cluster.


The normalisation factor, per k tested


A logical indicating if a progress bar should be displayed. Default is TRUE.


An optional random number generation seed.


As there are typically many protein/annotation sets that may fit the data we order protein sets by best fit i.e. cluster tightness, by computing the mean normalised Euclidean distance for all instances per protein set.

For each protein set i.e. proteins that have been labelled with a specified term/information criteria, we find the best k cluster components for the set (the default is to testk = 1:5) according to the minimum mean normalised pairwise Euclidean distance over all component clusters. (Note: when testing k if any components are found to have less than n proteins these components are not included and k is reduced by 1).

Each component cluster is normalised by N^p (where N is the total number of proteins per component, and p is the power). Hueristally, p = 1/3 and normalising by N^1/3 has been found the optimum normalisation factor.

Candidates in the matrix are ordered according to lowest mean normalised pairwise Euclidean distance as we expect high density, tight clusters to have the smallest mean normalised distance.

This function is a wrapper for running clustDist, getNormDist, see the "Annotating spatial proteomics data" vignette for more details.


An updated MSnSet containing the newly ordered fcol matrix.


Lisa M Breckels

See Also

addGoAnnotations and example therein.

pRoloc documentation built on Nov. 8, 2020, 6:26 p.m.