stremo (version 0.2)

harnik: Effects of biological factors on extinction risk in fossil bivalves.

Description

Raw data and model syntaxes for the paper by Harnik (2011) on the direct and indirect effects of biological factors on extinction risk in fossil bivalves.

Usage

data(harnik) harnik.model.A harnik.model.B harnik.model.B

Format

A data frame with 108 observations on 8 variables and three character vectors containing model syntaxes.

Source

Article: http://www.pnas.org/cgi/doi/10.1073/pnas.1100572108 Supplementary information: http://www.pnas.org/content/suppl/2011/07/25/1100572108.DCSupplemental

References

Harnik P. 2011. Direct and indirect effects of biological factors on extinction risk in fossil bivalves. PNAS early view. doi: 10.1073/pnas.1100572108.

Examples

## Not run: 
# data(harnik)
# harnik.t <- harnik
# 
# harnik.t$abundance <- asin(sqrt(harnik.t$abundance))
# harnik.t$range.extent <- asin(sqrt(harnik.t$range.extent))
# 
# harnik.t$abundance <- log(harnik.t$abundance)
# harnik.t$body.size <- log(harnik.t$body.size)
# 
# harnik.t$abundance <- scale(harnik.t$abundance, center = TRUE, scale = TRUE)
# harnik.t$range.extent <- scale(harnik.t$range.extent, center = TRUE, scale =
#                                TRUE)
# harnik.t$body.size <- scale(harnik.t$body.size, center = TRUE, scale = TRUE)
# 
# harnik.model.A.1.fit <- sem(harnik.model.A[[1]], data = harnik.t)
# harnik.model.A.2.fit <- sem(harnik.model.A[[2]], data = harnik.t)
# harnik.model.A.3.fit <- sem(harnik.model.A[[3]], data = harnik.t)
# harnik.model.B.fit <- sem(harnik.model.B, data = harnik.t)
# harnik.model.C.fit <- sem(harnik.model.C, data = harnik.t, group =
#         "superfamily", 
#         group.equal = "regressions", 
#         group.partial = "duration ~ body.size") 
# 
# summary(harnik.model.A.fit)
# summary(harnik.model.B.fit)
# summary(harnik.model.C.fit)
# ## End(Not run)