4 Data Visualization
ApplyPolygenicScore provides three plotting functions designed to
operate on the output of apply.polygenic.score:
create.pgs.density.plotcreate.pgs.with.continuous.phenotype.plotcreate.pgs.rank.plot
These visualizations are intended to provide quick quality assessments of PGS in a cohort. They really shine when combined with phenotype data to produce at-a-glance summaries of the relationship between PGS and phenotypes.
4.1 Common plotting arguments
Plotting functions are built on the BoutrosLab.plotting.general
package and produce lattice multipanel plot objects. Each function can
be provided with an output.directory, a
filename.prefix, and a file.extension to write
the plots to a file. A filename will be generated using the current date
and the provided prefix, and the plot will be written in the format
specified by the extension.
Note: in the examples below, all plots are written to a temporary directory from which they are read and displayed inline in the vignette.
Plotting functions automatically detect and plot data from the
standardized PGS output columns produced by
apply.polygenic.score: PGS,
PGS.with.normalized.missing, and
PGS.with.replaced.missing. Plots are automatically titled
with the name of the PGS column being plotted, and the titles can be
“tidied up” to remove period characters and replace them with spaces
using the tidy.titles argument.
Multipanel plot dimensions can be controlled with width
and height arguments, which are interpreted as inches. If
your plots are looking squished, try increasing these values. Axis and
title label font sizes can be controlled with xaxis.cex,
yaxis.cex, and titles.cex arguments. Border
padding can be controlled with border.padding and is
applied to all four sides of the plot at once.
4.2 PGS Density
create.pgs.density.plot plots the PGS distribution in
the cohort.
4.2.1 Basic plot
create.pgs.density.plot(
pgs.data = pgs.results.with.phenotype.analysis$pgs.output,
output.dir = temp.dir,
filename.prefix = 'vignette-example-basic',
file.extension = 'png'
)## NULL
4.2.2 Add phenotypes
If provided with phenotype column names, this function will plot individual PGS density curves for each category in a categorical phenotype. - Non-categorical phenotypes are ignored. - Colors are automatically assigned up to 12 categories. - If more than 12 categories are present, all categories are plotted in black. - When multiple PGS columns are present and/or multiple phenotype columns are specified, each combination is plotted in a separate panel. - Panels are arranged in a grid with unique PGS inputs on the x-axis and unique phenotype inputs on the y-axis.
create.pgs.density.plot(
pgs.data = pgs.results.with.phenotype.analysis$pgs.output,
output.dir = temp.dir,
filename.prefix = 'vignette-example-phenotype',
file.extension = 'png',
tidy.titles = TRUE,
phenotype.columns = c('binary.factor.phenotype', 'categorical.phenotype', 'continuous.phenotype')
)## NULL
4.3 PGS Correlation
create.pgs.with.continuous.phenotype.plot plots a
scatterplot between the PGS and a continuous phenotype, and reports a
correlation. Non-continuous phenotypes are ignored. This function can
only be used if phenotype data is provided.
4.3.1 Basic plot
create.pgs.with.continuous.phenotype.plot(
pgs.data = pgs.results.with.phenotype.analysis$pgs.output,
output.dir = temp.dir,
filename.prefix = 'vignette-example-correlation',
file.extension = 'png',
tidy.titles = TRUE,
phenotype.columns = c('continuous.phenotype', 'categorical.phenotype')
)## NULL
Unsurprisingly, our randomly generated phenotype data is not correlated with our fake PGS.
- When multiple PGS columns are present and/or multiple phenotype columns are specified, each combination is plotted in a separate panel.
- Panels are arranged in a grid with unique PGS inputs on the x-axis and unique phenotype inputs on the y-axis.
4.4 PGS Percentile Rank
create.pgs.rank.plot plots several panels of
per-individual data, each in the order of increasing PGS percentile
rank.
At the top, a barplot of missing PGS variant genotype counts or percents for each individual (only if at least one individual is missing a variant). Next, a barplot of the percentile rank of each individual. Next, a covariate bar of the PGS quartile and decile of each individual. Next, a covariate bar of categories for each provided categorical phenotype in each individual. Finally, a heatmap of continuous phenotype values for each provided continuous phenotype in each individual.
4.4.1 Basic plot
# Introducing some missing genotypes to demonstrate the missing genotype barplot
pgs.results.with.phenotype.analysis$pgs.output$n.missing.genotypes <- rep(c(0, 1, 0, 2, 1), 2)
create.pgs.rank.plot(
pgs.data = pgs.results.with.phenotype.analysis$pgs.output,
output.dir = temp.dir,
filename.prefix = 'vignette-example-basic',
file.extension = 'png'
)## NULL
4.4.2 Add phenotypes
create.pgs.rank.plot(
pgs.data = pgs.results.with.phenotype.analysis$pgs.output,
output.dir = temp.dir,
filename.prefix = 'vignette-example-phenotype',
file.extension = 'png',
phenotype.columns = c('binary.factor.phenotype', 'binary.phenotype', 'categorical.phenotype', 'continuous.phenotype')
)## [1] "springgreen3"
## [1] "darkturquoise"
## [1] "deepskyblue"## NULL
- Legend titles are inherited from phenotype column names.
- For cohorts exceeding 50 individuals, sample labels are automatically removed.
- If missing phenotypes exist, they are plotted in gray and accompanied by a missingness legend.
4.4.3 Optional arguments
This function comes with a handful of additional arguments to control the appearance of the plot.
- Use
missing.genotype.styleto toggle between displaying counts and percentages. - Use
categorical.paletteto specify a custom color palette for categorical phenotype bars. - Use
binary.palletteto specify a custom color palette for binary and continuous phenotype bars.