bhmaxSubgraph {apComplex} | R Documentation |
Given an adjacency matrix of bait-hit AP-MS protein data, this function finds the maximal BH-complete subgraphs and reports them as an affiliation matrix.
bhmaxSubgraph(adjMat,unrecip=1)
adjMat |
adjMat is an N by (N+M) adjacency matrix with N equal
to the number of bait proteins and M equal to the number of hit-only
proteins. adjMat should have row and column names corresponding to
the proteins in the experiment. An entry of "1" in the ith row and jth
column of adjMat corresponds to bait protein i finding protein j as a
hit. All other entries should be 0. |
unrecip |
By default set to 1 so that unreciprocated bait-bait edges are treated as present. If set to 0, unreciprocated bait-bait edges will be treated as absent. |
A BH-complete subgraph with n bait nodes and m hit-only nodes for AP-MS data is defined as a subgraph for which all n*(n-1)+nm directed edges exist. A maximal BH-complete subgraph is a BH-complete subgraph which is not contained in any other BH-complete subgraph.
This function first forms a symmetric matrix out of the first N columns of
adjMat
. By default, unreciprocated bait-bait observations will be
treated as present. If unrecip
is set to 0, they will be treated as
absent. If the sensitivity of the AP-MS technology is believed to be less
than the specificity, then it is suggested that unrecip
=1.
A matrix with (N+M) rows and a number of columns equal to the number of maximal BH-complete subgraphs. An entry of "1" in the ith row and jth column indicates membership of the ith protein in the jth maximal BH-complete subgraph.
Denise Scholtens
Scholtens D and Gentleman R. Making sense of high-throughput protein-protein interaction data. Statistical Applications in Genetics and Molecular Biology 3, Article 39 (2004).
Scholtens D, Vidal M, and Gentleman R. Local modeling of global interactome networks. Submitted.
data(apEX) PCMG0 <- bhmaxSubgraph(apEX) PCMG1 <- mergeComplexes(PCMG0,apEX,sensitivity=.7,specificity=.75)