Description Format Details Source References Examples
Effectiveness of 3 antibiotics against 16 bacterial species.
A data frame with 16 observations on the following 5 variables.
bacteria
bacterial species, 16 levels
penicillin
MIC for penicillin
streptomycin
MIC for streptomycin
neomycin
MIC for neomycin
gramstain
Gram staining (positive or negative)
The values reported are the minimum inhibitory concentration (MIC) in micrograms/milliliter, which represents the concentration of antibiotic required to prevent growth in vitro.
Will Burtin (1951). Scope. Fall, 1951.
Wainer, H. (2009). A Centenary Celebration for Will Burtin: A Pioneer of Scientific Visualization. Chance, 22(1), 51-55. https://chance.amstat.org/2009/02/visrev221/
Wainer, H. (2009). Visual Revelations: Pictures at an Exhibition. Chance, 22(2), 46–54. https://chance.amstat.org/2009/04/visrev222/
Wainer, H. (2014). Medical Illuminations: Using Evidence, Visualization and Statistical Thinking to Improve Healthcare.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 | data(antibiotic)
lucid(antibiotic)
## Not run:
# Plot the data similar to Fig 2.14 of Wainer's book, "Medical Illuminations"
require(lattice)
require(reshape2)
# Use log10 transform
dat <- transform(antibiotic,
penicillin=log10(penicillin),
streptomycin=log10(streptomycin),
neomycin=log10(neomycin))
dat <- transform(dat, sgn = ifelse(dat$gramstain=="neg", "-", "+"))
dat <- transform(dat,
bacteria = paste(bacteria, sgn))
dat <- transform(dat, bacteria=reorder(bacteria, -penicillin))
dat <- melt(dat)
op <- tpg <- trellis.par.get()
tpg$superpose.symbol$pch <- toupper(substring(levels(dat$variable),1,1))
tpg$superpose.symbol$col <- c("darkgreen","purple","orange")
trellis.par.set(tpg)
dotplot(bacteria ~ value, data=dat, group=variable,
cex=2,
scales=list(x=list(at= -3:3,
labels=c('.001', '.01', '.1', '1', '10', '100', '1000'))),
main="Bacterial response to Neomycin, Streptomycin, and Penicillin",
xlab="Minimum Inhibitory Concentration (mg/L)")
trellis.par.set(op)
## End(Not run)
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.