knitr::opts_chunk$set( collapse = TRUE, comment = "#>", fig.path = "man/figures/README-", out.width = "100%" )
The goal of flowers is to provide a lodgepole pine flowering phenology dataset in an accessible format. This package contains 15 years of lodgepole pine (Pinus contorta var. latifolia) phenology data for pollen shed and cone receptivity. The data were collected in seed orchards in British Columbia and represents multiple BC provenances.
You can install the development version of flowers from GitHub with:
# install.packages("devtools") devtools::install_github("scisus/flowers")
Read in phenology data sets from the package.
library(flowers) ## basic example code state <- flowers::lodgepole_phenology_state head(state) event <- flowers::lodgepole_phenology_event head(event)
Lodgepole pine in BC generally flowers in May and June.
library(ggplot2) library(dplyr) library(lubridate) flowering <- state %>% filter(State == 2) %>% mutate(Year = as.factor(Year)) ggplot(flowering, aes(x=Site, y=DoY, group=interaction(Site, Year), colour = Site)) + geom_line( alpha=0.5, size=3) + #geom_point(alpha=0.1) + scale_colour_brewer(type="div", palette = 8) + facet_wrap("Sex") + theme_dark() + theme(axis.text.x = element_text(angle = 90, vjust = 0.5, hjust=1)) + geom_abline(intercept = 122, slope = 0) + geom_abline(intercept = 153, slope = 0) + coord_flip() + ggtitle("Day of year trees observed flowering at 7 sites", subtitle = "Lodgepole pine observed between 1997 and 2012. Vertical lines at May and June 1")
Data is from r length(unique(state$Clone))
clones in
r length(unique(state$Orchard))
orchards grown at
r length(unique(state$Site))
BC seed orchard sites and observed between r min(state$Year, na.rm=TRUE)
and r max(state$Year, na.rm=TRUE)
.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.