pbaylis

program define gen_relative_humidity
	* Calculate relative humidity given dry bulb temperature and dew point
	syntax, DRYbulbtemp(varname) DEWpoint(varname) [GENerate(string)]
	* Inputs in F
	local T `drybulbtemp'
	local TD `dewpoint'

	if "`generate'" == "" {
		local generate "rel_humidity"
	}

pbaylis

program define gen_heat_index
	* Calculat heat index
	* See https://en.wikipedia.org/wiki/Heat_index. Note that there is no exact heat index formula
	* since it is defined by the table. There are approximations, but this code uses the table itself.
	* The only modification is that I cap at 140.
	syntax, DRYbulbtemp(varname) RELativehumidity(varname) [GENerate(string)]
	
	if "`generate'" == "" {
		local generate "heat_index"
	}

pbaylis

program define makebins, rclass
    syntax varname, [PREfix(name) min(real 0) max(real 100) step(real 5) SEQuential DROPbin(real 70)]
    local var `varlist'

    if mod(`max' - `min',`step') != 0 {
        di as error "Warning: Step size `step' does not divide evenly into range. Top bin may extend beyond `max'."
    }

    * Declare prefix
    if "`prefix'" != "" {

pbaylis

#!/bin/bash
DATE=`date "+%A, %b %d %Y"`
YESTERDAY=`date -v-1d "+%A, %b %d %Y"`
open "bear://x-callback-url/create?title=$DATE&text=## Beginning of day thoughts

## Goals
  -
## End of day thoughts
[[$YESTERDAY]]&tags=dailylog"

pbaylis

## Logging into server
ssh <server-address>
# Login
tmux
# CTRL-b d to detach tmux
exit

## When returning to server
tmux attach

pbaylis

x <- c(5, 15, 2015)
sprintf(sprintf("%03d", x))

pbaylis

points.shp <- points.shp[apply(sf::st_intersects(points.shp, polys), 1, any),]

pbaylis

multipetal <- function(df, n) {
    varname <- paste("petal", n , sep=".")
    mutate(df, !!varname := Petal.Width * n)
}

pbaylis

Reduce(function(x, y) merge(x, y, all=TRUE), list(df1, df2, df3))
# OR #
Reduce(function(...) merge(..., all=TRUE), list(df1, df2, df3))

pbaylis

# Create sensitivity curve of coefficient estimates
library(tidyverse)
library(cowplot)
library(fastDummies)

# Setup ----
rm(list = ls())
theme_set(theme_cowplot())
set.seed(42)