plot_density: Create density plots for mixture models

In tidyLPA: Easily Carry Out Latent Profile Analysis (LPA) Using Open-Source or Commercial Software

Description

Creates a faceted plot of density plots for an object of class 'tidyLPA'. For each variable, a Total density plot will be shown, along with separate density plots for each latent class, where cases are weighted by the posterior probability of being assigned to that class.

Usage

plot_density(
  x,
  variables = NULL,
  bw = FALSE,
  conditional = FALSE,
  alpha = 0.2,
  facet_labels = NULL
)

Arguments

x	Object to plot.
variables	Which variables to plot. If NULL, plots all variables that are present in all models.
bw	Logical. Whether to make a black and white plot (for print) or a color plot. Defaults to FALSE, because these density plots are hard to read in black and white.
conditional	Logical. Whether to show a conditional density plot (surface area is divided amongst the latent classes), or a classic density plot (surface area of the total density plot is equal to one, and is subdivided amongst the classes).
alpha	Numeric (0-1). Only used when bw and conditional are FALSE. Sets the transparency of geom_density, so that classes with a small number of cases remain visible.
facet_labels	Named character vector, the names of which should correspond to the facet labels one wishes to rename, and the values of which provide new names for these facets. For example, to rename variables, in the example with the 'iris' data below, one could specify: facet_labels = c("Pet_leng" = "Petal length").

Value

An object of class 'ggplot'.

Author(s)

Caspar J. van Lissa

Examples

## Not run: 
results <- iris %>%
  subset(select = c("Sepal.Length", "Sepal.Width",
    "Petal.Length", "Petal.Width")) %>%
  estimate_profiles(1:3)

## End(Not run)
## Not run: 
plot_density(results, variables = "Petal.Length")

## End(Not run)
## Not run: 
plot_density(results, bw = TRUE)

## End(Not run)
## Not run: 
plot_density(results, bw = FALSE, conditional = TRUE)

## End(Not run)
## Not run: 
plot_density(results[[2]], variables = "Petal.Length")

## End(Not run)

Example output

You can use the function citation('tidyLPA') to create a citation for the use of {tidyLPA}.
Mplus is not installed. Use only package = 'mclust' when calling estimate_profiles().
Warning message:
In system2("type", args = "mplus", stdout = FALSE, stderr = FALSE) :
  error in running command

tidyLPA documentation built on Nov. 17, 2021, 5:07 p.m.