In PhilipPallmann/toxbox: Boxplots for Toxicological Data

library("knitr")
library("toxbox")

We use the rat pup data from package WWGbook and reorder the treatment factor:

library("WWGbook")
data(ratpup)
ratpup$treatment <- factor(ratpup$treatment, levels=c("Control", "Low", "High"))

We plot the data assuming the Dunnett p-values for comparisons to the negative control were 0.07 and 0.04.

boxclust(data=ratpup, outcome="weight", treatment="treatment", cluster="litter",
         covariate="sex", xlabel="", ylabel="Pup weight [g]\n", psize=1.5,
         hjitter=0.01, pneg=c(0.07, 0.04), stars=TRUE, white=TRUE)

We plot the data assuming high were the positive control, and the p-values for comparisons to the negative and positive control were 0.009 and 0.12, respectively, and the comparison between both controls yielded a p-values of 0.0004.

boxclust(data=ratpup, outcome="weight", treatment="treatment", cluster="litter",
         covariate="sex", xlabel="", ylabel="Pup weight [g]\n", psize=1.5,
         hjitter=0.01, pneg=0.009, ppos=0.12, pposneg=0.0004, white=TRUE)

PhilipPallmann/toxbox documentation built on May 8, 2019, 1:35 a.m.