Return to extinction_risk_1.R CVS log | Up to [local] / ext |
version 1.1, 2015/05/25 11:45:07 | version 1.2, 2015/05/25 12:21:55 | ||
---|---|---|---|
|
|
||
# 61:115-143, 1991. | # 61:115-143, 1991. | ||
# | # | ||
# library(pracma) | |||
ext1 <- function(dat, t = 100, ne = 1, alpha = 0.05, verbose = FALSE) { | ext1 <- function(dat, t = 100, ne = 1, alpha = 0.05, verbose = FALSE) { | ||
yr <- ts(dat[, 1], start = c(dat[, 1][1])) # Year | yr <- ts(dat[, 1], start = c(dat[, 1][1])) # Year | ||
ps <- ts(dat[, 2], start = c(dat[, 1][1])) # Population size | ps <- ts(dat[, 2], start = c(dat[, 1][1])) # Population size |