plotOccupancyProfile: Plot Occupancy Profiles

Description Usage Arguments Details Value Author(s) References Examples

View source: R/plotOccupancyDev.R

Description

plotOccupancyProfile plots the predicted profiles. If provided, this functions will also plot ChIP-seq profiles, PWMScores (or Occupancy), chromatin States, Goodness of Fit estimates and gene information.

Usage

1
2
3
  plotOccupancyProfile(predictedProfile, ChIPScore = NULL,chromatinState = NULL
      ,occupancy = NULL,goodnessOfFit = NULL,PWM=FALSE,
      geneRef = NULL,axis=TRUE,...)

Arguments

predictedProfile

predictedProfile is a either GRanges containing the predicted profiles for one loci, all loci selected for one paramter, or all loci selected for all parameter combinations selected. (see searchSites)

ChIPScore

ChIPScore is a ChIPscore object containing ChIPscore (or a list of numeric values representing ChIP scores (Experimental ChIP))

chromatinState

chromatinState is a GRanges containing accesible DNA sites or chromatin States.

occupancy

occupancy is a GRanges or a genomicProfiles object contaning PWM scores and Occupancy ( see computeOccupancy)

goodnessOfFit

goodnessOfFit results of the profileAccuracyEstimate function.

PWM

PWM is a logical value that in the case occupancy is provided which of occupancy scores of PWM scores hsould be plotted. Default set at FALSE

geneRef

geneRef is a GRanges containing gene information on exons,introns, UTR's, enhancers or any other genetic element to be plotted.

axis

axis is a logical value indicating if axes should be plotted or not.

...

Any other graphical Parameter of the following : col, density, border, lty, lwd, cex, cex.axis, xlab, ylab, xlim, ylim, las and axislables. Each Parameter will be parsed in the same order as the arguments to this function. If the name of the argument is specified, the argument value will be parsed to the correct internal plotting. See examples below

Details

Once the predicted ChIP-seq like profiles have been computed, it is possible to plot these profiles.

Value

Returns a profile plot with "Occupancy" on the y axis and DNA position on the the X- axis. If the orange line is the predicted profile, the grey shaded area represents real ChIP-seq data. The yellow boxes represent regions on NON-accessible DNA. On the lower part of the plot, gene information is plotted with respect to the strand they are localised on. Finally, the blue vertical lines represent sites of either high Occupancy or high PWM score depending on which had been selected. The minimal plot will only contain the predicted profile. The more data is provided the more will be plotted. IMPORTANT: the colours described above represent the default colour values. If colours are provided by the user, they will be updated in plotting order or if the argument has been speicified assigned to said argument.

Author(s)

Patrick C.N. Martin <pm16057@essex.ac.uk>

References

Zabet NR, Adryan B (2015) Estimating binding properties of transcription factors from genome-wide binding profiles. Nucleic Acids Res., 43, 84–94.

Examples

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
#Data extraction
data(ChIPanalyserData)
# path to Position Frequency Matrix
PFM <- file.path(system.file("extdata",package="ChIPanalyser"),"BCDSlx.pfm")
#As an example of genome, this example will run on the Drosophila genome

if(!require("BSgenome.Dmelanogaster.UCSC.dm3", character.only = TRUE)){
    if (!requireNamespace("BiocManager", quietly=TRUE))
        install.packages("BiocManager")
    BiocManager::install("BSgenome.Dmelanogaster.UCSC.dm3")
    }
library(BSgenome.Dmelanogaster.UCSC.dm3)
DNASequenceSet <- getSeq(BSgenome.Dmelanogaster.UCSC.dm3)
#Building data objects
GPP <- genomicProfiles(PFM=PFM,BPFrequency=DNASequenceSet)


# Computing Genome Wide
GenomeWide <- computeGenomeWideScores(DNASequenceSet = DNASequenceSet,
    genomicProfiles = GPP)

#Compute PWM Scores
PWMScores <- computePWMScore(DNASequenceSet = DNASequenceSet,
    genomicProfiles = GenomeWide,
    loci = eveLocus, chromatinState = Access)
#Compute Occupnacy
Occupancy <- computeOccupancy(genomicProfiles = PWMScores)

#Compute ChIP profiles
chipProfile <- computeChipProfile(loci = eveLocus,
    genomicProfiles = Occupancy)


#Plotting Profile
plotOccupancyProfile(predictedProfile=chipProfile,
    ChIPScore = eveLocusChip,
    chromatinState = Access,
    occupancy = Occupancy,
    geneRef =geneRef)

## Changing graphical Parameters

#### In this examples ####
### predictedProfile will be red
### chipProfile will be blue
### and DNAAccessibility will be green
plotOccupancyProfile(predictedProfile=chipProfile,
    ChIPScore = eveLocusChip,
    chromatinState = Access,
    occupancy = Occupancy,
    geneRef =geneRef,col=c("red","blue","green"))

### If name is specified
plotOccupancyProfile(predictedProfile=chipProfile,
    ChIPScore = eveLocusChip,
    chromatinState = Access,
    occupancy = Occupancy,
    geneRef =geneRef,col=c("DNAAccessibility"="red","blue","green"))

patrickCNMartin/ChIPanalyserDev documentation built on Nov. 14, 2019, 4:39 p.m.