In brouwern/DRmencia: Data and code to replicate analyses for Latta, Brouwer, et al. 2018. PeerJ.

knitr::opts_chunk$set(echo = TRUE)

Libraries

library(reshape2)

Load cleaned data

filename <- "aceit_cleaned.csv"
filename <- here::here("data-raw",
                  "data-raw-aceitillar_96-02",
                  filename)

aceit <- read.csv(file = filename)

head(aceit)

Subset focal spp

aceit.focals <- aceit[which(aceit$status.focals %in% c("res","mig")),]

Put in order

focal.mig <- c("OVEN","BAWW","COYE","AMRE","CMWA",
               "BTBW","PAWA","PRAW")

focal.res <- c("HILC" # formerly HLCU
,"STOF"
,"RLTH"
,"NOMO"
,"GTGT"   # formerly GRWA
,"BANA"
,"BCPT"
,"YFGR"
,"BFGR"
,"GABU")

focals <- c(focal.mig,focal.res)

aceit.focals$spp.code <- factor(aceit.focals$spp.code,
                             levels = focals)

Cast

Calcualte sample size

aceit.N <- dcast(data = aceit.focals,
                   formula = spp.code + status.focals ~ site,
                   value.var = "band",
                   fun.aggregate = length)

Calcualte sample size

aceit.sex <- dcast(data = aceit.focals,
                   formula = spp.code + status.focals ~ sex,
                   value.var = "band",
                   fun.aggregate = length)

brouwern/DRmencia documentation built on May 6, 2019, 12:24 p.m.