jsks

(defmacro defn-memoise [name args & body]
  (with-syms [$cache $fun]
    ~(def ,name (let [,$cache @{}
                      ,$fun (fn ,name ,args ,;body)]
                  (fn [& lst]
                    (if (nil? (,$cache lst))
                      (set (,$cache lst) (,$fun ;lst))
                      (,$cache lst)))))))

jsks

#!/usr/bin/env -S julia -t auto
#
# Metropolis-within-Gibbs sampler for a poisson regression
#
#   y_i ~ poisson(λ_i)
#   λ_i = exp(X_i' β)
#   β_j ~ normal(0, 5)
#
###

jsks

# Analytical dirichlet model with flat priors --- see BDA3 ch23.
###

using Distributions, Plots, StatsBase

"""
    bayes_hist(x)

Bayesian parametric histogram using the Freedman-Diaconis rule for
number of bins.

jsks

using Distributions, LinearAlgebra, LogExpFunctions

function gaussians(μ, Σ)
    @views [MvNormal(μ[:, i], Σ[:, :, i]) for i in 1:size(μ, 2)]
end


function E_step!(Y, γ, π, μ, Σ)
    dists = gaussians(μ, Σ)
    lpi = log.(π)

jsks

#!/usr/bin/env Rscript
#
# Gaussian mixture model fit with the expectation-maximization algorithm
###

library(mclust)
library(mvtnorm)

set.seed(4323)

jsks

#!/usr/bin/env Rscript

library(cmdstanr)
library(dplyr)

options(mc.cores = parallel::detectCores() - 1)

load("./day5.rda")

# Actor list

jsks

data/
stan/*
!stan/*.stan

.DS_Store
.Rhistory
.Rprof.user

\#*\#
.\#*

jsks

# Makefile template for Rmarkdown projects
#
# Includes implicit rules to generate pdf and html versions of Rmd
# files as well as several commands for development workflows.
#
# Invoke `make help` to get started.
#
# By default, Make will compile the pdf versions of all Rmd files
# found in the current working directory.
#

jsks

/*
 * Copyright (c) 2020, Joshua Krusell
 *
 * Permission is hereby granted, free of charge, to any person obtaining
 * a copy of this software and associated documentation files (the
 * "Software"), to deal in the Software without restriction, including
 * without limitation the rights to use, copy, modify, merge, publish,
 * distribute, sublicense, and/or sell copies of the Software, and to
 * permit persons to whom the Software is furnished to do so, subject to
 * the following conditions:

jsks

From b7b0801c8b9a97a67b1e329d25fe090fe61a818d Mon Sep 17 00:00:00 2001
From: Jamie Magee <[email protected]>
Date: Thu, 31 Oct 2019 09:57:34 +0100
Subject: [PATCH] Python 3.8 support (#10777)

* Python 3.8 support

time.clock() was deprecated in Python 3.3, and removed in Python 3.8. Python 3.8 will be released on 14th October. See [PEP 596](https://www.python.org/dev/peps/pep-0569/) for more info.

Tested locally against Python 3.8.0rc1.