thoughtfulbloke’s gists

thoughtfulbloke / Excel_arrangement

Created February 10, 2025 20:34

NZ covid data organised for excel graphing

	"Week_Ending","Border_Cases","Community_Cases","Deaths"
	2020-02-23,0,0,0
	2020-03-01,1,0,0
	2020-03-08,4,2,0
	2020-03-15,3,0,0
	2020-03-22,56,22,0
	2020-03-29,273,180,1
	2020-04-05,160,350,0
	2020-04-12,52,232,4
	2020-04-19,16,81,7

thoughtfulbloke / gist:7aeb5174d6ccc6e09ae170d4de0f16e9

Created January 18, 2025 03:25

proof of concept of batch add info from EXIF as text on photos using magick and exiftoolr in R

	library(magick)
	library(exiftoolr)

	# Hard-coded folder path
	folder <- "~/Desktop/photos"

	# Function to process images
	process_images <- function(folder) {
	# Create the "texted" subfolder if it doesn't exist
	texted_folder <- file.path(folder, "texted")

thoughtfulbloke / phototext.R

Created January 17, 2025 09:07


	library(shiny)
	library(magick)

	ui <- fluidPage(
	titlePanel("Image Upload and Text Overlay"),
	sidebarLayout(
	sidebarPanel(
	fileInput("file1", "Choose an Image File",
	accept = c('image/png', 'image/jpeg')),

thoughtfulbloke / adventcodeday12.R

Created December 13, 2024 03:20

	library(dplyr)
	asLines <- readLines("d12.txt")
	height <- length(asLines)
	width <- nchar(asLines[1])
	chars <- strsplit(paste(asLines, collapse = ""), "")[[1]]

	aCropMatrix <- matrix(chars, nrow = height, byrow = TRUE)
	aFieldMatrix <- matrix(1:length(chars), nrow = height, byrow = TRUE)
	priorFieldMatrix <- matrix(rep(1, length(chars)), nrow = height, byrow = TRUE)
	while (!identical(priorFieldMatrix, aFieldMatrix)) {

thoughtfulbloke / Advent_of_code_2024_day07_part2.R

Created December 7, 2024 21:20

My very tidy solution to day7 in the advent of code

	library(gtools)
	library(stringr)
	library(dplyr)
	library(tidyr)
	d7 <- readLines("day07.txt")
	#assumes you kept the resut of part1 to save work
	d7p1 <- as.numeric(readLines("d7p1result.txt"))
	p2data <- d7[d7p1 == 0]
	p2sums <- as.numeric(gsub(":.*","",p2data))
	p2input <- gsub(".*: ","",p2data)

thoughtfulbloke / gist:4c87c97db569f821eb134eb876e79c99

Created November 24, 2024 18:14

Day 25 2024 #30daymapchallenge code

	library(rvest)
	library(dplyr)
	library(RcppRoll)
	library(tidygeocoder)
	library(maps)
	library(mapdata)
	library(ggplot2)
	library(viridis)
	library(patchwork)
	library(ggtext)

thoughtfulbloke / MortalityRates.csv

Last active December 15, 2024 01:13

NZ Weekly Mortality Rates by Age(<30,30-59,60-79,80+). Date of Death not Registration. Weekly population interpolated from as at Quarter end. source StatsNZ and DIA (via file emailed to me)

Period	Ages	RawDeaths	InterpPop	DeathRateMillion
2000-01-09	30 to 59	74	1569863.62637363	47.1378524585218
2000-01-09	60 to 79	189	495217.472527473	381.650508079589
2000-01-09	80 and over	194	104575.054945055	1855.12692392971
2000-01-09	Under 30	27	1681940.65934066	16.0528850111658
2000-01-16	30 to 59	64	1570169.78021978	40.7599234211741
2000-01-16	60 to 79	196	495324.395604396	395.700275898667
2000-01-16	80 and over	226	104711.208791209	2158.31717166629
2000-01-16	Under 30	14	1681754.50549451	8.32463950847774
2000-01-23	30 to 59	73	1570475.93406593	46.4827243872527

thoughtfulbloke / KnKplus.R

Created March 11, 2024 08:07

KnKplusGibson

	library(dplyr)
	library(tidyr)
	library(lubridate)
	NZstmf <- read.csv("NZL_NPstmfout.csv") \|>
	filter(Sex == "b")
	NZQTR <- read.csv("TotalQuarterlyPopulationAsAt.csv", skip=3,
	col.names = c("YRQR","Pop")) \|>
	filter(!is.na(Pop)) \|>
	separate(YRQR, into=c("YR", "QR"), sep="Q", convert = TRUE) \|>
	mutate(asDate = ceiling_date(ISOdate(YR,QR*3,25),

thoughtfulbloke / AgeRates.csv

Created December 3, 2023 03:49

We can make this file beautiful and searchable if this error is corrected: It looks like row 8 should actually have 20 columns, instead of 18 in line 7.

	"Age-specific death rates by sex, December years (total population) (Annual-Dec)","","","","","","","","","","","","","","","","","","",""
	"","Total Population","","","","","","","","","","","","","","","","","",""
	" ","0-4 Years","5-9 Years","10-14 Years","15-19 Years","20-24 Years","25-29 Years","30-34 Years","35-39 Years","40-44 Years","45-49 Years","50-54 Years","55-59 Years","60-64 Years","65-69 Years","70-74 Years","75-79 Years","80-84 Years","85-89 Years","90 Years and Over"
	1971,4.31,0.39,0.35,0.96,1.06,0.97,1.29,1.81,2.93,4.70,7.50,12.19,19.17,29.34,44.81,71.51,118.74,192.62,280.73
	1972,4.04,0.37,0.33,1.16,1.13,0.89,1.19,1.71,2.65,4.51,7.33,11.72,19.32,29.64,45.89,74.35,118.35,187.90,309.62
	1973,3.99,0.42,0.35,1.22,1.24,0.98,1.17,1.86,2.70,4.76,7.26,12.57,18.57,30.31,45.92,72.05,114.84,190.45,320.67
	1974,3.68,0.38,0.43,1.16,1.06,1.04,1.15,1.76,2.96,4.74,7.51,11.68,19.02,28.10,46.32,70.63,111.84,190.96,314.86
	1975,3.78,0.46,0.34,1.02,1.11,1.01,1.17,1.75,2.64,4.60,7.39,11.07,19.04,27.29,45.13,65.78,106

thoughtfulbloke / scale graph code.R

Created November 12, 2022 04:22

	library(rmapshaper)
	library(sf)
	library(ggplot2)
	library(dplyr)
	library(patchwork)

	# geographic data from koordinates.com
	start_map <- read_sf("new-zealand-2018-estimated-resident-population-statistical-g/new-zealand-2018-estimated-resident-population-statistical-g.shp")
	coast <- read_sf("nz-coastlines-and-islands-polygons-topo-150k/nz-coastlines-and-islands-polygons-topo-150k.shp")
	# simplify the graph to what is visually useful

David Hood thoughtfulbloke