spaMM (version 2.1.6)

salamander: Salamander mating data

Description

Data from a salamander mating experiment discussed by McCullagh and Nelder (1989, Ch. 14). Twenty males and twenty females from two populations (Rough Butt and Whiteside) were each paired with 6 individuals from their own or from the other population. The experiments were later published by Arnold et al. (1996).

Usage

data(salamander)

Format

The data frame includes 360 observations on the following variables:

Female

Index of the female;

Male

Index of the male;

Mate

Whether the pair successfully mated or not;

TypeF

Population of origin of female;

TypeM

Population of origin of male;

Cross

Interaction term between TypeF and TypeM;

Season

A factor with levels Summer and Fall;

Experiment

Index of experiment

References

Arnold, S.J., Verrell, P.A., and Tilley S.G. (1996) The evolution of asymmetry in sexual isolation: a model and a test case. Evolution 50, 1024-1033.

McCullagh, P. and Nelder, J.A. (1989). Generalized Linear Models, 2nd edition. London: Chapman & Hall.

Molas, M., Lesaffre, E. (2011) Hierarchical Generalized Linear Models: The R Package HGLMMM. Journal of Statistical Software 39, 1-20.

Examples

# NOT RUN {
data(salamander)

# }
# NOT RUN {
<!-- % in test-salamander.R -->
# }
# NOT RUN {
 # < 1 s.
 HLfit(cbind(Mate,1-Mate)~TypeF+TypeM+TypeF*TypeM+(1|Female)+(1|Male),
    family=binomial(),data=salamander,HLmethod="ML")
# }
# NOT RUN {
# In this example fitme() is slower
if (spaMM.getOption("example_maxtime")>1.2) {
 fitme(cbind(Mate,1-Mate)~TypeF+TypeM+TypeF*TypeM+(1|Female)+(1|Male),
    family=binomial(),data=salamander)
}
# }