data(SwissMarbledWhite)
str(SwissMarbledWhite)
# Data preparation for AHM2 section 4.13
y <- as.matrix(SwissMarbledWhite[,14:24]) # Grab detection/nondetection data
DATE <- as.matrix(SwissMarbledWhite[,3:13]) # Grab survey dates
for(t in 1:11) { # Mean-impute date (but don't transform)
DATE[is.na(DATE[,t]),t] <- mean(DATE[,t], na.rm=TRUE)
}
year <- SwissMarbledWhite$year
nsites <- length(unique(SwissMarbledWhite$site))
nyears <- length(unique(SwissMarbledWhite$year))
nsurveys <- ncol(y)
nobs <- nrow(y)
Run the code above in your browser using DataLab