Skip to content

Instantly share code, notes, and snippets.

View davipatti's full-sized avatar

David Pattinson davipatti

  • UW-Madison
  • New York, USA
View GitHub Profile
@davipatti
davipatti / matplotlib-error.txt
Created January 13, 2023 16:47
Matplotlib error - fixed by pip installing reportlab
pattinson@MacBook-Pro-82 ~/r/MapDeduce (master)> python (limix)
Python 3.8.15 (default, Nov 24 2022, 09:04:07)
[Clang 14.0.6 ] :: Anaconda, Inc. on darwin
Type "help", "copyright", "credits" or "license" for more information.
>>> import matplotlib.pyplot as plt
Traceback (most recent call last):
File "<stdin>", line 1, in <module>
File "/Users/pattinson/opt/anaconda3/envs/limix/lib/python3.8/site-packages/matplotlib/pyplot.py", line 32, in <module>
import matplotlib.colorbar
File "/Users/pattinson/opt/anaconda3/envs/limix/lib/python3.8/site-packages/matplotlib/colorbar.py", line 32, in <module>
@davipatti
davipatti / least-absolute-deviations.ipynb
Created January 4, 2023 16:56
Demonstrating least absolute deviation regression #numpy #regression
Sorry, something went wrong. Reload?
Sorry, we cannot display this file.
Sorry, this file is invalid so it cannot be displayed.
@davipatti
davipatti / model-index-mixing.py
Created November 21, 2022 20:31
[Model index not mixing] pymc example code to illustrate model index parameter not mixing #pymc #bayesian
#!/usr/bin/env python
import pymc as pm
import numpy as np
import arviz as az
import matplotlib.pyplot as plt
for module in pm, np, az:
print(f"{module.__name__} {module.__version__}")
@davipatti
davipatti / error.txt
Created November 8, 2022 20:43
[another aesara error]
ERROR (aesara.graph.rewriting.basic): Rewrite failure due to: save_mem_new_scan
ERROR (aesara.graph.rewriting.basic): node: for{cpu,scan_fn}(Shape_i{0}.0, Elemwise{Cast{int8}}.0, IncSubtensor{Set;:int64:}.0)
ERROR (aesara.graph.rewriting.basic): TRACEBACK:
ERROR (aesara.graph.rewriting.basic): Traceback (most recent call last):
File "/Users/pattinson/.virtualenvs/mfsera-env/lib/python3.10/site-packages/aesara/graph/rewriting/basic.py", line 1933, in process_node
replacements = node_rewriter.transform(fgraph, node)
File "/Users/pattinson/.virtualenvs/mfsera-env/lib/python3.10/site-packages/aesara/graph/rewriting/basic.py", line 1092, in transform
return self.fn(fgraph, node)
File "/Users/pattinson/.virtualenvs/mfsera-env/lib/python3.10/site-packages/aesara/scan/rewriting.py", line 1445, in save_mem_new_scan
nw_input = expand_empty(_nw_input, tmp_idx)
@davipatti
davipatti / aesara-error.txt
Created October 20, 2022 19:18
[aesara error] Error when trying to use aesara.scan with aesara==2.8.7 pymc==4.2.2 #aesara #pymc
---------------------------------------------------------------------------
NotImplementedError Traceback (most recent call last)
File ~/.virtualenvs/mfsera-env/lib/python3.10/site-packages/aesara/link/basic.py:111, in Container.__set__(self, value)
109 try:
110 # Use in-place filtering when/if possible
--> 111 self.storage[0] = self.type.filter_inplace(
112 value, self.storage[0], **kwargs
113 )
114 except NotImplementedError:
@davipatti
davipatti / DefaultKeyBinding.dict
Last active November 16, 2022 16:24
[DefaultKeyBinding.dict] ~/Library/KeyBindings/DefaultKeyBinding.dict. Getting 'home' and 'end' keys to do what they're meant to on a mac. From https://damieng.com/blog/2015/04/24/make-home-end-keys-behave-like-windows-on-mac-os-x.
{
"\UF729" = "moveToBeginningOfLine:"; /* Home */
"\UF72B" = "moveToEndOfLine:"; /* End */
"$\UF729" = "moveToBeginningOfLineAndModifySelection:"; /* Shift + Home */
"$\UF72B" = "moveToEndOfLineAndModifySelection:"; /* Shift + End */
"^\UF729" = "moveToBeginningOfDocument:"; /* Ctrl + Home */
"@\UF729" = "moveToBeginningOfDocument:"; /* Command + Home */
"^\UF72B" = "moveToEndOfDocument:"; /* Ctrl + End */
"@\UF72B" = "moveToEndOfDocument:"; /* Command + End */
"$^\UF729" = "moveToBeginningOfDocumentAndModifySelection:"; /* Shift + Ctrl + Home */
#!usr/bin/env python3
def attendance(price: int):
"""Number of people that will attend for a given price."""
formula = 120 - ((price - 5) / 0.1) * 15
return max((formula, 0))
def cost_of_performance(n_people: int) -> float:
# Fish function for converting jupyter notebooks to html with and without code cells
function jnc
for file in $argv
set root (basename "$file" .ipynb)
jupyter nbconvert --to html "$file"
jupyter nbconvert --to html --TemplateExporter.exclude_input=True "$file" --output "$root"-nocode.html
end
end
@davipatti
davipatti / reactivity-adjust-nth-optim.R
Created August 5, 2022 14:39
Run optimizeAgReactivity on every nth optimisation
library(Racmacs)
library(progress)
# Make a map with 1000 optimisations
n_optim <- 1000
titers <- read.titerTable("hi.csv")
map <- make.acmap(titer_table = titers, number_of_dimensions = 2, number_of_optimizations = n_optim)
#!/usr/bin/env python3
import time
import random
from typing import Generator, Iterable
def extend_loop(boxes: tuple[int, ...], loop: list[int], pointer: int) -> tuple[int, ...]:
"""