Learn R Programming

LowMACA (version 1.4.2)

lmPlot: Draw a comprehensive LowMACA plot

Description

LowMACA comprehensive plot is a four layers plot that summarize the entire LowMACA output

Usage

lmPlot(object , conservation=NULL, splitLen=NULL)

Arguments

object
a LowMACA class object
conservation
a numeric value in the range of 0-1 that defines the threshold of trident conservation score to include the specified position. The default value is inherited from the slot entropy, whose default is 0.1
splitLen
An integer, defines after how many amino acids the plot should be split By default this parameter is set to NULL, that mean that the plot is not split.

Value

NULL

Details

The method returns a plot, which is divided into four layers. The LowMACA object must have been passed through the methods alignSequences , getMutations , mapMutations and entropy. The four layers of the plot are:
  1. The bar plot visualized by bpAll
  2. The distribution of mutations against the 95% confidence interval superior limit of the null hypothesis (dotted line) with orange bars representing a position with a pvalue <0.05 and="" a="" red="" star="" for="" qvalue<0.05="" <="" li="">
  3. The Trident score distribution
  4. The logo plot representing the consensus sequence

If this plot is used on a LowMACA object with a single protein, the result is formed by three layers only:

  1. The bar plot visualized by bpAll
  2. The Pfam domains structure inside the protein
  3. The distribution of mutations against the 95% confidence interval superior limit of the null hypothesis (dotted line) with orange bars representing a position with a pvalue <0.05 and="" a="" red="" star="" for="" qvalue<0.05="" <="" li="">

See Also

alignSequences getMutations mapMutations entropy bpAll

Examples

Run this code
#Load homeobox example and draw the plot
data(lmObj)
#Calculate statistics for nullProfile
lmObj <- entropy(lmObj)
lmPlot(lmObj)

Run the code above in your browser using DataLab