Incidental Findings

Cessna 172S Performance Calculations

Predictive modeling of performance characteristics, including takeoff, climb, cruise, and landing. FOR EDUCATIONAL PURPOSES ONLY.

Delaying R Shiny textInput processing

Delay textInput processing until submit button or enter key is pressed, with example Shiny app to fetch METAR data.

Julia computation within Quarto

Quick test of executing a Julia code block from within a Quarto blog post.

Leaflet, GeoTIFF, and Colormaps

Data visualization

Using {leaflet} with a GeoTIFF background layer and resolving difficulties with RGB colormaps

METAR analysis, part 2

Updated to an interactive shinylive client-side web application to analyze historical METAR weather data for flight conditions by month and time of day.

Gestational age calculator

Newborn gestational age calculations in R. Given any set of LMP, EDC, DOB, date, GA, age, DOL, PMA, calculate as many of the others as possible.

PediTools bilirubin webtool via Alfred keyboard shortcut

Creating an Alfred app workflow to trigger the PediTools 2022 AAP hyperbilirubinemia guidelines webtool via zsh script

Gestational age parsing

Helper functions for vectorized parsing of gestational age strings to numeric and back

Future-proofing with PediTools Universal

Preparing PediTools for the future: support for 120+ distinct growth curves and beyond, including bulk calculations, plotting, and custom growth chart generation.

WHO growth standards, expanded

Data visualization

WHO Growth Standard expanded chart parameters – mildly automated data scraping and metadata parsing.

Night flight training

Data visualization

How to spend hours down the rabbit hole of finding precise definitions for night flying and analyzing exactly how accurate the usual rule of thumb is.

Custom checklists

Creating customized, precisely-sized, double-sided, laminated, easy-to-update checklists.

Crafting my own VFR Navlog

Creating a navigation log with heading, time, fuel consumption, and weight & balance calculations, also suitable for printing and folding for use on a kneeboard.

Folding aviation charts and internet archaeology

Retrieving long lost instructions on how to fold aviation charts, both the original Duncan McKillop method as well as a later improvement by Steve Hirlehey

Aviation winds aloft forecast locations

Data visualization

Update the original poor quality 30 kb GIF of aviation wind and temperature aloft forecast locations with R-generated vector graphics

Magic Labyrinth – generating fair (or unfair) game boards

Generate Magic Labyrinth random game boards and analyze fairness for different player start locations.

Plot histogram with overlaid normal curve

Given a vector of values, create a ggplot histogram with overlaid best-fitting normal curve, with prettified caption of numerics

From flight plan to Google Earth video, now via Shiny App

Developed an R Shiny app to facilitate conversions of flight plans to Google Earth video tours.

From flight plan to Google Earth video tour

Using R to convert a .FPL flight plan into a Google Earth video tour

geom_bar plot ordered by value

Create a ggplot geom_bar plot with bars ordered by value, x axis labels rotated, and y axis formatted as percent

METAR analysis

Analyzing historical METAR weather data for flight conditions by month and time of day, including calculation of crosswinds for available runways.

Interpolating missing values

Linear interpolation of missing dependent values based on an independent variable, not permitting extrapolation

First flight

After decades of dreaming of flying, I scheduled and experienced a ‘discovery flight’.

2022 CDC extended BMI-for-age

Data visualization

Extrapolation of the original CDC 2000 BMI-for age growth charts beyond the 97th percentile is a really, really bad idea – graphical demonstration.

Splitting a delimited text column

Convert a delimited text column of substrings in a dataframe into multiple columns specifying presence or absence of that substring

Preterm bilirubin treatment thresholds

Data visualization

Hyperbilirubinemia management recommendations for infants at less than 35 weeks gestation, based on Stanford Premie BiliRecs

Deterministically anonymize identifiers and dates

Code snippets to deterministically anonymize identifiers and dates, while maintaining intra-individual time series relationships, to facilitate sharing examples of clinical…

Unexpected Pandoc failure after RStudio upgrade

I have a fragile pipeline to generate monthly clinical dashboards that fails regularly in fairly typical ways. However, this month it failed in a completely new and…

Data science goals, shorter TLDR version

The last post was way too long. Here’s the short version.

Clinician data scientists – musings in response to an email

I received an email from a PediTools user that led me to think about what it means to be a clinician who also codes and enjoys data science. I really think that more…

Disable Citrix Workspace automatic startup

In a recent update, Citrix decided it would be a good idea to force automatic startup of the Citrix Workspace with every reboot, with no way of disabling it.

PediTools: improving ease of use of the AAP 2022 bilirubin guidelines

PediTools clinical decision support tool for managing neonatal hyperbilirubinemia, including newly developed treatment threshold plots to facilitate clinical care.

Bilirubin treatment thresholds

Data visualization

Comparing bilirubin phototherapy and exchange transfusion thresholds between the AAP 2004 and 2022 guidelines

From Blogdown to Quarto

When this blog was created in 2018, it used the paradigm of R with blogdown to GitHub through Netlify with Hugo to continuous deployment. I set up a custom domain, found a Hu…

Learning Julia with Advent of Code

Advent of Code is a great way to learn a new programming language.

R for Animal Crossing flower genetics

Animal Crossing New Horizons uses Mendelian inheritance of bi-allelic polygenic traits to determine flower colors. I wrote some tools in R to help plan flower crosses.

Finding Minecraft Strongholds with Math

Find a Minecraft stronghold in just three Eye of Ender throws

Using the medianeffect package for dose-effect analysis

Vignette for the medianeffect R package for dose-effect analysis and quantifying synergism and antagonism via the Median Effect Principle.

Cheating at anagrams with Python

How to cheat at anagram word puzzles (or Scrabble) with Python.

RISK battles: Simulating the distribution of outcomes

Playing RISK recently, I was annoyed when I lost a battle despite overwhelming numbers. ‘What are the odds of that?!?’ Which led me to write this simulation.

Tidy Tuesday – NLP on Wine Ratings Dataset

Machine learning

Natural language processing of wine descriptions to predict wine review scores.

Analysis of Pokemon Go Spawn Frequencies to Identify Possible Biomes

Data visualization

Using Pokemon spawns and geocoding to identify biomes.

Generating a union of overlapping intervals

This function takes a set of overlapping intervals and returns the union of the intervals.

Using R in an exploratory data visualization for predicting tree cover

Data visualization

I share here the full code and output for an exploratory data visualization in R Markdown submitted for a contest in January 2017.

Up and running with blogdown, GitHub, Netlify, and tranquilpeak

I am primarily a MacOS user who does most of my analytics and presentation work in R / RStudio / RMarkdown, so it was natural to take advantage of using the R blogdown…

Hello world!

This blog will serve as a place to share data science analyses, professional and personal projects, snippets of useful code, and other various geekery.