TraMineR (version 2.2-0.1)

plot.stslist.meant: Plot method for objects produced by the seqmeant function

Description

This is the plot method for objects of class stslist.meant produced by the seqmeant function.

Usage

# S3 method for stslist.meant
plot(x, cpal = NULL, ylab = NULL, yaxis = TRUE,
  xaxis = TRUE, cex.axis = 1, ylim = NULL, cex.plot, ...)

Arguments

x

an object of class stslist.meant as produced by the seqmeant function.

cpal

alternative color palette to use for the states. If user specified, a vector of colors with number of elements equal to the number of states in the alphabet. By default, the 'cpal' attribute of the 'seqdata' sequence object is used (see seqdef).

ylab

an optional label for the y axis. If set to NA, no label is drawn.

yaxis

controls whether the y axis is plotted. Default is TRUE.

xaxis

if TRUE (default) the xaxis is plotted.

cex.axis

Axis annotation magnification. See par.

ylim

an optional vector setting the limits for the y axis. If NULL (default), limits are set to (0, max. sequence length).

cex.plot

Deprecated. Use cex.axis instead.

further graphical parameters. For more details about the graphical parameter arguments, see barplot and par.

Details

This is the plot method for the output produced by the seqmeant function, i.e., objects of class stslist.meant. It produces a plot showing the mean times spent in each state of the alphabet.

When the "se" attribute of x is TRUE, i.e., when x contains also the standard errors of the mean times, error bars are automatically displayed on the plot. See the serr argument of seqmeant.

This method is called by the generic seqplot function (if type="mt") that produces more sophisticated plots, allowing grouping and automatic display of the states legend. The seqmtplot function is a shortcut for calling seqplot with type="mt".

Examples

Run this code
# NOT RUN {
## Loading the mvad data set and creating a sequence object
data(mvad)
mvad.labels <- c("employment", "further education", "higher education",
                    "joblessness", "school", "training")
mvad.scodes <- c("EM","FE","HE","JL","SC","TR")
mvad.seq <- seqdef(mvad, 15:86, states=mvad.scodes, labels=mvad.labels)

## Computing the mean times
mvad.meant <- seqmeant(mvad.seq)

## Plotting
plot(mvad.meant, main="Mean durations in each state of the alphabet")

## Changing the y axis limits
plot(mvad.meant, main="Mean durations in each state of the alphabet",
 ylim=c(0,40))

## Displaying error bars
mvad.meant.e <- seqmeant(mvad.seq, serr=TRUE)
plot(mvad.meant.e, main="Mean durations in each state of the alphabet",
 ylim=c(0,40))

# }

Run the code above in your browser using DataCamp Workspace