Plot one gene / one peak scatter plot.

plotPeakGene(peak, gene, atac_matrix, rna_matrix, color = NA, legend = TRUE)

Arguments

peak

The peak coordinate you want to plot.

gene

The gene ID you want to plot.

atac_matrix

The normalized ATAC-seq data frame obtained by **quantification**. Same as **getPeak2Gene** parameter.

rna_matrix

The normalized RNA-seq dataframe. Same as **getPeak2Gene** parameter.

color

The color to plot.

legend

Wheather to plot legend. default: TRUE.

Examples

  plotPeakGene("1:39081-39341","gene","ATAC_norm_quant.tsv","RNA_norm_quant.tsv",legend=F)
#> Warning: cannot open file 'ATAC_norm_quant.tsv': No such file or directory
#> Error in file(file, "rt"): cannot open the connection