This data set describes the habitat preference of two species of
lizards, Anolis grahami and A. opalinus, on the
island of Jamaica and is originally from Schoener (1970). McCullagh and
Nelder (1989) and Burnham and Anderson (2002) reanalyzed the data. Note
that a typo occurs in table 3.11 of Burnham and Anderson (2002).
Usage
data(lizards)
Format
A data frame with 48 rows and 6 variables.
Insolation
position of perch, either shaded or
sunny.
Diameter
diameter of the perch, either < 2 in
or >= 2 in.
Height
perch height, either < 5 or
>= 5.
Time
time of day, either morning,
midday, or afternoon.
Species
species observed, either grahami or
opalinus.
Counts
number of individuals observed.
Details
Burnham and Anderson (2002, p. 137) use this data set originally from
Schoener (1970) to illustrate model selection for log-linear models.
Source
Burnham, K. P., Anderson, D. R. (2002) Model Selection and
Multimodel Inference: a practical information-theoretic
approach. Second edition. Springer: New York.
McCullagh, P., Nelder, J. A. (1989) Generalized Linear
Models. Second edition. Chapman and Hall: New York.
Schoener, T. W. (1970) Nonsynchronous spatial overlap of lizards in
patchy habitats. Ecology51, 408–418.
Examples
data(lizards)
## Not run:
##log-linear model as in Burnham and Anderson 2002, p. 137
##main effects
m1 <- glm(Counts ~ Insolation + Diameter + Height + Time + Species,
family = poisson, data = lizards)
##main effects and all second order interactions = base
m2 <- glm(Counts ~ Insolation + Diameter + Height + Time + Species +
Insolation:Diameter + Insolation:Height + Insolation:Time +
Insolation:Species + Diameter:Height + Diameter:Time +
Diameter:Species + Height:Time + Height:Species +
Time:Species, family = poisson, data = lizards)
##base - DT
m3 <- glm(Counts ~ Insolation + Diameter + Height + Time + Species +
Insolation:Diameter + Insolation:Height + Insolation:Time +
Insolation:Species + Diameter:Height + Diameter:Species +
Height:Time + Height:Species + Time:Species,
family = poisson, data = lizards)
##base + HDI + HDT + HDS
m4 <- glm(Counts ~ Insolation + Diameter + Height + Time + Species +
Insolation:Diameter + Insolation:Height + Insolation:Time +
Insolation:Species + Diameter:Height + Diameter:Time +
Diameter:Species + Height:Time + Height:Species +
Time:Species + Height:Diameter:Insolation +
Height:Diameter:Time + Height:Diameter:Species,
family = poisson, data = lizards)
##base + HDI + HDS + HIT + HIS + HTS + ITS
m5 <- glm(Counts ~ Insolation + Diameter + Height + Time + Species +
Insolation:Diameter + Insolation:Height + Insolation:Time +
Insolation:Species + Diameter:Height + Diameter:Time +
Diameter:Species + Height:Time + Height:Species +
Time:Species + Height:Diameter:Insolation +
Height:Diameter:Species + Height:Insolation:Time +
Height:Insolation:Species + Height:Time:Species +
Insolation:Time:Species, family = poisson, data = lizards)
##base + HIT + HIS + HTS + ITS
m6 <- glm(Counts ~ Insolation + Diameter + Height + Time + Species +
Insolation:Diameter + Insolation:Height + Insolation:Time +
Insolation:Species + Diameter:Height + Diameter:Time +
Diameter:Species + Height:Time + Height:Species +
Time:Species + Height:Insolation:Time +
Height:Insolation:Species + Height:Time:Species +
Insolation:Time:Species, family = poisson, data = lizards)
##base + HIS + HTS + ITS
m7 <- glm(Counts ~ Insolation + Diameter + Height + Time + Species +
Insolation:Diameter + Insolation:Height + Insolation:Time +
Insolation:Species + Diameter:Height + Diameter:Time +
Diameter:Species + Height:Time + Height:Species +
Time:Species + Height:Insolation:Species +
Height:Time:Species + Insolation:Time:Species,
family = poisson, data = lizards)
##base + HIT + HIS + HTS + ITS - DT
m8 <- glm(Counts ~ Insolation + Diameter + Height + Time + Species +
Insolation:Diameter + Insolation:Height + Insolation:Time +
Insolation:Species + Diameter:Height + Diameter:Species +
Height:Time + Height:Species + Time:Species +
Height:Insolation:Time + Height:Insolation:Species +
Height:Time:Species + Insolation:Time:Species,
family = poisson, data = lizards)
##base + HIT + HIS + ITS - DT
m9 <- glm(Counts ~ Insolation + Diameter + Height + Time + Species +
Insolation:Diameter + Insolation:Height + Insolation:Time +
Insolation:Species + Diameter:Height + Diameter:Species +
Height:Time + Height:Species + Time:Species +
Height:Insolation:Time + Height:Insolation:Species +
Insolation:Time:Species,
family = poisson, data = lizards)
##base + HIT + HIS - DT
m10 <- glm(Counts ~ Insolation + Diameter + Height + Time + Species +
Insolation:Diameter + Insolation:Height + Insolation:Time +
Insolation:Species + Diameter:Height + Diameter:Species +
Height:Time + Height:Species + Time:Species +
Height:Insolation:Time + Height:Insolation:Species,
family = poisson, data = lizards)
##set up in list
Cands <- list(m1, m2, m3, m4, m5, m6, m7, m8, m9, m10)
Modnames <- paste("m", 1:length(Cands), sep = "")
##model selection
library(AICcmodavg)
aictab(Cands, Modnames)
## End(Not run)
Results
R version 3.3.1 (2016-06-21) -- "Bug in Your Hair"
Copyright (C) 2016 The R Foundation for Statistical Computing
Platform: x86_64-pc-linux-gnu (64-bit)
R is free software and comes with ABSOLUTELY NO WARRANTY.
You are welcome to redistribute it under certain conditions.
Type 'license()' or 'licence()' for distribution details.
R is a collaborative project with many contributors.
Type 'contributors()' for more information and
'citation()' on how to cite R or R packages in publications.
Type 'demo()' for some demos, 'help()' for on-line help, or
'help.start()' for an HTML browser interface to help.
Type 'q()' to quit R.
> library(AICcmodavg)
> png(filename="/home/ddbj/snapshot/RGM3/R_CC/result/AICcmodavg/lizards.Rd_%03d_medium.png", width=480, height=480)
> ### Name: lizards
> ### Title: Habitat Preference of Lizards
> ### Aliases: lizards
> ### Keywords: datasets
>
> ### ** Examples
>
> data(lizards)
> ## Not run:
> ##D ##log-linear model as in Burnham and Anderson 2002, p. 137
> ##D ##main effects
> ##D m1 <- glm(Counts ~ Insolation + Diameter + Height + Time + Species,
> ##D family = poisson, data = lizards)
> ##D
> ##D
> ##D ##main effects and all second order interactions = base
> ##D m2 <- glm(Counts ~ Insolation + Diameter + Height + Time + Species +
> ##D Insolation:Diameter + Insolation:Height + Insolation:Time +
> ##D Insolation:Species + Diameter:Height + Diameter:Time +
> ##D Diameter:Species + Height:Time + Height:Species +
> ##D Time:Species, family = poisson, data = lizards)
> ##D
> ##D ##base - DT
> ##D m3 <- glm(Counts ~ Insolation + Diameter + Height + Time + Species +
> ##D Insolation:Diameter + Insolation:Height + Insolation:Time +
> ##D Insolation:Species + Diameter:Height + Diameter:Species +
> ##D Height:Time + Height:Species + Time:Species,
> ##D family = poisson, data = lizards)
> ##D
> ##D ##base + HDI + HDT + HDS
> ##D m4 <- glm(Counts ~ Insolation + Diameter + Height + Time + Species +
> ##D Insolation:Diameter + Insolation:Height + Insolation:Time +
> ##D Insolation:Species + Diameter:Height + Diameter:Time +
> ##D Diameter:Species + Height:Time + Height:Species +
> ##D Time:Species + Height:Diameter:Insolation +
> ##D Height:Diameter:Time + Height:Diameter:Species,
> ##D family = poisson, data = lizards)
> ##D
> ##D ##base + HDI + HDS + HIT + HIS + HTS + ITS
> ##D m5 <- glm(Counts ~ Insolation + Diameter + Height + Time + Species +
> ##D Insolation:Diameter + Insolation:Height + Insolation:Time +
> ##D Insolation:Species + Diameter:Height + Diameter:Time +
> ##D Diameter:Species + Height:Time + Height:Species +
> ##D Time:Species + Height:Diameter:Insolation +
> ##D Height:Diameter:Species + Height:Insolation:Time +
> ##D Height:Insolation:Species + Height:Time:Species +
> ##D Insolation:Time:Species, family = poisson, data = lizards)
> ##D
> ##D ##base + HIT + HIS + HTS + ITS
> ##D m6 <- glm(Counts ~ Insolation + Diameter + Height + Time + Species +
> ##D Insolation:Diameter + Insolation:Height + Insolation:Time +
> ##D Insolation:Species + Diameter:Height + Diameter:Time +
> ##D Diameter:Species + Height:Time + Height:Species +
> ##D Time:Species + Height:Insolation:Time +
> ##D Height:Insolation:Species + Height:Time:Species +
> ##D Insolation:Time:Species, family = poisson, data = lizards)
> ##D
> ##D ##base + HIS + HTS + ITS
> ##D m7 <- glm(Counts ~ Insolation + Diameter + Height + Time + Species +
> ##D Insolation:Diameter + Insolation:Height + Insolation:Time +
> ##D Insolation:Species + Diameter:Height + Diameter:Time +
> ##D Diameter:Species + Height:Time + Height:Species +
> ##D Time:Species + Height:Insolation:Species +
> ##D Height:Time:Species + Insolation:Time:Species,
> ##D family = poisson, data = lizards)
> ##D
> ##D ##base + HIT + HIS + HTS + ITS - DT
> ##D m8 <- glm(Counts ~ Insolation + Diameter + Height + Time + Species +
> ##D Insolation:Diameter + Insolation:Height + Insolation:Time +
> ##D Insolation:Species + Diameter:Height + Diameter:Species +
> ##D Height:Time + Height:Species + Time:Species +
> ##D Height:Insolation:Time + Height:Insolation:Species +
> ##D Height:Time:Species + Insolation:Time:Species,
> ##D family = poisson, data = lizards)
> ##D
> ##D ##base + HIT + HIS + ITS - DT
> ##D m9 <- glm(Counts ~ Insolation + Diameter + Height + Time + Species +
> ##D Insolation:Diameter + Insolation:Height + Insolation:Time +
> ##D Insolation:Species + Diameter:Height + Diameter:Species +
> ##D Height:Time + Height:Species + Time:Species +
> ##D Height:Insolation:Time + Height:Insolation:Species +
> ##D Insolation:Time:Species,
> ##D family = poisson, data = lizards)
> ##D
> ##D ##base + HIT + HIS - DT
> ##D m10 <- glm(Counts ~ Insolation + Diameter + Height + Time + Species +
> ##D Insolation:Diameter + Insolation:Height + Insolation:Time +
> ##D Insolation:Species + Diameter:Height + Diameter:Species +
> ##D Height:Time + Height:Species + Time:Species +
> ##D Height:Insolation:Time + Height:Insolation:Species,
> ##D family = poisson, data = lizards)
> ##D
> ##D ##set up in list
> ##D Cands <- list(m1, m2, m3, m4, m5, m6, m7, m8, m9, m10)
> ##D Modnames <- paste("m", 1:length(Cands), sep = "")
> ##D
> ##D ##model selection
> ##D library(AICcmodavg)
> ##D aictab(Cands, Modnames)
> ## End(Not run)
>
>
>
>
>
> dev.off()
null device
1
>