This article puts the hypothesized causal DAG from Part 1 to the test by fitting various Multiple Linear Regression models to ~4,000 large BC wildfires (>100 hectares) and examining whether the atmospheric relationships we drew actually hold up in the data. Before the main analysis, we take a pedagogical detour through Bayesian Multiple Linear Regression using PyMC to demonstrate what happens when weakly informative priors meet a large dataset and why that finding justifies switching to simpler Frequentist methods for the rest of the analysis. The honest result: Our final model explains roughly 11% of variance in fire size, which sounds underwhelming until you consider what that 11% actually represents. Two of our hypothesized mediators don’t survive contact with the data, and our revised DAG is far better for it.

This article makes the case for causal inference as a framework for thinking more deliberately about the models we build and the predictors we choose. Using multiple linear regression as our primary analytical lens, we introduce the core ideas behind causal reasoning, such as the data generating processes, mediators, confounding variables, and spurious correlations, and apply them to a real dataset of large BC wildfires (>100 hectares). The analysis culminates in a hypothesized causal DAG mapping the atmospheric conditions we believe drive wildfire size, which we’ll formally test with regression in Part 2. Spoiler: not every arrow in our DAG will survive contact with the data, and that’s precisely the point.

This notebook demonstrates the complete workflow for building a Bayesian Simple Linear Regression model using PyMC to predict wildfire sizes solely from wind speed data. In it we formalize the model using statistical notation, implement prior predictive simulations to validate our assumptions, and then generate posterior distributions through Markov Chain Monte Carlo (MCMC) sampling. Unfortunately, the analysis revealed that wind speed alone is a weak (quite honestly, a terrible) predictor of fire size. Nonetheless, the value-add in our work was that it provided practical insights about the importance of model diagnostics and the pitfalls of violating assumptions of linearity and heteroscedasticity when working with real-world data!

This notebook introduces the language of statistical model notation used to describe Bayesian models, demonstrating how mathematical notation translates directly into executable PyMC code. Using a real-world example of Vancouver Island Coastal Wolves (a recently discovered subspecies known for displaying distinguished behaviour from other wolf populations such the ability to swim great distances or the significant prevelance of marine organisms in their diet), we build a Bayesian model to estimate gender ratios from limited sample data to emphasize how posterior distributions reveal the full range of plausible outcomes rather than single point estimates.

This notebook introduces the basic idea behind Bayes’ Theorem and highlights some key differences between Bayesian Statistics and the widely taught traditional branch of statistics, commonly known as Frequentist Statistics. Furthermore we challenge the traditional medical testing example used in countless textbooks to introduce Bayesian Inference, arguing that using fixed constants (also known as point estimates) misrepresents the true nature of Bayesian Statistics. Rather than inputting single values in Bayes’ Theorem, we utilize a more faithful Bayesian approach by considering probability distributions of possible outcomes, thus revealing how disease prevalence uncertainty affects diagnostic accuracy.

This notebook demonstrates how to augment wildfire datasets with historical weather data using Python and Google Earth Engine’s ERA5 dataset. Here we transform basic fire records (coordinates and timestamps) into enriched datasets containing comprehensive environmental context including temperature, wind patterns, humidity, and soil conditions which are all critical information for fire risk modeling and analysis. The workflow includes automatic retry logic for handling API timeouts and incremental batch saving to prevent data loss during long processing runs.

This notebook provides a step-by-step guide on wildfire mapping techniques using Python and Google Earth Engine. We explore how to access historical geospatial data from the Sentinel-2 and GOES-16/17 satellites in order to build an interactive map that offers real-time insights on the wildfire perimeters within British Columbia, Canada and beyond.