Plot PCA

Plot multi-dimensional scaling plot using algorithm of BiocSingular::runPCA(). It is recommended this be done with the log-methylation-ratio matrix generated by bsseq_to_log_methy_ratio().

Usage

plot_pca(
  x,
  plot_dims = c(1, 2),
  labels = colnames(x),
  groups = NULL,
  legend_name = "group"
)

Arguments

x

the log-methylation-ratio matrix.

plot_dims

the numeric vector of the two dimensions to be plotted.

labels

the character vector of labels for data points. By default uses column names of x, set to NULL to plot points.

groups

the character vector of groups the data points will be coloured by.

legend_name

the name for the legend.

Value

ggplot object of the MDS plot.

Examples

nmr <- load_example_nanomethresult()
bss <- methy_to_bsseq(nmr)
#> [2024-08-22 04:06:06] creating intermediate files...
#> [2024-08-22 04:06:06] parsing chr11...
#> [2024-08-22 04:06:06] parsing chr12...
#> [2024-08-22 04:06:06] parsing chr18...
#> [2024-08-22 04:06:06] parsing chr5...
#> [2024-08-22 04:06:06] parsing chr7...
#> [2024-08-22 04:06:06] parsing chrX...
#> [2024-08-22 04:06:06] samples found: B6Cast_Prom_3_cast B6Cast_Prom_3_bl6 B6Cast_Prom_2_cast B6Cast_Prom_2_bl6 B6Cast_Prom_1_cast B6Cast_Prom_1_bl6 
#> [2024-08-22 04:06:06] creating bsseq object...
#> [2024-08-22 04:06:06] reading in parsed data...
#> [2024-08-22 04:06:06] constructing matrices...
#> [2024-08-22 04:06:07] done
lmr <- bsseq_to_log_methy_ratio(bss)
plot_pca(lmr)