{{ message }}

Instantly share code, notes, and snippets.

# Nicholas J. Matzkenmatzke

Last active May 27, 2019
View prt_alone_v1.R
 require("ape") # print tree in hierarchical format ####################################################### # prt ####################################################### #' Print tree in table format #' #' Learning and using APE's tree structure can be difficult and confusing because much of the information is #' implicit. This function prints the entire
Created Apr 18, 2019
Revising Chris Rackauckas's "exp_euler_test.jl" from Julia 0.6 so it works in Julia 1.1
View _ex_Eulers_WORKS.jl
 ######################################################## # Revising Chris Rackauckas's "exp_euler_test.jl" from # Julia 0.6 so it works in Julia 1.1 ######################################################## # # Original 0.6 code from: # # ChrisRackauckas/exp_euler_test.jl # https://gist.github.com/ChrisRackauckas/cf91f4575b0587a45c72620384612b82
Created Apr 22, 2019
BDQ
View BDQ.jl
 import Pkg using Pkg # for Pkg.add, Pkg.PackageSpec using DifferentialEquations using Plots using BenchmarkTools # for @benchmark using LSODA # for lsoda() #Pkg.resolve() using Sundials # for CVODE_BDF() e.g.
Created Apr 22, 2019
BDQ v5 w help from Chris
View BDQ_v5.jl
 using DifferentialEquations using BenchmarkTools # for @benchmark using LSODA # for lsoda() using Sundials # for CVODE_BDF() # State-dependent birth-death process: # Problem setup n = 500 # Number of states two = 2.0 # Constant (but missing in some published equations)
Last active Apr 24, 2019
time results
View BDQ_time_results.jl
 # Time results posted to: # https://gist.github.com/nmatzke/34ba6b3f8b548af92724bcb3f41bbf3d # # ODE problem setup here: # https://gist.github.com/nmatzke/b6332845747d7452b6e3b45564f460e8 @benchmark sol_CVODE = solve(prob, CVODE_BDF(linear_solver=:GMRES), dense=false, save_everystep=false, save_end=true, save_start=true) # memory estimate: 2.72 GiB # allocs estimate: 5711652 # --------------
Created Feb 27, 2020
Takes 4.6 seconds for 2000 exp operations
View exp_example.jl
 using LinearAlgebra tmpstr = "Any[2 9 0.08; 2 10 0.15; 2 11 0.22; 2 12 0.29; 2 13 0.36; 2 14 0.43; 3 9 0.02; 3 15 0.16; 3 16 0.23; 3 17 0.3; 3 18 0.37; 3 19 0.44; 4 10 0.03; 4 15 0.1; 4 20 0.24; 4 21 0.31; 4 22 0.38; 4 23 0.45; 5 11 0.04; 5 16 0.11; 5 20 0.18; 5 24 0.32; 5 25 0.39; 5 26 0.46; 6 12 0.05; 6 17 0.12; 6 21 0.19; 6 24 0.26; 6 27 0.4; 6 28 0.47; 7 13 0.06; 7 18 0.13; 7 22 0.2; 7 25 0.27; 7 27 0.34; 7 29 0.48; 8 14 0.07; 8 19 0.14; 8 23 0.21; 8 26 0.28; 8 28 0.35; 8 29 0.42; 9 30 0.31; 9 31 0.45; 9 32 0.59; 9 33 0.73; 9 34 0.87; 10 30 0.18; 10 35 0.45999999999999996; 10 36 0.6; 10 37 0.74; 10 38 0.88; 11 31 0.19; 11 35 0.32999999999999996; 11 39 0.61; 11 40 0.75; 11 41 0.89; 12 32 0.2; 12 36 0.33999999999999997; 12 39 0.48; 12 42 0.76; 12 43 0.8999999999999999; 13 33 0.21000000000000002; 13 37 0.35; 13 40 0.49; 13 42 0.63; 13 44 0.9099999999999999; 14 34 0.22000000000000003; 14 38 0.36; 14 41 0.5; 14 43 0.6399999999999999; 14 44 0.78; 15 30 0.05; 15 45 0.47; 15 46 0.61; 15 47 0.75; 15 48 0.89; 16 3
Last active Mar 30, 2020
BIOSCI220: Lab 4 exercise and questions
View lab04_exercise_v5.R
 # Module 2: "Model-based Inference" # https://canvas.auckland.ac.nz/courses/48087/pages/module-2-model-based-inference?module_item_id=747842 # You can find the Lab 4 Worksheet questions at: # # Canvas -> BIOSCI 220 -> Quizzes -> Lab 4 Worksheet # The below Lab 4 code is posted at: #
Last active Apr 6, 2020
BIOSCI220: Week 5 lab exercise
View week_05_exercise_v3.R
 ####################################################### # BIOSCI220: Week 5 Lab Exercise: Modelling epidemics ####################################################### ####################################################### # Setup ####################################################### #
Last active May 2, 2020
week_06_Discussion_example_option2.R
View week_06_Discussion_example_option2.R
 ####################################################### # DISCUSSION ASSIGNMENT: FIT R0s FOR A COUNTRY OF INTEREST # # For the final part of your assignment, you will # pretend you are an epidemiologist writing a report # on the spread of Covid-19 in a particular country, # and commenting on what you can tell about the effectiveness # of interventions at this point. #
Last active Apr 30, 2020
week_06_Discussion_example_option1_v2.R
View week_06_Discussion_example_option1_v2.R
 ####################################################### # DISCUSSION ASSIGNMENT: FIT R0s FOR A COUNTRY OF INTEREST # # For the final part of your assignment, you will # pretend you are an epidemiologist writing a report # on the spread of Covid-19 in a particular country, # and commenting on what you can tell about the effectiveness # of interventions at this point. #
You can’t perform that action at this time.