FlexParamCurve (version 1.5-5)

penguin.data: Field data on growth of little penguins \(Eudyptula minor\)

Description

The penguin.data data frame has 2244 rows and 11 columns of records of the measured masses for little penguin chicks between 13 and 74 days of age collected at Philip Island, Victoria, Australia in 2000 and 2002 (see Chiaradia & Nisbet 2006).

Usage

penguin.data

Format

This object of class c("nfnGroupedData", "nfGroupedData", "groupedData", "data.frame") containing the following columns:

site

Three character factor for the site (only one site in dataset).

year

A factor specifying the year of measurement.

bandid

an ordered factor indicating unique id of each individual: the union of the laying date of the nest relative to the colony and the band combination

siteyear

A factor specifying levels of year for different sites (only one site in dataset).

weight

a numeric vector of chick masses (g).

ckage

a numeric vector of chick ages (days).

Jdate

a numeric vector of first egg-laying date of the nest(days), relative to the mean laying date for all nests in that year.

nest

A factor of unique codes that identify each nest.

ck

A factor of hatching order for each chick (A = first hatched, B = second hatched).

outcome

A factor of codes for fate of each chick (F = fledged; only fledged chicks included).

clutch

A factor of size of clutch/brood that each chick comes from (either 1- or 2-chick brood).

Details

Data were collected as outlined in Chiaradia & Nisbet (2006). Penguin chicks are generally considered to exhibit a double-Gompertz growth form.Please contact Andre Chiaradia ([e-mail]) for use in collaborations.

Examples

require(stats); require(graphics)
#view data
penguin.data
modpar(penguin.data$ckage, penguin.data$weight, pn.options = "myoptions")
plot(weight ~ ckage, data = penguin.data, subset = bandid == penguin.data$bandid[1],
     xlab = "Chick age (day)", las = 1,
     ylab = "Chick mass (g)",
     main = "penguin.data and fitted curve (Chick #307 only)")
fm1 <- nls(weight ~ SSposnegRichards(ckage,Asym=Asym,K=K,Infl=Infl, RAsym=RAsym,
           modno=31, pn.options= "myoptions"),
           data = penguin.data, subset = bandid == penguin.data$bandid[1])
ckage <- seq(0, 74, length.out = 101)
lines(ckage, predict(fm1, list(ckage = ckage)))

Run the code above in your browser using DataCamp Workspace