Become an expert in R — Interactive courses, Cheat Sheets, certificates and more!
Get Started for Free

ossification

Teratogenic effects of phenytoin and trichloropropene oxide


Description

The data come from a 2x2 factorial design with 81 pregnant mice. In the experiment each pregnant mouse was randomly allocated to an control group and three treated groups, which received daily, by gastric gavages, 60 mg/kg of phenytoin, 100 mg/kg of trichloropropene oxide, or 60 mg/kg phenytoin and 100 mg/kg of trichloropropene oxide. On day 18 of gestation, fetuses were recovered, stained, and cleared. Then, by visual inspection, the presence or absence of ossification was determined for the different joints of the right and left forepaws. The purpose of the experiment was to investigate the synergy of phenytoin and trichloropropene oxide to produce ossification at the phalanges, that is, teratogenic effects. See Morel and Nagaraj (2012, page 103).

Usage

data(ossification)

Format

A data frame with 81 rows and 4 variables:

fetuses

a numeric vector giving the number of fetuses showing ossification on the left middle third phalanx.

litter

a numeric vector giving the litter size.

pht

a factor giving the dose (mg/kg) of phenytoin: "0 mg/kg" or "60 mg/kg".

tcpo

a factor giving the dose (mg/kg) of trichloropropene oxide: "0 mg/kg" or "100 mg/kg".

References

Morel J.G. and Neerchal N.K. (1997) Clustered binary logistic regression in teratology data using a finite mixture distribution. Statistics in Medicine 16, 2843-2853.

Morel J.G. and Nagaraj N.K. (2012) Overdispersion Models in SAS. SAS Institute Inc., Cary, North Carolina, USA.

Examples

boxplot(100*fetuses/litter ~ pht, data=subset(ossification, tcpo=="0 mg/kg"),
        at=c(1:2)-0.2, col="yellow", boxwex=0.25, outline=FALSE, xaxt="n",
        xlab="Dose of PHT", ylab="% of fetuses showing ossification")
boxplot(100*fetuses/litter ~ pht, data=subset(ossification, tcpo=="100 mg/kg"),
        add=TRUE, at=c(1:2)+0.2, col="blue", boxwex=0.25, outline=FALSE, xaxt="n")
axis(1, at=1:2, labels=levels(ossification$pht))
legend(0.2, 45, legend=c("0 mg/kg","100 mg/kg"), y.intersp=0.1, x.intersp=0.3, pch=15,
       col=c("yellow","blue"), bty="n", cex=0.9, title="Dose of TCPO", title.adj=0.2)

glmtoolbox

Set of Tools to Data Analysis using Generalized Linear Models

v0.1.0
GPL-2 | GPL-3
Authors
Luis Hernando Vanegas [aut, cre], Luz Marina Rondón [aut], Gilberto A. Paula [aut]
Initial release

We don't support your browser anymore

Please choose more modern alternatives, such as Google Chrome or Mozilla Firefox.