lmPlot: Draw a comprehensive LowMACA plot

In ste-depo/LowMACA: LowMACA - Low frequency Mutation Analysis via Consensus Alignment

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.

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

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

Value

NULL

Author(s)

Stefano de Pretis , Giorgio Melloni

See Also

alignSequences getMutations mapMutations entropy bpAll

Examples

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

ste-depo/LowMACA documentation built on Oct. 15, 2022, 11:53 p.m.