| plot.MethComp {MethComp} | R Documentation |
Plots the pairwise conversion formulae between methods from a
MethComp object.
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,
... )
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. |
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 standard error.
data( hba1c )
str( hba1c )
hba1c <- transform( subset( hba1c, type=="Ven" ),
meth = dev,
repl = d.ana )
# options( bugs.directory="c:/Program Files/WinBUGS14/" )
# options( bugs.directory="c:/Stat/Bugs/WinBUGS14/")
# hb.res <- MethComp( hba1c, n.iter=50 )
data( hba.MC )
str( hba.MC )
par( ask=TRUE )
plot( hba.MC )
plot( hba.MC, pl.obs=TRUE )