Plotting the output of pairwiseCI
Easy method for plotting estimates and confidence bounds calculated using pairwiseCI
.
## S3 method for class 'pairwiseCI' plot(x, CIvert=NULL, CIlty = 1, CIlwd=1, CIcex=1, H0line=NULL, H0lty=1, H0lwd=1, main=NULL, ylab="", xlab="", ... )
x
: an object of class "pairwiseCI", the output of function \link{pairwiseCI}CIvert
: logical, whether confidence intervals shall be plotted vertical if CIvert=TRUE and horizontal if CIvert=FALSECIlty
: integer, giving the line type of the CI, as documented for cex in ?parCIlwd
: integer, giving the line width of the CI, as documented for lwd ?parCIcex
: numerical value giving the size of CIsymbols relative to the default value, see cex in ?parH0line
: Value to be plotted as vertical or horizontal line, depending on the value of CIvertH0lty
: integer, giving the line type of the CI, as documented for lty in ?parH0lwd
: integer, giving the line width of the CI, as documented for lwd in ?parmain
: as main in plotylab
: label of y-axis as ylab in plot , default is no labelxlab
: label of x-axis as ylab in plot , default is no label...
: Further arguments to be passed to axis . Note, that arguments las , at , labels are defined internally and can not be set via ...
Frank Schaarschmidt
data(Oats) output <- pairwiseCI(yield ~ Block, data=Oats, by="nitro",method="Param.diff", control="I") # default plot for difference methods: plot(output) # some small changes: plot(output, CIvert=TRUE, H0line=c(-2,0,2), H0lty=c(2,1,2)) output <- pairwiseCI(yield ~ Block, data=Oats, by="nitro", method="Param.ratio", control="I") # default plot for ratio methods: plot(output) # some small changes: plot(output, CIvert=FALSE, H0line=c(0.7, 1, 1/0.7), H0lty=c(3,2,3))
Useful links