bland.altman.plot: Produce Bland-Altman Plot Bland-AltmanPlots for assessing...

In BlandAltmanLeh: Plots (Slightly Extended) Bland-Altman Plots

Description

Produce Bland-Altman Plot

Bland-AltmanPlots for assessing agreement between two measuring methods or repeatability (test-retest agreement) of measurements. Using either base graphics or ggplot2.

Usage

bland.altman.plot(group1, group2, two = 1.96, mode = 1,
  graph.sys = "base", conf.int = 0, silent = TRUE, sunflower = FALSE,
  geom_count = FALSE, ...)

Arguments

group1	Measurements with first method or first measurement
group2	Measurements with second method or second measurement
two	Lines are drawn "two" standard deviations from mean differences. This defaults to 1.96 for proper 95 percent confidence interval estimation but can be set to 2.0 for better agreement with e. g. the Bland Altman publication.
mode	if 1 then difference group1 minus group2 is used, if 2 then group2 minus group1 is used. Defaults to 1.
graph.sys	Graphing system within R. This defaults to "base" but can be one out of c("base", "ggplot2"), providing ggplot2 is installed.
conf.int	Defaults to 0 which draws the usual Bland Altman plot which contains no confidence intervalls. Change to .95 for 95 percent confidence intervalls to be drawn.
silent	logical. If graph.sys=="base" and silent==TRUE then no return value. If graph.sys=="base" and silent==FALSE then returns statistics.
sunflower	logical. If TRUE, the plot will be based on a sunflower plot and ties will be marked accordingly. Try with data with ties. Works only with graph.sys=="base".
geom_count	logical. If TRUE, the dots will get larger the more frequent given pair is. Use in presence of ties. Works only with graph.sys=="ggplot2" version >=2.0.0.
...	passed on to graphics functions if graph.sys=="base"

Value

Depends on graphic system chosen. In case of "base" depending on whether silent==TRUE. If silent==TRUE then no returns. If silent==FALSE than returns list of statistics as returned by bland.altman.stats(). In case the graphics system is "ggplot2" than the graphic object is returned so that it can be printed or altered.

Author(s)

Bernhard Lehnert <[e-mail]>

See Also

bland.altman.stats

Examples

bland.altman.plot(rnorm(20), rnorm(20), xlab="mean measurement", 
                  ylab="differences", main="Example plot")
                  
bland.altman.plot(rnorm(20), 2+.8*rnorm(20), xlab="mean measurement", 
                  ylab="differences", conf.int=.95)
                  
bland.altman.plot(rnorm(200), 2+.8*rnorm(200), xlab="mean measurement", 
                  ylab="differences", conf.int=.95)
                  
# this is what fig.2 in Bland&Altman1986 would have looked like
PEFR1 <- bland.altman.PEFR[,1]
PEFR2 <- bland.altman.PEFR[,3]
bland.altman.plot(PEFR1, PEFR2, silent=TRUE, xlim=c(0,800),
                  xlab="Average PEFR by two meters",
                  ylab="Difference in PEFR (large-mini)")

# and this is the same but with additional 95 percent CIs
data(bland.altman.PEFR)
bland.altman.plot(PEFR1, PEFR2, silent=TRUE, conf.int=.95, xlim=c(0,800))
                  
# an example with many ties and the 'sunflower'-option
a <- rep(c(1,1,2,2,2,2,2,2,2,2,3,3,3,3,3,3,3,3,4,4,4,5,6,6),2)
b <- rep(c(1,1,1,2,2,2,3,1,4,2,5,3,3,3,3,3),3)
bland.altman.plot(a,b,sunflower=TRUE, xlab="Mean",ylab="Difference",
                  main="discrete values lead to ties")
                  
library(ggplot2)
a <- bland.altman.plot(rnorm(20), rnorm(20), graph.sys="ggplot2", conf.int=.9)
print(a + xlab("you can change this later") + ggtitle("Title goes here"))                  

Example output

NULL
NULL
NULL
NULL
Warning message:
In data(bland.altman.PEFR) : data set 'bland.altman.PEFR' not found
NULL
NULL

BlandAltmanLeh documentation built on May 1, 2019, 8:04 p.m.