zuidhof.broiler: Daily weight, feed, egg measurements for a broiler In agridat: Agricultural Datasets

Description

Daily weight, feed, egg measurements for a broiler

Format

A data frame with 59 observations on the following 6 variables.

`bw`

Body weight, grams

`targetbw`

Target body weight, grams

`adfi`

Average daily feed intake, grams

`adg`

Average daily gain, grams

`eggwt`

Egg weight, grams

`age`

Age, days

Details

Using graphs like the one in the examples section, the authors discovered that a drop in body weight commonly occurs around the time of first egg production.

Source

Martin J. Zuidhof and Robert A. Renema and Frank E. Robinson, (2008). Understanding Multiple, Repeated Animal Measurements with the Help of PROC GPLOT. SAS Global Forum 2008, Paper 250-2008.

Used with permission of Martin Zuidhof.

Examples

 ``` 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``` ```data(zuidhof.broiler) dat <- zuidhof.broiler dat <- transform(dat, age=age/7) # Change days into weeks # Reproducing figure 1 of Zuidhof et al. # Plot using left axis op <- par(mar=c(5,4,4,4)) plot(bw~age, dat, xlab="Age (weeks)", ylab="Bodyweight (g)", main="zuidhof.broiler", xlim=c(20,32), ylim=c(0,4000), pch=20) lines(targetbw~age, subset(dat, !is.na(targetbw)), col="black") # Now plot using the right axis par(new=TRUE) plot(adfi~age, subset(dat, !is.na(adfi)), xlab="", ylab="", xlim=c(20,32), xaxt="n",yaxt="n", ylim=c(-50,175), type="s", lty=2) axis(4, at=c(-50,-25,0,25,50,75,100,125,150,175), col="red", col.axis="red") mtext("Weight (g)", side=4, line=2, col="red") lines(adg~age, subset(dat, !is.na(adg)), col="red", type="s", lty=1, lwd=2) abline(h=c(0,52), col="red") with(dat, segments(age, 0, age, eggwt, col="red")) legend(20, -40, c("Body weight", "Target BW", "Feed/day", "Gain/day", "Egg wt"), bty="n", cex=.5, ncol=5, col=c("black","black","red","red","red"), lty=c(-1,1,2,1,1), lwd=c(1,1,1,2,1), pch=c(20,-1,-1,-1,-1)) par(op) ```

agridat documentation built on May 2, 2019, 4:01 p.m.