MDplot (version 1.0.1)

timeseries: General timeseries plotting function

Description

Plot one (or more) timeseries plots.

Usage

timeseries( tsData,
            printLegend = TRUE,
            snapshotsPerTimeInt = 1000,
            timeUnit = "ns",
            valueName = NA,
            valueUnit = NA,
            colours = NA,
            names = NA,
            legendPosition = "bottomright",
            barePlot = FALSE,
            ... )

Arguments

tsData

List of (alternating) indices and response values, as also produced by load_timeseries() for example.

printLegend

Boolean, which triggers plotting of the legend.

snapshotsPerTimeInt

Number, specifying how many snapshots are within one timeUnit.

timeUnit

Specifies, which unit the x-axis is given in.

valueName

Name of response variable.

valueUnit

Specifies, which unit the y-axis is given in.

colours

Vector of colours used for plotting.

names

Vector of the names of the trajectories.

legendPosition

Indicate position of legend: either "bottomright", "bottomleft", "topleft" or "topright".

barePlot

Boolean, indicating whether the plot is to be made without any additional information.

...

Additional arguments (ellipsis).

Value

Returns a list of list, the latter each holding for every data input series:

• maxValue The maximum value over the whole set.

• meanValue The mean value over the whole set.

• sd The standard deviation over the whole set.

Examples

# GROMOS
timeseries( load_timeseries( c( system.file( "extdata/timeseries_example_1.txt.gz",
                                             package = "MDplot" ),
                                system.file( "extdata/timeseries_example_2.txt.gz",
                                             package = "MDplot" ) ) ),
            snapshotsPerTimeInt = 100 )

# GROMACS
timeseries( load_timeseries( c( system.file( "extdata/timeseries_example_GROMACS.txt.gz",
                                             package = "MDplot" ) ),
                             mdEngine = "GROMACS" ),
            ylim = c( 0.545, 0.7 ), valueName = "Area per lipid", valueUnit = "nm^2" )