ppiPermTest: ppiPermTest

Description Usage Arguments Value Author(s) Examples

Description

Performs a permutation test for enrichment of PPI edges given a database. Peptides are selected from each module and mapped to potential protein parents in the mass tag database. We check if these proteins are found in the PPI network, and record any edges between them. This is compared to edges found using randomly selected proteins (taken from the mass tag database). A p-value is computed as the number of times the randomly sampled proteins incurred more edges than the observed proteins, divided by the number of iterations.

Usage

1
2
ppiPermTest(pnet, pepdat, pepinfo, pepColName, pi_colName, pi_edges,  
   			   threshold, iterations)

Arguments

pnet

procona network object

pepdat

the data matrix with peptides as columns.

pepinfo

Maps peptides to proteins ... same format as in ppiTable

pepColName

The column in pepinfo with peptide IDs... as in pepdat (the peptide data matrix)

pi_colName

The column in pepinfo that maps peptides to unit found in pi_edges

pi_edges

Must be two columns A-B ... sort out evidence levels (in vivo or in vitro) in advance

threshold

Minimum peptide correlation with module eigenvector.

iterations

Number of repititions

Value

returns list of test results.

Author(s)

David L Gibbs

Examples

1
2
3
4
data(ProCoNA_Data)
#net1 <- buildProconaNetwork("peptide network", peptideData, pow=12)
ppis <- data.frame(A=sample(masstagdb$Reference, 50), B=sample(masstagdb$Reference, 50))
ppiPermTest(net1, peptideData, masstagdb, "Mass_Tag_ID", "Reference", ppis, 0.33, 100)

Gibbsdavidl/ProCoNA documentation built on May 8, 2019, 7:51 p.m.