plot.MethComp {MethComp}R Documentation

Plot estimated conversion lines and formulae.

Description

Plots the pairwise conversion formulae between methods from a MethComp object.

Usage

plot.MethComp( x,
             axlim = range( attr(x,"data")$y ),
             which,
          lwd.line = c(3,1), col.line = rep("black",2), lty.line=rep(1,2),
               eqn = TRUE, digits = 2,
              grid = FALSE, col.grid=gray(0.8),
            pl.obs = FALSE,
           col.pts = "black", pch.pts = 16, cex.pts = 0.8,
               ... )

Arguments

x A MethComp object
axlim The limits for the axes in the panels
which Numeric vector or vector of method names. Which of the methods should be included in the plot?
lwd.line Numerical vector of length 2. The width of the conversion line and the prediction limits. If the second values is 0, no prediction limits are drawn.
col.line Numerical vector of length 2. The color of the conversion line and the prediction limits.
lty.line Numerical vector of length 2. The line types of the conversion line and the prediction limits.
eqn Should the conversion equations be printed on the plot?. Defaults to TRUE.
digits How many digits after the decimal point shoudl be used when printing the conversion equations.
grid Should a grid be drawn? If a numerical vector is given, the grid is drawn at those values.
col.grid What color should the grid have?
pl.obs Logical or character. Should the points be plotted. If TRUE or "repl" paired values of single replicates are plotted. If "perm", replicates are randomly permuted within (item, method) befor plotting. If "mean", means across replicates within item, method are formed and plotted.
col.pts What color should the observation have.
pch.pts What plotting symbol should be used.
cex.pts What scaling should be used for the plot symbols.
... Parameters to pass on. Currently not used.

Value

Nothing. The lower part of a (M-1) by (M-1) matrix of plots is drawn, showing the pairwise conversion lines. In the corners of each is given the two conversion equations together with the prediction standar error.

See Also

MethComp, print.MethComp

Examples

data( hba1c )
str( hba1c )
hba1c <- transform( hba1c, meth = interaction(dev,type),
                           repl = d.ana )
# options( bugs.directory="c:/Program Files/WinBUGS14/" )
options( bugs.directory="c:/Stat/Bugs/WinBUGS14/")
hb.res <- MethComp( hba1c, n.iter=100 )
str( hb.res )
par( ask=TRUE )
plot( hb.res )
plot( hb.res, pl.obs=TRUE, which=1:4 )

[Package MethComp version 0.1.15 Index]