PLNmodels (version 1.2.2)

plot.PLNmixturefit: Mixture visualization of a PLNmixturefit object

Description

Represent the result of the clustering either by coloring the individual in a two-dimension PCA factor map, or by representing the expected matrix of count reorder according to the clustering.

Usage

# S3 method for PLNmixturefit
plot(x, type = c("pca", "matrix"), main = NULL, plot = TRUE, ...)

Value

a ggplot2::ggplot graphic

Arguments

x

an R6 object with class PLNmixturefit

type

character for the type of plot, either "pca", for or "matrix". Default is "pca".

main

character. A title for the plot. If NULL (the default), an hopefully appropriate title will be used.

plot

logical. Should the plot be displayed or sent back as ggplot2::ggplot object

...

Not used (S3 compatibility).

Examples

data(trichoptera)
trichoptera <- prepare_data(trichoptera$Abundance, trichoptera$Covariate)
myPLN <- PLNmixture(Abundance ~ 1 + offset(log(Offset)),
           data = trichoptera, control = PLNmixture_param(smoothing = "none"))  %>% getBestModel()
if (FALSE) {
plot(myPLN, "pca")
plot(myPLN, "matrix")
}