Bios2cor (version 1.2)

corr_contact: Calculates the number of contacts for each element of the top pairs

Description

Given a correlation object, calculates the number of pairs (contacts) each element in the top X pairs are involved in

Usage

corr_contact(corr, top= 25, contact_filepath= NULL, score_top_filepath= NULL)

Arguments

corr

An object created by a correlation function (omes, mip, elsc, mcbasc, rotamer_circular, rotamer_omes, rotamer_mip).

top

A integer indicating the number of top pairs used for this analysis

contact_filepath

The full path name of the output contact file

score_top_filepath

The full path name of the output top score file

Value

A vector containing top elements names

Details

This function sorts element pairs by correlation score and analyses the top X pairs to determine the number of pairs (contacts) each element of the top X pairs is involved in.

Examples

# NOT RUN {
  msf <- system.file("msa/toy_align.msf", package = "Bios2cor")
  align <- import.msf(msf)

  #Creating OMES object
  omes <- omes(align,fileHelix= NULL , diag= 0, fileCSV= NULL, gap_val= 0.8, z_score= TRUE)

  omes <-omes$normalized
   
  contact <- corr_contact(omes, top= 25, "contact_file.txt","top_scores.txt")
# }