# Data: Counts of Lucy's warblers on 36 plots established along a gradient of shrub invasion # Plots were surveyed 4 times each, all in a single year # Site-level (plot) covariates, all of which have been natural-log transformed: # LBareCov is percent cover of plot that is bare gound # LForbCov " " covered with forbs # LLitterCov " " covered with litter # LGrassCov " " covered with grass # LShrubCov " " covered with shrubs (<= 2 m height) # LTreeCov " " covered with trees (> 2 m height) # Survey covariate: # week is week of the year (numbered 1 for 1-6 January, and so on) library(unmarked) # Import data and create unmarkedFrame bird.data <- read.csv("http://ag.arizona.edu/~steidl/LUWA_count.csv", header=TRUE) luwa.umf <- unmarkedFramePCount( y=bird.data[,1:4], siteCovs=bird.data[,c("LBareCov", "LForbCov", "LLitterCov", "LGrassCov", "LShrubCov", "LTreeCov")], obsCovs=list(week=bird.data[,5:8]))