Statistical analysis:
Initial graphical exploration of refuge use behavior over the course of the study revealed persistent oscillations in the mean refuge use observed each night with approximately a 14-day period. These oscillations in behavior appeared to be negatively correlated with the mean tidal level at Oyster Landing (the collection site of crabs) at the time of observation in the lab. Crabs used the refuge most while it was low tide (when they are generally inactive in the field), indicative of a circatidal rhythm in refuge use. We tested for this influence of the tidal cycle on refuge use behavior, among other factors affecting crab refuge use behavior, in the following analysis. To explore factors influencing crab refuge use behavior (pre-release), we tested the effects of toadfish predation threat, crab carapace width, an interaction between threat and carapace width, and mean tidal level during observation on refuge use with generalized linear mixed models (GLMM, lme4 package in the statistical software R). Female crabs were smaller than male crabs (Welch two-sample t -test: t = -8.267, p <0.001), which confounded crab sex with crab size. Therefore, we tested the effects of these factors on refuge use separately for males and females. Observational block was modeled as a random factor in both GLMM. Because crab refuge use (the response variable) was proportional, we modeled this behavior using a binomial distribution and logit link (Bolker et al. 2009).
After recapturing a portion of these original crabs (108 crabs recaptured), we calculated the repeatability of their refuge use behavior using pre-release and post-recapture behavioral measurements. Repeatability (r) is defined as the proportion of the total variation that occurs within individuals as opposed to between individuals and is calculated as r =sA 2 / (s 2+sA 2 ), where sA 2 is the among-individual variance and s 2 is the within-individual variance (Bell et al. 2009; Nakagawa and Schielzeth 2010). Thus, repeatability provides a metric of the amount of behavioral variation between relative to within individuals, where a higher repeatability value indicates a higher level of individual behavioral consistency between measurements. Again, due to the proportional behavioral measure (refuge use), we used GLMM-based repeatability estimation (rptR package in R, Nakagawa and Schielzeth 2010) with a binomial distribution and logit link. Repeatability was calculated separately for crab refuge use in the absence and presence of the toadfish cue, and confidence intervals (95 %) and statistical significance (p-values) were estimated using parametric bootstrapping with 1,000 resamplings.
Next, we explored factors driving change in refuge use behavior (i.e., deviation from perfect repeatability) of recaptured crabs after time in the field. We calculated behavioral change by subtracting the value of the first behavioral measurement (pre-release) from the second behavioral measurement (post-recapture). Behavioral change was logtransformed to meet assumptions of linear regression. To test for a predominant direction in behavioral change, we first tested whether behavioral change was significantly different than zero in the absence and presence of toadfish predation threat using one-sample t -tests. We then used general linear models testing the fixed effects of duration (days) in the field, crab carapace width, and crab sex on individual behavioral change of recaptured crabs. Two separate linear models were used to test the effects of these factors in the absence and presence of toadfish predation threat.
As previously mentioned, crab refuge use behavior oscillated with a circatidal rhythm over the course of the study. This means that if a recaptured crab was originally observed (pre-release) during one tidal height and observed for a second time (post-recapture) at a different tidal height, then behavioral change would be generated. We tested for this tidal influence on behavioral change as follows. We first subtracted for each individual crab the tidal height when the pre-release observation was made from the tidal height when the postrecapture observation was made. We then used the absolute value of this difference as a factor (termed ‘tidal influence’) in general linear models testing the effects of duration in the field, carapace width, sex, and the tidal influence on the absolute value of behavioral change in the absence and presence of predation threat. This analysis allowed us to explore the relative influences of these factors on the overall magnitude of behavioral change.
Lastly, we tested for the differential recapture of crabs with low versus high refuge use (i.e., a sampling bias).We did this by comparing the recapture rate of crabs from the lower and upper quartiles of refuge use behavior using Fisher's exact tests. We conducted this analysis separately for crabs with refuge use measured in the absence and presence of toadfish predation threat
BCO-DMO Processing:
- added conventional header with dataset name, PI name, version date, reference information
- renamed parameters to BCO-DMO standard
- reduced number of digits for refuge use and activity level from 9 to 2 places after decimal