ext/extinction_risk_1.R - diff

Return to extinction_risk_1.R CVS log

Up to [local] / ext

Diff for /ext/extinction_risk_1.R between version 1.10 and 1.11

version 1.10, 2015/08/06 11:05:29

version 1.11, 2015/08/06 11:08:24

Line 245 print.ext1 <- function(obj, digits = 5) {

# Examples

# Yellowstone grizzly bears (from Dennis et al., 1991)

dat <- data.frame(Year = c(1959, 1960, 1961, 1962, 1963, 1964, 1965, 1966, 1967, 1968, 1969, 1970, 1971, 1972, 1973, 1974, 1975, 1976, 1977, 1978, 1979, 1980, 1981, 1982, 1983, 1984, 1985, 1986, 1987),

# dat <- data.frame(Year = c(1959, 1960, 1961, 1962, 1963, 1964, 1965, 1966, 1967, 1968, 1969, 1970, 1971, 1972, 1973, 1974, 1975, 1976, 1977, 1978, 1979, 1980, 1981, 1982, 1983, 1984, 1985, 1986, 1987),

Population = c(44, 47, 46, 44, 46, 45, 46, 40, 39, 39, 42, 44, 41, 40, 33, 36, 34, 39, 35, 34, 38, 36, 37, 41, 39, 51, 47, 57, 47))

# Population = c(44, 47, 46, 44, 46, 45, 46, 40, 39, 39, 42, 44, 41, 40, 33, 36, 34, 39, 35, 34, 38, 36, 37, 41, 39, 51, 47, 57, 47))

# The probability of extinction (of decline to population size 1) within 100 years

# ext1(dat, t = 100)

# The probability of decline to 10 individuals within 100 years