raw.plot: Raw plot for descriptive statistics

In BivRegBLS: Tolerance Interval and EIV Regression - Method Comparison Studies

Description

Display a plot with the raw data in an (X,Y) plot or (M,D) plot with or witout error bars.

Usage

raw.plot(data.plot = NULL, xname = "X", yname = "Y", graph = "XY.means",
         col.ID = NULL, pch.ID = NULL, ...)

Arguments

data.plot	a desc.stat class object (obtained by using the desc.stat function).
xname	a character string for the name of the X device.
yname	a character string for the name of the Y device.
graph	a character string for the kind of graph. Options available are: XY.means, XY.points, XY.bar.range, XY.bar.SEM, XY.bar.SD, MD.means.
col.ID	a numeric or character vector describing the color of the points per ID.
pch.ID	a numeric vector describing the type of points per ID to plot the data.
...	the common arguments from 'plot' or 'par' that may be used additionaly, such as xlim, ylim, xlab, ylab.

Details

The data.plot argument is mandatory. The labels of the X and Y axes are built by default with xname and yname, and with the type of graph. The arguments col.ID and pch.ID are useful if the argument IDcol was used in the function desc.stat. col.ID and pch.ID are recycled if shorter than the number of IDs, i.e. if 3 colors are specified as col.ID=c(1,2,3) while there are 7 IDs, then the colors will be used as c(1,2,3,1,2,3,1).
With the argument graph, the option XY.means plots the mean measures in a (X,Y) plot, XY.points plots all the X (Y) replicated data centered on the means of Y (X), XY.bar.range plots the error bars from the minimum to the maximum (the range) of the replicates, XY.bar.SEM plots the error bars with the standard errors of the means, XY.bar.SD plots the error bars with the standard deviations, MD.means plots the mean measures in a (M,D) plot.

Value

A plot in a new window.

Note

The limits of the axes and the labels are set automatically by the function. To compare different plots with fixed limits, use xlim and ylim. To write customized labels, use xlab and ylab.

Author(s)

Bernard G FRANCQ

References

Francq BG, Govaerts BB. Measurement methods comparison with errors-in-variables regressions. From horizontal to vertical OLS regression, review and new perspectives. Chemometrics and Intelligent Laboratory Systems 2014; 134:123-139.

See Also

XY.plot, MD.plot

Examples

library(BivRegBLS)
data(Aromatics)
# Calculate the descriptive statistics
res=desc.stat(data=Aromatics,xcol=3,ycol=4,IDcol="Type")
# Plot the mean or single measures (it is the same for unreplicated data)
raw.plot(data.plot=res,xname="HPLC",yname="GC MS",graph="XY.means")
raw.plot(data.plot=res,xname="HPLC",yname="GC MS",graph="XY.points")
# Plot with customized colours and type of points per type of samples
raw.plot(data.plot=res,xname="HPLC",yname="GC MS",graph="XY.points",pch.ID=c(19,5,8),col.ID=c(1,2))
raw.plot(data.plot=res,xname="HPLC",yname="GC MS",graph="MD.means")
raw.plot(data.plot=res,xname="HPLC",yname="GC MS",graph="MD.means",col.ID=c(1,2,4))
data(SBP)
# Calculate the descriptive statistics
res=desc.stat(data=SBP,xcol=c("J1","J2","J3"),ycol=8:10)
# Different plots to show the measurement uncertainties
raw.plot(data.plot=res,xname="J",yname="S",graph="XY.means")
raw.plot(data.plot=res,xname="J",yname="S",graph="XY.points")
raw.plot(data.plot=res,xname="J",yname="S",graph="XY.bar.range")
raw.plot(data.plot=res,xname="J",yname="S",graph="XY.bar.SD")
raw.plot(data.plot=res,xname="J",yname="S",graph="XY.bar.SEM")
raw.plot(data.plot=res,xname="J",yname="S",graph="MD.means")

Example output

Loading required package: ellipse

Attaching package: ‘ellipse’

The following object is masked from ‘package:graphics’:

    pairs

BivRegBLS documentation built on Oct. 11, 2019, 1:05 a.m.