DDPNA (version 0.4.1)

ID_match: homolog protein Uniprot ID transformation

Description

homolog protein Uniprot ID match

Usage

ID_match(data, db1.path = NULL, db2.path = NULL,out.folder = NULL,
         blast.path = NULL,evalue = 0.1, verbose = 1)

Value

a data.frame included 4 columns: ori.ID, ENTRY.NAME, new.ID, match.type.

Arguments

data

dataset of protein information.Column Names should contain "ori.ID" and "ENTRY.NAME". "ori.ID" is Uniprot ID

db1.path

fasta file, database of transfered species

db2.path

fasta file, database of original species

out.folder

blast result output folder, the folder path should be the same with db1.path

blast.path

blast+ software install path

evalue

blast threshold, the lower means more rigorous

verbose

integer level of verbosity. Zero means silent, 1 means have Diagnostic Messages.

Author

Kefu Liu

Details

homolog protein Uniprot ID match is based on the ENTRY.NAME, gene name and sequence homophyly in two different species or different version of database.

Examples

# suggested to install blast+ software
# \donttest{
# it will take a long time without blast+ software
data(Sample_ID_data)
if(requireNamespace("Biostrings", quietly = TRUE)){
  out.folder = tempdir();
  write.table(Sample_ID_data$db1,file.path(out.folder,"db1.fasta"),
              quote = FALSE,row.names = FALSE, col.names = FALSE);
  write.table(Sample_ID_data$db2,file.path(out.folder,"db2.fasta"),
              quote = FALSE,row.names = FALSE, col.names = FALSE);
  data <- ID_match(Sample_ID_data$ID_match_data,
                   db1.path = file.path(out.folder,"db1.fasta"),
                   db2.path = file.path(out.folder,"db2.fasta"),
                   out.folder = out.folder,
                   blast.path = NULL,
                   evalue = 0.1, verbose = 1)
 file.remove( file.path(out.folder,"db1.fasta"),
              file.path(out.folder,"db2.fasta"))
}
# }