**Contact Information:** Michael W. Belitz (mbelitz@ufl.edu)

The R package ‘phenesse’ provides functions to calculate phenological estimates for multiple phenological metrics. It can be especially useful when estimating the onset or offset of a phenological event using presence-only data.

We provide example incidental observations from iNaturalist for four species and a small extent of the United States. These data are for the year 2019 up until mid October and are not scored by phenological phases. The four species are Speyeria cybele, Danaus plexippus, Rudbeckia hirta, and Asclepias syriaca.

example iNaturalist data:

Estimate the onset, 10% and 50% of when Speyeria cybele has been observed in 2019 across the entire extent.

```
s_cybele <- subset(inat_examples, scientific_name == "Speyeria cybele")
# calculate onset, note iterations are much lower than desired for quick
# vignette knit. Increase interations for robust calculations
weib_percentile(observations = s_cybele$doy, percentile = 0, iterations = 20)
#> [1] 110.1495
#calculate 1st percentile
weib_percentile(observations = s_cybele$doy, percentile = 0.1, iterations = 20)
#> [1] 155.4506
#calculate 50th percentile
weib_percentile(observations = s_cybele$doy, percentile = 0.5, iterations = 20)
#> [1] 194.3791
```

Estimate the 50% of when Speyeria cybele were observed in 2019 and calculate CI

```
s_cybele <- subset(inat_examples, scientific_name == "Speyeria cybele")
# calculate onset - note ignoring for vignette test
#weib_percentile_ci(observations = s_cybele$doy, iterations = 50,
# percentile = 0.5, bootstraps = 50)
# note warnings of extreme order points used as endpoints is due to the
# low number of bootsraps used. Please use higher number of bootstraps if
# using for analyses.
```

Estimate the 10% and 50% phenometrics and confidence intervals for a quantile estimate of Rudbeckia hirta.

```
r_hirta <- subset(inat_examples, scientific_name == "Rudbeckia hirta")
# calculate 50% quantile and CIs
quantile_ci(observations = r_hirta$doy, percentile = 0.5, bootstraps = 100)
#> Warning in norm.inter(t, adj.alpha): extreme order statistics used as endpoints
#> estimate low_ci high_ci
#> 50% 185.5 173 193
```

Calculate the mean estimate and confidence intervals of the estimate of Rudbeckia hirta.