Description Usage Arguments Value Author(s) Examples

Plots latent and observed trajectories for a list of growth mixture models, specified using Mplus' '|' notation (e.g., i s | ).

1 2 3 4 |

`modelList` |
A list object of Mplus models, or a single Mplus model. This function additionally requires the models to be growth mixture models. |

`bw` |
Logical. Should the plot be black and white (for print), or color? |

`rawdata` |
Logical. Should raw data (observed trajectories) be plotted in the background? Setting this to TRUE might result in long plotting times. Requires including the Mplus syntax 'SAVEDATA: FILE IS "filename"; SAVE = cprobabilities' in the Mplus input. |

`estimated` |
Logical. Should the Mplus estimates growth trajectories be displayed? Defaults to TRUE. |

`poly` |
Logical. Should polynomial smooth lines be displayed for each group? Defaults to FALSE. If set to TRUE, the order of the polynomial is determined by taking the number of observed timepoints - 1 (e.g., 3 time points results in a quadratic polynomial). The data are weighted according to their posterior class probabilities. Note that rawdata must be TRUE in order to obtain polynomial smooth lines, because these are calculated on the raw data. |

`alpha_range` |
Numeric vector. The minimum and maximum values of alpha (transparancy) for the raw data. Minimum should be 0; lower maximum values of alpha can help reduce overplotting. |

`growth_variables` |
Character vector. Indicate the names of the latent variables for the growth trajectory to plot. If NULL (default), all latent growth variables are used. Use this option to plot one trajectory when a model contains multiple latent growth trajectories. |

`time_scale` |
Numeric vector. In case some of the loadings of the growth model are freely estimated, provide the correct time scale here (e.g., c(0, 1, 2)). |

`jitter_lines` |
Numeric. Indicate the amount (expressed in fractions of a standard deviation of all observed data) by which the observed trajectories should be vertically jittered. Like alpha_range, this parameter helps control overplotting. |

An object of class 'ggplot'.

Caspar J. van Lissa

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 | ```
## Not run:
mydat <- read.table("http://statmodel.com/usersguide/chap8/ex8.2.dat", header = FALSE)[,-6]
createMixtures(classes = 1:3, filename_stem = "ex8.2",
model_overall = "i s | [email protected] [email protected] [email protected] [email protected]; i s on V5;",
rdata = mydat,
OUTPUT = "tech11 tech14;", usevariables = c("V1", "V2", "V3",
"V4", "V5"))
runModels(replaceOutfile = "modifiedDate")
mplus_output <- readModels(filefilter = "ex8.2")
plotGrowthMixtures(mplus_output, estimated = TRUE, rawdata = TRUE, time_scale = c(0, 1, 2, 3))
plotGrowthMixtures(mplus_output, estimated = FALSE, rawdata = TRUE,
poly = TRUE)
createMixtures(classes = 1:3, filename_stem = "ex8.2_free",
model_overall = "i s | [email protected] V2* V3* [email protected]; i s on V5;",
rdata = mydat,
OUTPUT = "tech11 tech14;", usevariables = c("V1", "V2", "V3",
"V4", "V5"))
runModels(replaceOutfile = "modifiedDate")
mplus_output_free <- readModels(filefilter = "ex8.2_free")
plotMixtureDensities(mplus_output_free, variables = c("V1", "V2", "V3",
"V4"))
plotGrowthMixtures(mplus_output_free, estimated = TRUE, rawdata = TRUE,
bw = TRUE, time_scale = c(0, 1, 2, 3),
alpha_range = c(0, .05))
plotGrowthMixtures(mplus_output, estimated = FALSE, rawdata = TRUE,
poly = TRUE, bw = TRUE, time_scale = c(0, 1, 2, 3))
## End(Not run)
``` |

Embedding an R snippet on your website

Add the following code to your website.

For more information on customizing the embed code, read Embedding Snippets.