plotOccupancyProfile {ChIPanalyser} | R Documentation

## Plot Occupancy Profiles

### 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

      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

    #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"))

[Package ChIPanalyser version 1.9.1 Index]