Skip to content

Instantly share code, notes, and snippets.

@aaronberdanier
aaronberdanier / popsim.r
Created December 10, 2012 17:55
Simulation
time <- 100
N <- Nhat <- dN <- K <- rep(NA,time)
P <- rlnorm(100,7,0.7)
r <- 0.05
Kmax <- 500
alpha <- mean(P)
N[1] <- 100
Nhat[1] <- rpois(1,N[1])
K[1] <- NA