Plot one gene / one peak scatter plot.
plotPeakGene(peak, gene, atac_matrix, rna_matrix, color = NA, legend = TRUE)
The peak coordinate you want to plot.
The gene ID you want to plot.
The normalized ATAC-seq data frame obtained by **quantification**. Same as **getPeak2Gene** parameter.
The normalized RNA-seq dataframe. Same as **getPeak2Gene** parameter.
The color to plot.
Wheather to plot legend. default: TRUE.
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