Instantly share code, notes, and snippets.

# Mike Birdgeneaumikebirdgeneau

• Sort options
Created Sep 19, 2017
Keybase
View keybase.md

### Keybase proof

I hereby claim:

• I am mikebirdgeneau on github.
• I am mikebirdgeneau (https://keybase.io/mikebirdgeneau) on keybase.
• I have a public key ASBq5JoSE_gmzc7c_R7AsWNZv4X0nd6hN2-Y1qZDrNai8wo

To claim this, I am signing this object:

Created Dec 22, 2016
Multivariate Beta Pert Distribution
View mvtbetapert.R
 #' Multivariate Beta PERT distributions #' #' @description Generates random deviates from correlated (modified) pert distributions. #' this is performed by remapping correlated normal distributions to the beta pert #' distributions using quantiles. #' #' @param n Number of observations. If length(n) > 1, the length is taken to be the number required. #' @param min Vector of minima. #' @param mode Vector of modes.
Created Oct 1, 2016
spatial_rotation.R
View Rotation of Arbitrary Spatial Data
 library(data.table) library(sp) library(ggplot2) library(broom) library(gridExtra) set.seed(1) # Set-up Dummy Data (for example) n <-
Created Sep 21, 2016
View 3d_kriging_cgs_cosimulation.R
 library(gstat) library(sp) library(lattice) library(data.table) library(ggplot2) # Create Data Points (Random) n <- 50 data3D <- data.frame(x = runif(n), y = runif(n), z = runif(n), v = rnorm(n)) coordinates(data3D) = ~x+y+z
Created Sep 20, 2016
3D Kriging with Conditional Gaussian Simulation in R
View 3d_kriging_cgs.R
 library(gstat) library(lattice) # Create Data Points (Random) n <- 50 data3D <- data.frame(x = runif(n), y = runif(n), z = runif(n), v = rnorm(n)) coordinates(data3D) = ~x+y+z # Create empty grid to krige range1D <- seq(from = 0, to = 1, length = 20)
Created Apr 6, 2016
When is your birthday if you're born on Feb 29th of a leap year?
View leapYearBirthday.py
 from skyfield.api import load from scipy import optimize from datetime import datetime from dateutil import tz import pandas as pd import numpy as np import matplotlib.pyplot as plt def earthPosition(t1):
Last active Feb 5, 2016
Prediction of Baby Timing & What days of work I'll likely miss...
View BirdgeneauBabyPredictor.R
 library(ggplot2) library(data.table) library(gridExtra) library(zoo) babyTiming<-function(ExpectedDate="2016/02/29",DateofInterest="2016/02/2015",VacationDuration=14){ DateofInterest<-as.Date(DateofInterest) ExpectedDate<-as.Date(ExpectedDate) currentDate<-Sys.Date()
Last active Oct 10, 2015
Federal Spending on Science & Technology
View federal_expenditures_science_technology.R
 library(ggplot2) library(data.table) library(lubridate) library(stringr) options(stringsAsFactors = FALSE) # Data Source: Statistics Canada (Table 358-163) #(http://www5.statcan.gc.ca/cansim/a26?lang=eng&retrLang=eng&id=3580163&&pattern=&stByVal=1&p1=1&p2=31&tabMode=dataTable&csid=)
Last active Oct 10, 2015
Plot GDP per Capita for G7 countries, normalized to when the Harper government took power
View GDP_per_capita_G7.R