Skip to content

Instantly share code, notes, and snippets.

View StaffanBetner's full-sized avatar

Staffan Betnér StaffanBetner

23 followers · 89 following
View GitHub Profile
@StaffanBetner
StaffanBetner / iterfuns.hpp
Last active March 8, 2024 20:56
Sample Bootstrap Weights within Stan (once for every iteration)
#ifndef DIRICHLET_RNG_WRAPPER_HPP
#define DIRICHLET_RNG_WRAPPER_HPP
#include <stan/math.hpp>
#include <boost/random/mersenne_twister.hpp>
#include <chrono>
#include <Eigen/Dense>
#include <iostream>
// Declare an integer to keep track of the iteration count
@StaffanBetner
StaffanBetner / 2023-10-04 Bortom linjäritet - en introduktion till splines.html
Created September 28, 2023 13:01
This file has been truncated, but you can view the full file.
<!DOCTYPE html>
<html lang="en"><head>
<script>/*!
* clipboard.js v2.0.11
* https://clipboardjs.com/
*
* Licensed MIT © Zeno Rocha
*/
!function(t,e){"object"==typeof exports&&"object"==typeof module?module.exports=e():"function"==typeof define&&define.amd?define([],e):"object"==typeof exports?exports.ClipboardJS=e():t.ClipboardJS=e()}(this,function(){return n={686:function(t,e,n){"use strict";n.d(e,{default:function(){return b}});var e=n(279),i=n.n(e),e=n(370),u=n.n(e),e=n(817),r=n.n(e);function c(t){try{return document.execCommand(t)}catch(t){return}}var a=function(t){t=r()(t);return c("cut"),t};function o(t,e){var n,o,t=(n=t,o="rtl"===document.documentElement.getAttribute("dir"),(t=document.createElement("textarea")).style.fontSize="12pt",t.style.border="0",t.style.padding="0",t.style.margin="0",t.style.position="absolute",t.style[o?"right":"left"]="-9999px",o=window.pageYOffset||document.documentElement.scrollTop,t.style.top="".concat(o,"px"),t.setAttribute("readonly",""),t.value=n,t);return e.contain
@StaffanBetner
StaffanBetner / custom pak functions.R
Last active September 28, 2023 12:12
if(!("pak" %in% as.data.frame(installed.packages())$Package)){install.packages("pak")}
require(pak)
pkg_load <- function(...) {
if(!("pak" %in% as.data.frame(installed.packages())$Package)){install.packages("pak")}
pack <- sapply(match.call(), as.character)[-1]
pack[!(pack %in% as.data.frame(installed.packages())$Package)] -> to_install
if(length(to_install) > 0) pak(to_install) else ""
lapply(X = pack, FUN = function(x){tryCatch(require(x, character.only = T), error = function(e){NULL})}) -> X
@StaffanBetner
StaffanBetner / qs caching assignment pipe.R
Last active February 16, 2022 14:49
library(qs)
`%q>%` <- function(lhs, rhs){as.character(substitute(rhs)) -> name
qcache(substitute(lhs), name = name) -> value
assign(x = name, value = value, envir = .GlobalEnv)}
@StaffanBetner
StaffanBetner / deso_valdeltagande_kommunval2018.csv
Created February 8, 2022 16:27
We can't make this file beautiful and searchable because it's too large.
county_code,county_name,county_name_long,municipality_code,municipality_name,municipality_name_long,regso_code,regso_name,deso,valdeltagande
01,Stockholm,Stockholms län,0114,Upplands Väsby,Upplands Väsby kommun,0114R001,Bollstanäs,0114C1010,0.88
01,Stockholm,Stockholms län,0114,Upplands Väsby,Upplands Väsby kommun,0114R001,Bollstanäs,0114C1030,0.88
01,Stockholm,Stockholms län,0114,Upplands Väsby,Upplands Väsby kommun,0114R001,Bollstanäs,0114C1050,0.91
01,Stockholm,Stockholms län,0114,Upplands Väsby,Upplands Väsby kommun,0114R002,Hammarby-Vaxmyra,0114C1180,0.84
01,Stockholm,Stockholms län,0114,Upplands Väsby,Upplands Väsby kommun,0114R002,Hammarby-Vaxmyra,0114C1220,0.72
01,Stockholm,Stockholms län,0114,Upplands Väsby,Upplands Väsby kommun,0114R003,Odenslunda norra-Frestaby-Ekeby,0114C1100,0.83
01,Stockholm,Stockholms län,0114,Upplands Väsby,Upplands Väsby kommun,0114R003,Odenslunda norra-Frestaby-Ekeby,0114C1140,0.9
01,Stockholm,Stockholms län,0114,Upplands Väsby,Upplands Väsby kommun,0114R004,Odenslunda södra
@StaffanBetner
StaffanBetner / to_clipboard.R
Created November 12, 2021 10:57
to_clipboard <- function(x, row.names = F){
write.table(x, "clipboard", row.names=row.names, sep='\t', na = "")
}
@StaffanBetner
StaffanBetner / ordered_beta_intercept_notebook.html
Created October 7, 2021 12:24
This file has been truncated, but you can view the full file.
<!DOCTYPE html>
<html>
<head>
<meta charset="utf-8" />
<meta name="generator" content="pandoc" />
<meta http-equiv="X-UA-Compatible" content="IE=EDGE" />
@StaffanBetner
StaffanBetner / induced_dirichlet_ordbeta.R
Last active October 7, 2021 12:04
library(tidyverse)
library(magrittr)
library(brms)
#library(betareg) # for example data
#### Code declaring ordered beta family for brms ####
# code from Robert Kubinec
# custom family
ord_beta_reg <- custom_family("ord_beta_reg",
@StaffanBetner
StaffanBetner / names_born_sweden_2020.csv
Last active August 23, 2021 14:46
We can't make this file beautiful and searchable because it's too large.
# Data source: Statistics Sweden
# No data on county and municipality level if count < 3
# No data on country level if count < 2
"sex","county_code","county","municipality","municipality_code","name","count"
"male",NA,NA,NA,NA,"Aadam",2
"male",NA,NA,NA,NA,"Aadhiran",2
"female",NA,NA,NA,NA,"Aadya",3
"male",NA,NA,NA,NA,"Aahil",2
"female",NA,NA,NA,NA,"Aaliya",2
"female","12","Skåne län","Malmö","1280","Aaliyah",3
@StaffanBetner
StaffanBetner / mgcv ocat dharma.R
Created June 23, 2021 19:53
model.check.ocat <- function(model) {
# Some R 4.1 syntax here!
simulatedResponse = model %>%
predict(type="response") %>%
apply(1, \(x) sample(1:length(x), prob = x, size=1000, replace = TRUE)) %>%
t()
createDHARMa(
simulatedResponse = simulatedResponse,
observedResponse = model$y,
fittedPredictedResponse = model %>% predict(), # the linear predictor