parGOSim: Protein/DNA Sequence Similarity Calculation based on Gene...

Description Usage Arguments Details Value Author(s) See Also Examples

Description

Protein/DNA Sequence Similarity Calculation based on Gene Ontology (GO) Similarity

Usage

1
2
parGOSim(golist, type = c("go", "gene"), ont = "MF",
  organism = "human", measure = "Resnik", combine = "BMA")

Arguments

golist

A character vector, each component contains a character vector of GO terms or one Entrez Gene ID.

type

Input type of golist, 'go' for GO Terms, 'gene' for gene ID.

semData

GOSemSimDATA object

measure

Default is 'Resnik', could be one of 'Resnik', 'Lin', 'Rel', 'Jiang' or 'Wang'.

combine

Default is 'BMA', could be one of 'max', 'average', 'rcmax' or 'BMA' for combining semantic similarity scores of multiple GO terms associated with protein.

Details

This function calculates protein/DNA sequence similarity based on Gene Ontology (GO) similarity.

Value

A n x n similarity matrix.

Author(s)

Min-feng Zhu <wind2zhu@163.com>, Nan Xiao <http://nanx.me>

See Also

See twoGOSim for calculating the GO semantic similarity between two groups of GO terms or two Entrez gene IDs. See parSeqSim for paralleled protein/DNA similarity calculation based on Smith-Waterman local alignment.

Examples

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
19
20
# Be careful when testing this since it involves GO similarity computation
# and might produce unpredictable results in some environments

require(GOSemSim)
require(org.Hs.eg.db)

# by GO Terms
go1 = c('GO:0005215', 'GO:0005488', 'GO:0005515', 'GO:0005625', 'GO:0005802', 'GO:0005905')  # AP4B1
go2 = c('GO:0005515', 'GO:0005634', 'GO:0005681', 'GO:0008380', 'GO:0031202')  # BCAS2
go3 = c('GO:0003735', 'GO:0005622', 'GO:0005840', 'GO:0006412')  # PDE4DIP
glist = list(go1, go2, go3)
d = godata('org.Hs.eg.db', ont="MF", computeIC = TRUE)
gsimmat1 = parGOSim(glist, type = 'go', semData = d, measure = 'Wang')
print(gsimmat1)

# by Entrez gene id
genelist = list(c('150', '151', '152', '1814', '1815', '1816'))
d = godata('org.Hs.eg.db', ont="CC", computeIC = TRUE)
gsimmat2 = parGOSim(genelist, type = 'gene', semData = d)
print(gsimmat2)

wind22zhu/BioMedR documentation built on Oct. 21, 2019, 12:51 p.m.